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PREFACE  TO  THE  FIRST  EDITION. 


OF  all  the  celestial  bodies  whose  motions  have  formed 
the  subject  of  the  investigations  of  astronomers,  the  Moon 
has  always  been  regarded  as  that  which  presents  the 
greatest  difficulties,  on  account  of  the  number  of  inequali- 
ties to  which  it  is  subject;  but  the  frequent  and  important 
applications  of  the  results  render  the  Lunar  Problem  one 
of  the  highest  interest,  and  we  find  that  it  has  occupied 
the  attention  of  the  most  celebrated  astronomers  from  the 
earliest  times. 

Newton's  discovery  of  Universal  Gravitation,  suggested, 
it  is  supposed,  by  a  rough  consideration  of  the  motions  of 
the  Moon,  led  him  naturally  to  examine  its  application  to 
a  more  severe  explanation  of  her  disturbances;  and  his 
Eleventh  Section  is  the  first  attempt  at  a  theoretical 
investigation  of  the  Lunar  inequalities.  The  results  he 
obtained  were  found  to  agree  very  nearly  with  those 
determined  by  observation,  and  afforded  a  remarkable 
confirmation  of  the  truth  of  his  great  principle ;  but 
the  geometrical  methods  which  he  had  adopted  seem 
inadequate  to  so  complicated  a  theory,  and  recourse  has 
been  had  to  analysis  for  a  complete  determination  of  the 
disturbances,  and  for  a  knowledge  of  the  true  orbit. 
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The  following  pages  will,  it  is  hoped,  form  a  proper 
introduction  to  more  recondite  works  on  the  subject:  the 
difficulties  which  a  person  entering  upon  this  study  is  most 
likely  to  stumble  at  have  been  dwelt  upon  at  considerable 
length,  and  though  different  methods  of  investigation  have 
been  employed  by  different  astronomers,  the  difficulties  met 
with  are  nearly  the  same,  and  the  principle  of  successive 
approximation  is  common  to  all.  In  the  present  work, 
the  approximation  is  carried  to  the  second  order  of  small 
quantities,  and  this,  though  far  from  giving  accurate  values, 
is  amply  sufficient  for  the  elucidation  of  the  method. 

The  differences  in  the  analytical  solutions  arise  from 
the  various  ways  in  which  the  position  of  the  moon  may 
be  indicated  by  altering  the  system  of  coordinates  to  which 
it  is  referred,  or  again,  in  the  same  system,  by  choosing 
different  quantities  for  independent  variables. 

D'Alembert  and  Clairaut  chose  for  coordinates  the  pro- 
jection of  the  radius  vector  on  the  plane  of  the  ecliptic  and 
the  longitude  of  this  projection.  To  form  the  differential 
equations,  the  true  longitude  was  taken  for  independent 
variable. 

To  determine  the  latitude,  they,  by  analogy  to  Newton's 
method,  employed  the  differential  variations  of  the  motion 
of  the  node  and  of  the  inclination  of  the  orbit. 

Laplace,  Damoiseau,  Plana,  and  also  Herschel  and  Airy, 
in  their  more  elementary  works,  have  found  it  more  con- 
venient to  express  the  variations  of  the  latitude  directly, 
by  an  equation  of  the  same  form  as  that  of  the  radius 
vector. 

Lubbock  and  Pontdcoulant,  taking  the  same  coordinates 
of  the  Moon's  position,  make  the  time  the  independent 
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variable ;  and  when  it  is  desired  to  carry  the  approximation 
to  a  high  order,  this  method  offers  the  advantage  of  not 
requiring  any  reversion  of  series. 

Poisson  proposed  the  method  used  in  the  planetary 
theory,  that  is,  to  determine  the  variation  in  the  elements 
of  the  Moon's  orbit,  and  thence  to  conclude  the  corre- 
sponding variations  of  the  radius  vector,  the  longitude, 
and  the  latitude. 

The  selection  of  the  method  followed  in  the  present  work, 
which  is  the  same  as  that  of  Airy,  Herschel,  &c.,  was  made 
on  account  of  its  simplicity;  moreover,  it  is  the  method 
which  has  obtained  sanction  in  this  University,  and  it  is 
hoped  that  it  may  prove  of  service  to  the  student  in  his 
reading  for  the  examination  for  Honours.  In  furtherance 
of  this  object,  one  of  the  chapters  (the  sixth)  contains  the 
physical  intrepretation  of  the  various  important  terms  in 
the  radius  vector,  latitude,  and  longitude.* 

The  seventh  chapter,  or  appendix,  contains  some  of  the 
most  interesting  results  in  the  terms  of  the  higher  orders, 
among  which  will  be  found  the  values  of  c  and  g  com- 
pletely obtained  to  the  third  order. 

The  last  chapter  is  a  brief  historical  sketch  of  the  Lunar 
Problem  before  Newton,  containing  an  account  of  the  dis- 
coveries of  the  several  inequalities  and  of  the  methods  by 
which  they  were  represented,  those  only  being  mentioned 
which,  as  the  theory  has  since  verified,  were  real  onward 
steps.  The  perusal  of  this  chapter  will  shew  to  what  extent 
we  are  indebted  to  our  great  philosopher;  at  the  same 
time  we  cannot  fail  being  impressed  with  reverence  for  the 

*  See  the  Report  of  the  " Board  of  Mathematical  Studies"  for  1850. 


viij  PREFACE. 

genius  and  preservance  of  the  men  who  preceded  him, 
and  whose  elaborate  and  multiplied  hypotheses  were  in 
some  measure  necessary  to  the  discovery  of  his  simple 
and  single  law. 

I  take  this  opportunity  of  acknowledging  my  obligations 
to  several  friends,  whose  valuable  suggestions  have  added 
to  the  utility  of  the  work. 

HUGH  GODFEAY. 

St.  John's  College,  Cambridge, 
April,  1853. 


PREFACE  TO  THE  TRIED  EDITION, 

In  the  present  edition  several  important  alterations  have 
been  made,  many  of  the  parts  have  been  entirely  re-written, 
and  a  general  re-arrangement  of  the  whole  subject  has 
taken  place.  In  Chap.  IV.  Sect.  3,  the  necessity  for  the 
computation  of  the  elements  in  the  order  s,  w,  9  is  shewn. 
In  the  Quarterly  Journal,  Mr.  Walton  pointed  out  the 
logical  necessity  for  having  the  value  of  u  completely  correct 
to  the  second  order  before  proceeding  to  find  6,  although 
the  result  is  the  same  as  when  the  preparation  for  both 
values  is  carried  on  simultaneously.  I  have  here  proved 
that  there  is  the  same  necessity  for  making  the  value  of  s 
precede  that  of  w,  although,  again,  the  change  has  no  effect 
on  the  result. 

The  term  of  the  third  order  "  Parallactic  Inequality  "  will 
be  found  completely  investigated  in  the  appendix ;  and  at 
the  end  of  the  volume  a  selection  of  questions  from  University 
and  College  Examination  Papers  has  been  introduced. 

August,  1871. 
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LUNAR   THEORY. 


CHAPTER  1. 

INTRODUCTION. 

BEFORE  proceeding  to  the  consideration  of  the  moon's 
motion,  it  will  be  desirable  to  say  a  few  words  on  the  law 
of  attractions,  and  on  the  peculiar  circumstances  which  enable 
us  to  simplify  the  present  investigation. 

1.  The  law  of  universal  gravitation,  as  laid  down  by 
Newton,  is  that  u  Every  particle  in  the  universe  attracts  every 
other  particle,  with  a  force  varying  directly  as  the  mass  of  the 
attracting  particle  and  inversely  as  the  square  of  the  distance 
between  them" 

The  truth  of  this  law  cannot  be  established  by  abstract 
reasoning ;  but  as  it  is  found  that  the  motions  of  the  heavenly 
bodies,  calculated  on  the  assumption  of  its  truth,  agree  more 
and  more  closely  with  the  observed  motions  the  more  strictly 
our  calculations  are  performed,  we  have  every  reason  to 
consider  the  law  as  an  established  truth,  and  to  attribute 
any  slight  discrepancy  between  the  results  of  calculation  and 
observation  to  instrumental  errors,  to  an  incomplete  analysis, 
or  to  our  ignorance  of  the  existence  of  some  of  the  dis- 
turbing causes. 

Of  the  last  source  of  error  there  is  a  remarkable  ex- 
ample connected  with  the  discovery  of  the  planet  Neptune, 
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which  became  known  to  us,  as  one  of  the  bodies  of  our 
system,*  solely  by  means  of  the  perturbations  it  produced  in 
the  calculated  places  of  the  planet  Uranus.  These  pertur- 
bations were  too  great  to  be  attributed  wholly  to  errors  of 
instruments  or  of  calculation ;  and  therefore,  either  the  law 
of  universal  gravitation  ceased  to  hold  for  a  body  so  remote 
as  Uranus,  or  else  some  unknown  cause  was  disturbing  the 
path  of  the  planet.  The  first  supposition  was  too  repugnant 
to  an  astronomer  of  the  nineteenth  century  to  be  entertained 
until  all  others  had  failed, — and  the  second  supposition  led 
to  the  detection  of  Neptune.  The  distinguished  names  of 
Adams  and  Le  Verrier  will  be  for  ever  connected  with  the 
history  of  this  planet,  and  their  solution  of  the  difficult  in- 
verse problem : — "  Given  the  perturbations  caused  by  an 
unknown  planet,  determine,  on  the  assumption  of  the  truth 
of  Newton's  law,  the  orbit  and  position  of  this  disturbing 
body," — will  always  be  considered,  one  of  the  triumphs  of 
Mathematical  Science. 

Another  remarkable  instance  of  the  vindication  of  Newton's 
law,  where  it  seemed  at  first  to  be  at  fault,  will  be  found 
mentioned  in  Chap.  VI.  in  connection  with  the  motion  of 
the  Lunar  Apogee.  The  error  was  owing  to  an  incomplete 
analysis;  but,  when  the  calculations  were  carried  out  more 
fully,  the  strongest  confirmation  of  the  law  was  afforded  by 
the  exact  agreement  of  the  results  with  observation. 

We  shall  proceed  to  apply  the  law  of  gravitation  to  in- 
vestigate the  circumstances  of  the  moon's  motion ;  and  then 
see  how  it  will  enable  us  to  assign  her  position  at  any  time 
when  sufficient  data  have  been  obtained  by  observation. 

2.  The  problem  would  be  one  of  extreme,  if  not  insur- 
mountable difficulty,  if  we  had  to  take  into  account  simul- 

*  It  had  been  seen  by  Dr.  Lament  at  Munich,  one  year  before  its  being  known 
to  be  a  planet.  "  Solar  System,  by  J.  R.  Hind." 
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taneonsly  the  actions  of  the  earth,  sun,  planets,  &c.,  on  the 
moon ;  but  fortunately  the  planets  are  so  small  or  so  distant 
that  their  action  may  be  neglected— at  any  rate  to  the  order 
of  our  present  approximation — and  the  attraction  of  the 
earth,  on  account  of  its  proximity,  is  very  much  greater 
than  the  disturbing*  action  of  the  sun  notwithstanding  his 
enormous  mass.  We  may  therefore  treat  the  question  as 
that  of  one  body,  the  moon,  revolving  about  another  body, 
the  earth,  and  continually  disturbed  by  a  third  body,  the 
sun.  This  is  the  celebrated  problem  of  Three  Bodies. 

Stated  in  this  general  form  its  exact  solution  has  hitherto 
defied  the  powers  of  the  analyst.  The  disturbing  action 
of  the  sun  can  be  expressed  without  difficulty  in  terms  of 
the  masses  and  distances,  but  the  integration  of  the  diffe- 
rential equations  cannot  be  effected.  If  however,  instead 
of  retaining  the  whole  disturbing  force  of  the  sun,  which, 
as  we  have  said,  is  small  compared  with  the  earth's  action, 
we  expand  its  expression  in  a  series  and  neglect  all  very 
small  terms,  it  is  then  found  possible  to  obtain  a  solution. 

This  breaking  up  of  the  expression  for  the  sun's  action  is 
obviously  equivalent  to  a  breaking  up  of  the  sun's  force ; 
and  if  we  afterwards  wish  to  take  into  account  the  very 
small  forces  which  have  been  thus  neglected,  or  the  small 
disturbances  due  to  the  action  of  the  planets,  the  following 
principle  will  shew  that  we  may  do  so  by  considering  their 
effects  separately, — and  that  the  algebraical  sum  of  all  the 
disturbances  so  obtained  will,  to  an  order  of  approximation 
far  beyond  that  which  we  contemplate,  be  the  same  as  the 
disturbance  due  to  their  simultaneous  action.  We  shall  find 

*  Since  the  sun  attracts  both  the  earth  and  moon,  it  is  clear  that  its  effects  on 
the  moon's  motion  relatively  to  the  earth,  or  the  disturbing  force,  will  not  be  the 
same  as  the  absolute  force  on  either  body.  The  absolute  force  of  the  sun  on 
the  moon  is  more  than  double  that  of  the  earth  (Art.  99),  but  the  relative  force 
does  not  exceed  T£y-  of  the  earth's  action  (Art.  25). 
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an  application  of  this  principle  in  the  investigation  of  the 
parallactic  inequality. 

Principle  of  Superposition  of  Small  Disturbances. 

3.  Let  a  particle  be  moving  under  the  action  of  any 
number  of  independent  forces, 
some  of  which  are  very  small, 
and  let  A  be  the  position  of  the 
particle  at  any  instant.  Let  two  of  these  small  forces  m^  mz 
be  omitted,  and  suppose  the  path  of  the  particle  under  the 
action  of  the  remaining  forces  to  be  AP  in  any  given  time. 

Let  APl  be  the  path  which  would  have  been  described 
in  the  same  time  if  ml  also  had  acted ;  APl  differing  very 
slightly  from  AP,  the  disturbance  being  PPa. 

Similarly,  if  m2  instead  of  m^  had  acted,  suppose  AP2  to 
represent  the  disturbed  path,  PP2  being  the  disturbance, 
(AP,  AP^  APZ  are  not  necessarily  in  the  same  plane,  nor 
even  plane  curves). 

Lastly,  let  A  Q  be  the  actual  path  of  the  body  when  both 
m1  and  w2  act.  Join  P^  Q. 

The  two  disturbances  PP2  and  P^Q,  being  due  to  the 
action  of  the  same  small  force  mz  on  the  slightly  different 
paths  AP  and  AP^  must  be  very  nearly  equal  both  in 
magnitude  and  in  direction.  What  difference  exists  between 
them  must  be  a  very  small  fraction  of  either  disturbance, 
and  may  be  neglected  compared  with  the  original  path. 
We  may,  in  fact,  expect  that  this  difference  which  is  the 
disturbance  due  to  ml  of  a  disturbance  due  to  mz  will  be 
of  an  order  compounded  of  the  orders  of  the  two.  Thus, 
if  the  disturbance  due  to  ml  and  mz  be  respectively  of  the 
second  and  fourth  orders,  the  difference  between  P^Q  and 
PP2  would  probably  be  of  the  sixth  order  compared  with 
A  P.  Therefore  PtQ  may  be  considered  parallel  and  equal 
to  PP.. 
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Hence  the  projection  of  the  whole  disturbance  PQ  on  any 
straight  line,  being  equal  to  the  algebraical  sum  of  the  pro- 
jections of  PP}  and  Pt  Qj  will  be  equal  to  the  algebraical  sum 
of  the  projections  of  the  separate  disturbances  PP^  PPZ. 

Next,  let  there  be  three  small  disturbing  forces  01,,  m^  m9. 
We  may  consider  the  joint  action  of  the  two  m^  m^  as  one 
small  disturbing  force  ;  therefore,  by  what  precedes,  the  total 
disturbance  along  any  axis  will  be  the  sum  of  the  separate 
disturbances  of  ml  and  of  the  system  mzJ  ms  ;  but  this  last  is 
the  sum  of  the  separate  disturbances  of  mz  and  ms  :  therefore 
the  whole  disturbance  equals  the  sum  of  the  three  separate 
disturbances. 

4.  This  reasoning  can  evidently  be  extended  to  any 
number  of  forces  ;  and  if  a?,  y,  z  be  the  coordinates  of  the 
disturbed  particle,  <f>  (#,  y,  z)  any  function  of  #,  y,  z  ;  the 
disturbance  produced  in  0  (a?,  y^  z)  will  be 

a?,  y,  z]  -  ^-  Bx  4  —  By  +  -^  Sz,  omitting  (&z)2,  &c., 


where  Sic=5iCj  +  Sa?2+...=  sum  of  disturbances  along  axis  of  a? 

due  to  separate  forces, 
fysBfyj+fy  +...»  ........................  along  axis  of  y, 

50  =  52,  +50s+...=  ........................  along  axis  of  z; 

therefore  S<£  (ar,  y,  0)  =  —  5^  +       %,  +  ^-  ^ 


=  5^  (SB,  3/,  0)  4  52^>  (or,  y,  2)  +  &c.7 

or  total   disturbance   equals   sum   of  separate   disturbances, 
which  establishes  the  principle. 

Since  <£  (a?,  #,  2)  may  be  the  radius  vector,  the  latitude 
or  the  longitude  of  the  disturbed  body,  it  follows  that  the 
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total  disturbance  in  any  of  these  elements  is  the  sum  of  the 
partial  disturbances. 

5.  But  we  must  still  prove  another  proposition,  without 
which  the  problem,  though  so  far  simplified,  would  scarcely 
be  less  complicated  than  in  its  most  general  form. 

Newton's  law  refers  to  particles,  whereas  the  sun,  earth, 
and  moon  are  large  nearly  spherical  bodies,  and  it  becomes 
necessary  to  examine  the  mutual  action  of  such  bodies. 
Now,  it  happens  that,  with  this  law  of  force,  the  attraction 
of  one  sphere  on  another  can  be  correctly  obtained,  and 
leaves  the  question  in  exactly  the  same  state  as  if  they  were 
particles.  (Princip.  lib.  I.  prop.  75.) 

Attractions  of  Spherical  Bodies. 

6.  Let  P  be  a  particle  situated  at  a  distance  OP—  a  from 
the  centre  of  a  sphere  of  uniform 

density  p  and  radius  c.  The 
particle  being  without  the  sphere, 
a>  c. 

Let  the  whole  sphere  be  divi 
ded  into  circular  laminae  by  planes 
perpendicular  to  OP.     Let  SQ  be 
one  of  these.     PS=x,  PQ  —  z^  and  thickness  of  lamina  =  $x. 

Next,  let  this  lamina  be   divided  into  concentric  rings. 
Let  RS  =  r  be  the  radius  of  one  of  these  rings  and  Sr  its 
breadth,  LRPO^B*, 
therefore  r  —  x  tan  6, 

The  attraction  of  an  element  R  of  this  ring  on  the  particle 

.,          mass  of  element    . 
P  will  be  -      —  „„-,—     -  along  PZ2,  and  the  resolved  part 


PR* 
of  this  along  PO  will  be 


mass  of  element 
"  XA  sec'0 


COS0. 


ATTRACTIONS. 


But  the  resultant  attraction  of  the  whole  ring  will  clearly 
be  the  sum  of  the  resolved  parts  along  PO  of  the  attractions- 
of  its  constituent  elements ;  therefore, 

c\  51        £ 

attraction  of  ring  =  — ~- — ^—  cos  6  =  %7rp  sin  0Sx .  80 ; 
x  sec  c/ 

r-i- 
sin  d.dO 


z 

Again,    s2  =  a;2  +  c2  -  (a  -  x)*  =  2ax  -  (a2  -  c2) ; 
therefore  zBz  —  aSx, 

and  attraction  of  lamina  =  2jrp  I  -  2 $z }  j 

/.  attraction  of  whole  sphere  =  %7rp  \ — ^  — z      \ 

(from  z  —  a  —  c  to  2  = 
4-Trpc3 


where  M—  mass  of  sphere  =  — |—  . 

o 

Hence,  the  attraction  of  the  whole  sphere  is  precisely  the 
same  as  if  the  whole  mass  were  condensed  into  its  centre. 

COE.  1.  The  attraction  of  a  shell  (radius  c  and  thickness 
Sc)  will  be  obtained  from  the  preceding  expression  by  diffe- 
rentiating it  with  respect  to  c,  and  is 

, ,        .         ,,  ,    „      47iy)c2. 8c      mass  of  shell 

attraction  of  shell  =  —        -  = , 

a  a2 

the  same  as  if  the  mass  were  collected  at  its  centre. 
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COR.  2.  Therefore,  the  attraction  of  a  heterogeneous  sphere 
en  an  external  particle  will  be  the  same  as  if  the  whole  mass 
were  condensed  into  its  centre,  provided  the  density  le  the  same 
at  all  points  equally  distant  from  the  centre,  for  then  the 
whole  sphere  may  be  considered  as  the  aggregate  of  an 
infinite  number  of  uniform  shells. 

7.  Let  us  now  consider  the  case  of  one  sphere  attract- 
ing another.     Suppose  P  in  the  preceding  article  to  be  an 
elementary  particle  of  a  sphere  M' ,  whose  centre  0'  suppose 
at  a  distance  a  from    0.     Then,  since  action  and  reaction 
are  equal  and  opposite,  P  will  attract  the  whole  sphere  M 
just  as  it  would  do  a  particle  of  mass  M  placed  at  0.     The 
same  is  true  of  all  the  elementary  particles  which  compose 
the  sphere  M\  therefore  the  sphere  M1  will  attract  the  sphere 
M  as  if  the  whole  mass  of  the  latter  were  condensed  into  its 
centre  0 ;  but  the  attraction  of  the  sphere  M'  on  a  particle 
0  is  the  same  as  if  the  attracting  sphere  were  condensed  into 
its  centre  0' ;  therefore, 

Two  spheres  attract  one  another  as  if  the  whole  matter  of 
each  sphere  were  collected  at  its  centre. 

8.  This  remarkable  result,  which,  as  may  be  shewn,  holds 
only  when  the  law  of  attraction  is  that  of  the  inverse  square 
of  the  distance,  or  that  of  the  direct  distance,  or  a  com- 
bination of  these    by   addition    or   subtraction,   reduces   the 
problem  of  the  sun,  earth,  and  moon  to  that  of  three  par- 
ticles.    The  slight  error  due  to  the  bodies  not  being  perfect 
spheres  is  here  neglected,  being   of  an   order  higher  than 
that  to   which  we  intend   to   carry  the   present  investiga- 
tion :  this  error,  however,  though  very  small,  is  appreciable, 
and  when  a  nearer  approximation  is  required,  it  becomes 
necessary  to  have  regard  to  this  circumstance.     (See  Ap- 
pendix, Art.  109). 


CHAPTER    II. 


MOTION   RELATIVE  TO  THE  EARTH. 

9.  When  a  number  of  particles  are  in  motion  under  their 
mutual  attractions  or  other  forces,  and  the  motion  relatively 
to  one  of  them  is  required,  we  must  bring  that  one  to  rest 
and  then  keep  it  at  rest,  without  altering  the  relative  motions 
of  the  others  with  respect  to  it. 

Now,  first,  the  chosen  particle  will  be  brought  to  rest 
by  giving  it  at  any  instant  a  velocity  equal  and  opposite  to 
that  which  it  has  at  that  instant ;  secondly,  it  will  be  kept 
at  rest  by  applying  to  it  accelerating  forces  equal  and 
opposite  to  those  which  act  upon  it. 

Therefore  give  the  same  velocity  and  apply  the  same 
accelerating  forces  to  all  the  bodies  of  the  system,  and  their 
absolute  motions  about  the  chosen  body,  which  is  now  at  rest, 
will  be  the  same  as  their  relative  motions  previously. 

Problem  of  Two  Bodies. 

10.  As  the  sun  disturbs  the  moon's  motion  with  respect 
to  the  earth,  it  is  important  to  know  what  that  motion  would 
have  been  if  no  disturbance  had  existed,  or  generally : — 

Two  bodies  attracting  one  another  with  forces  varying 
directly  as  the  mass  and  inversely  as  the  square  of  the  dis- 
tance^ to  determine  the  orbit  of  one  relatively  to  the  other. 

Let  M)  M'  be  the  masses  of  the  bodies,  M'  being  the  body 
whose  motion  relatively  to  M  is  required,  r  the  distance 
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between  them  at  any  time  2,  and  6  the  angle  between  r  and 
some  fixed  prime  radius. 

The  accelerating  force  of  Mon  M'  equals  -^  acting  towards 

M* 

M)  while  that  of  M'  on  M  equals  -r  in  the  opposite  direction. 

Therefore,  by  the  principle  above  stated,  we  must  apply  to 
both  M  and  M  '  accelerating  forces  equal  and  opposite  to  this 
latter  force,  and  M'  will  move  about  M  fixed,  the  accelerating 

force  on  M'  being  -  —  -^  —  =  /*M*,  if  p  =  M+M'  and  r  =  -  . 

rr  d*u  IJL 

Hence,  -^  +  u  =  p  , 

Ct(7  /* 

7/1 

where  h  =  r*  -=-  =  twice  the  area  described  in  a  unit  of  time, 
and  integrating 


e  and  a  being  constants  to  be  determined  by  the  circumstances 
of  the  motion  at  any  given  time. 

This   is   the  equation  to  a  conic   section   referred  to  its 

72 

focus,  the  eccentricity  being  e,  the  semi-latus  rectum  —  ,  and 

the  angle  made  by  the  apse  line  with  the  prime  radius  a. 

By  observations  made  on  the  actual  path  of  the  moon, 
which  cannot  differ  very  widely  from  its  undisturbed  orbit, 
we  infer  that  this  latter  would  be  an  ellipse  with  an 
eccentricity  about  ^. 

In  the  same  manner,  if  the  sun  and  earth  be  the  two 
bodies  considered,  the  relative  orbit  would  be  an  ellipse  with 
an  eccentricity  about  ^G. 

11.  The  angle  6  —  a.  between  the  radius  vector  and  the 
apse  line  is  called  the  true  anomaly. 
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If  n  is  the  angular  velocity  of  a  radius  vector  which, 
moving  uniformly,  would  accomplish  its  revolution  in  the 
same  time  as  the  true  one,  both  passing  through  the  apse 
at  the  same  instant;  then  nt  +  s-a.  is  called  the  mean 
anomaly  where  s  is  a  constant  depending  on  the  instant 
from  which  the  time  is  reckoned,  its  value  being  the  angle 
between  the  prime  radius  and  the  uniformly  revolving  one 
when  £  =  0. 

Thus,  let  Jt/T  be  the  fixed  line  or  prime  radius, 
A  the  apse, 
M'  the  moving  body 
at  time  t, 

Mfj,  the  uniformly  re- 
volving radius  at  same  time,  the 
direction  of  motion  being  repre- 
sented by  the  arrow. 

And,  let  MD  be  the  position  of  M/J,  when  t  =  0, 

then  TJf-Z)  =  s  and  is  called  the  epoch,* 


TMA  =  a  =  longitude  of  the  apse  ; 
therefore,  mean  anomaly  =  AMp  =  nt  -f  e  -  a, 
true  anomaly  =  AMM'  =  TMM  '  -  tMA  =  0  -  a. 

12.    To  express  the  mean  anomaly  in  terms  of  the  true  in 
a  series  ascending  according  to  the  powers  of  e,  as  far  as  e2. 
2-7T  2  area 

27T-7-  -  ?  - 


periodic  time 
h 


*  We  shall  see  (Art.  36)  that  the  introduction  of  the  epoch  e  is  avoided 
in  the  Lunar  Theory  by  a  particular  assumption ;  but  in  the  Planetary  it  forms 
one  of  the  important  elements  of  the  orbit. 
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therefore  h  =  na*  (1  —  e2)*, 

dt       r*      aVl-e2)"  I 


dti      h  h  {l+ecos(0-a)J2 


(1  -  f  e2)  {1  -  2e  cos  (0  -  a)  +  3e2  cos2  (0  -  a)} 


-  {1  -  2e  cos  (0  -  a)  +  }e"  cos2  (0  -  a)}  ; 


therefore     nt  +  z  =  0-2e  sin  (0  -  a)  -f  fe2  sin  2  (0  -  a), 
or     (nt  +  s  -  a)  =  (6  -  a)  -  2e  sin.  (0  -  a)  +  f  e2  sin  2  (0  -  a), 
the  required  relation. 

13.  To  express  the  true  anomaly  in  terms  of  the  mean  to 
the  same  order  of  approximation. 

B  -  a.  =  nt  +  e  -  a  +  2e  sin  (0  -  a)  -  |e2  sin  2  (0  -  a)  ,  .  .(1)  ; 
.*.  6  —  a  =  nt  +  s  —  a  first  approximation. 
Substituting  this  in  the  first  small  terms  of  (1),  we  get 
Q  —  a  =  nt  +  s  —  a  +  2e  sin  (nt  -I-  s  —  a)...  a  second  approximation. 

Substitute  the  second  approximation  in  that  small  term 
of  (1  )  which  is  multiplied  by  e,  and  the  first  approximation 
in  that  multiplied  by  e2,  the  result  will  be  correct  to  that 
term,  and  gives 

0  -  a.  =  nt  +  e  —  a  +  2e  sin  {nt  -f  e  —  a  +  2e  sin  (nt  +  e  —  a)} 

-  |e*sin2(w«  +  e  -a) 
=  nt  +  e  -  a  -f  2e  sin  (nt  +  £  —  a) 

4-4e2  cos(w$+e  —  a)  sin(wZ+£—  a)—  |e2  sin2(w2-f  e-  a) 
=  nt  +  s  —  a  -f  2e  sin  (w£  +  e  —  a)  +  f  e2  sin  2  (nt  +  z  —  a), 
the  required  relation. 
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The  development  could  be  carried  on  by  the  same  process 
to  any  power  of  e,  the  coefficient  of  e3  would  be  found 

•}|  sin  3  (nt  +  £  —  a)  -  J  sin  (nt  +  e  —  a), 
but  in  what  follows  we  shall  not  require  anything  beyond  e2. 

Problem  of  Three  Bodies. 

14.  In  order  to  fix  the  position  of  the  moon  with  respect 
to  the  centre  of  the  earth,  which,  by  means  of  the  process 
described  in  Art.  (9),  is  supposed  brought  to  rest,  we  must 
have  some  determinate  invariable  plane  passing  through  the 
earth's  centre  to  which  the  motion  may  be  referred. 

If  the  sun  and  earth  were  the  only  bodies  in  the  universe,, 
then  the  plane  of  the  elliptic  orbit  (in  which,  according 
to  the  last  section,  the  motion  would  take  place)  would  be 
fixed,  and  might  be  taken  as  the  plane  of  reference;  buty 
as  soon  as  we  take  into  account  the  action  of  the  moon 
and  planets — especially  the  moon — the  plane  ceases  to  be 
fixed,  and  some  other  plane  must  be  found  not  affected  by 
these  disturbances. 

Theory  teaches  us  that  such  a  plane  exists,*  but  as  its 
exact  determination  can  only  be  the  work  of  time,  the 


*  See  Poinsot,  "  Theorie  et  determination  de  V equateur  du  systeme  solaire" 
where  he  proves  that  an  invariable  plane  exists  for  the  solar  system,  that  is,  a 
plane  whose  position  relatively  to  the  fixed  stars  will  always  be  the  same, 
whatever  changes  the  orbits  of  the  planets  may  experience;  but  as  its  position 
depends  on  the  moments  of  inertia  of  the  sun,  planets,  and  satellites,  and  there- 
fore on  their  internal  conformation,  it  cannot  be  determined  a  priori,  and  ages 
must  elapse  before  observation  can  furnish  sufficient  data  for  doing  so  a, 
posteriori. 

This  result  Poinsot  obtains  on  the  supposition  that  the  solar  system  is  a 
free  system ;  but  it  is  possible,  as  he  furthermore  remarks,  nay  probable,  that  the 
stars  exert  some  action  upon  it,  it  follows  that  this  invariable  plane  may  itself  be 
variable,  though  the  change  must,  according  to  our  ideas  of  time  and  space, 
be  indefinitely  slow  and  small. 
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following  theorem  will  supply  us  with  a  plane  whose  motion 
is  extremely  slow,  and  which  may  for  a  very  long  period, 
and  to  a  degree  of  approximation  far  beyond  that  to  which 
we  shall  carry  our  investigations,  be  considered  as  fixed  and 
coinciding  with  its  position  at  present. 

In  what  follows  the  action  of  the  planets  is  too  small 
to  be  taken  into  account. 

15.  The  centre  of  gravity  of  the  earth  and  moon  describes 
relatively  to  the  sun  an  orbit  very  nearly  in  one  plane  and 
elliptic ;  the  square  of  the  ratio  of  the  distances  of  the  moon 
and  sun  from  the  earth  being  neglected.* 


Let  $,  E,  M  be  the  centres  of  the  sun,  earth,  and  moon, 
G  the  centre  of  gravity  of  the  last  two.  Now  the  motion 
of  G  is  the  same  as  if  the  whole  mass  E+M  were  collected 
there  and  acted  on  by  forces  equal  and  parallel  to  the  moving 
forces  which  act  on  E  and  M.  The  whole  force  on  G  is 
therefore  in  the  plane  SEM;  join  SGf. 

Let  L  SGM=  o),  and  let  m  be  the  sun's  absolute  force. 
Moving  force  on  E—  0'     in  ES, 


f  **•    m.M  .     ,rc* 

moving  force  on  M—        a  in  Mb. 


*  This  ratio  is  about  ,^y,  and,  as  we  shall  see  Art.  (21),  such  a  quantity 
we  shall  consider  as  of  the  2nd  order  of  small  quantities,  and  its  square  therefore 
of  the  4th  order.  Our  investigations  are  carried  to  the  2nd  order  only. 
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These   applied  to    G-  parallel  to    themselves   are   equi- 
valent to 


m'.E.GE 

m'.M.GM 

mflirpotion  G-M 

VvTT      flip 

sirifl 

SE* 
m'.E.SG 

m'.M.SG 

,  triangle 
of  forces. 

Vint 

SE* 

SM5 
E        M 

M+E 

whence        E.  GE=  M.  GM=  (M+  E)  ^j^ 

Therefore   substituting,    and    dividing    by    M  +  E,*  we    get 
accelerating  forces* 


m  - 


,  GM.GE  (  1  1 

--  (OF  - 


and 


Again, 


\  . 

m  directl 


directlon 


E         M 


*  +  2SG.GEcosco 


IUC1CJLUIO 

similarly 

SE3  '    SG3  \ 
1           1     / 

Q/7  '  WOW  1 

ZGM        \ 

omitting  quantities 
to  be  neglected. 

SM*  ~  SG'  ( 

/S  (jr                I  . 

'} 

*  In  strictness  it  would  be  necessary,  since  we  have  brought  S  to  rest  to 
apply  to  both  M  and  JE,  and  therefore  to  G,  accelerating  forces  equal  and  opposite 
to  those  which  E  and  M  themselves  exert  on  S ;  but  the  mass  of  S  is  so  large 
compared  with  those  of  E  and  M,  that  we  may  safely  neglect  these  forces  in 
this  approximate  determination  of  the  path  of  G,  the  error  being  of  a  still 

higher  order  than  that  introduced  by  the  neglect  of  (--" 
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Therefore  the  accelerating  force  in  the  direction  GM 
3m.GM.GE 


SG* 


COS  £0 


n.  GM     GE 
—  3  (accelerating  force  of  sun  on  Or)  -=-„  .  -^     coso> 


=  0,  according  to  standard  of  approximation  adopted. 
And  the  accelerating  force  in  GS 

m'.SG  (E  +  M     Z(E.GE-M.GM) 


(E  +  M 
\SG* 


l* 


~M+ESG*  SG* 

m' 
—  -CTFT*  to  «W  same  approximation. 

OCr 

Hence,  the  force  on  G  is  a  central  force  tending  to  8 
and  varying  inversely  as  the  square  of  the  distance  ;  there- 
fore the  orbit  of  G  about  S  is  very  approximately  an  ellipse 
with  S  in  the  focus  ;  and  the  plane  of  this  ellipse  is,  as  far 
as  our  investigations  are  concerned,  a  fixed  plane. 

This  fixed  plane  is  called  the  plane  of  the  ecliptic,  or 
simply  the  ecliptic. 

16.  We  shall  now  suppose  a  plane  through  the  earth's 
centre  parallel  to  the  ecliptic;    and  this,  when  the  earth  is 
brought   to  rest,   will  give  us  the  fixed  plane  of  reference 
we  require  (14)  ;   the  ecliptic   then    making  small  monthly 
oscillations  from  one  side  to  the  other  of  our  fixed  plane. 

17.  Since   G  describes  an  ellipse  relatively  to  the  sun, 
the  sun  will  describe  the  same  ellipse  relatively  to   G  ;  but, 
as  seen  from  the  earth,  the  orbit  of  the  sun  will  be  slightly 
different,  as  we  shall  now  shew. 

First:  The  sun  will  have  a  latitude,  that  is,  will  be  out 
of  the  plane  of  reference;  and  this  latitude  will  be  of  the 
same  name  (north  or  south)  as  that  of  the  moon,  and 
deducible  from  it  when  the  ratio  of  the  distances  of  the  two 
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bodies  from  the  earth,  and  that   the  masses  of  the  earth 


and  moon  are  known.  For  if  S'EM'  be  the  fixed  plane 
through  E,  and  /S",  #',  M'  the  projections  of  £,  (7,  M,  then 
SES'  is  the  sun's  latitude,  and  MEM'  is  the  moon's;  and 
the  two  bodies  are  obviously  both  on  the  same  side  of 
the  plane. 

.    ,       .   ....      88'       GG      EG.  sin  (moon's  lat.) 


M       EM 

T5-M-si 

Now,  it  is  known  that  M  is  about  ^  of  E, 

andjEflf  .........  ?£0  of.E'tf; 

.,       .  .    ,       ,   i  XN      sin  (moon's  lat.) 

therefore     sin  (sun  s  lat.)  =  —  -^-^-  -—  —    -nearly. 

And  as  the  moon's  latitude  never  exceeds  5°  9',  the  sun's 
latitude  will  always  be  less  than  1". 

Secondly,  as  to  the  sun's  longitude:  let  ET  be  the 
direction  of  the  first  point  of  Aries,  —  that  is,  a  fixed  line 
in  our  plane  of  reference  from  which  the  longitudes  of  the 
bodies  are  reckoned.  Then  TES'  —  6'  is  the  sun's  longitude 
as  seen  from  E;  but,  as  seen  from  6r,  the  sun's  longitude 
would  be  the  angle  between  (r>S  and  a  line  through  G  (not 
drawn  in  the  figure)  parallel  to  ET. 

The  difference  between  these  two  longitudes  is  the  angle 
and,  if  EG'S'  =  co, 

tril     EG'  .  sinw 

G  =7  smo)=-  approximately; 
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therefore  smES'G  never  exceeds  53}^)  and  ES'G'  is  a 
small  angle  not  exceeding  7". 

Also  ES'  ~8'G<  EG'  <  „  J^  S'  G'. 

18.  Now,  by  assuming  the  longitude  and  distance  of  the 
sun  as  seen  from  E  to  be  the  same  as  when  seen  from  G, 
we  commit  the  above  small  errors  in  the  position  of  $; 
that  is,  we  assume  the  sun  to  be  at  8"  instead  of  S,  8'  S" 
being  drawn  equal  and  parallel  to  G'E.  If  our  object  were 
the  determination  of  the  sun's  position,  it  would  be  neces- 
sary to  take  this  into  account  ;  but  when  investigating  the 
motion  of  the  moon,  the  disturbing  action  of  the  sun  will 
not,  on  account  of  his  great  distance,  be  appreciably  altered 
by  supposing  him  to  be  at  8"  instead  of  8',  and  the  errors 
thus  introduced  are  much  too  small  to  be  considered. 

Hence  we  may  assume  that  the  motion  of  the  sun  about 
the  earth  at  rest  is  an  ellipse  having  the  earth  for  its  focus, 
and  its  equation 

tt'  =  a'{l+e'cos(0'-?)}, 

where  f  is  the  longitude  of  the  sun's  perigee  and  e  =  ^  is 
the  eccentricity,  and  we  are  safe  that  no  appreciable  error 
will  ensue  in  the  determination  of  the  moon's  place.* 

*  That  is,  as  far  as  the  three  bodies  alone  are  concerned  ;—  but,  since  the 
attractions  of  the  planets  may,  and  in  fact  do,  disturb  the  elliptic  orbit  of 
the  sun  about  G,  the  same  cause  will  disturb  the  assumed  orbit  about  E.  A 
remarkable  result  of  this  disturbance  is  noticed  in  Appendix,  Art.  (108). 


CHAPTER    III. 

RIGOROUS   DIFFERENTIAL  EQUATIONS   OF   THE  MOON'S  MOTION 
AND  APPROXIMATE  EXPRESSIONS  OF  THE  FORCES. 

19.  The  earth  having  been  reduced  to  rest  by  the  process 
described  in  Art.  (9),  let  us  take  its 
centre  E  for  origin.  Let  r,  6  be  the 
polar  coordinates  of  the  projection  M' 
of  the  moon  on  the  fixed  plane  of  re- 
ference, 6  being  the  longitude  TEM'  measured  from  a  fixed 
line  ET  in  that  plane.  Let  s  be  the  tangent  of  the  moon's 
latitude  MEM\  so  that  MM'  =  rs. 

Next,  let  the  accelerating  forces  which  act  upon  the  moon 
be  resolved  into  these  three : 

P  parallel  to  the  projected  radius  Jf'^and  towards  the  earth, 
T  parallel  to  the  fixed  plane,  perpendicular  to  P,  and  in  the 

direction  of  6  increasing  ^ 
S  perpendicular  to  the  fixed  plane  and  towards  it. 

The  equations  of  motion  will  be 

1  d 


u-±fL=-S (iii). 

20.   These   three   equations  for   determining  the  moon's 
motion  take  the  time  t  for  independent  variable ;  but  it  will 
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be  more  convenient  in  the  following  process  to  consider  the 
longitude  as  such,  and  our  next  step  will  be  to  change  the 
independent  variable  from  t  to  6. 
From  (i)  we  get 

7/3      ^7     /        ,-7/D\  ^7/3 

rffl  d_f  *w\         sdVf 
r  dt  dt  \    dt)  dt  > 

therefore  (V2  ^Y=  h*  +  2  /VrW, 

\     at  /  J 

JTT 

=  Z?2  suppose,  whence  H  -j^=  Trz. 

7/3  7T  -I 

Therefore  ._  =  -^  =  ^w2    if  w  =  - . 

at       r  r 


an  important  equation  connecting   the   time  t  with  the  longi- 
tude e. 

.     dr      dr    dO         TTdu 
' 


<?r       d  I        du\         „  ,  d  td 

u 


Substitute  in  (ii), 

therefore  ffu*        H        +  H'u*  =  P, 


•jT]  f- 

and  substituting  for  H3  and  for  S-^  their  values  #*-!-  2  l—d0 

T 

and  -3  ,  then  dividing  by  AV,  and  transposing, 

d*u          ,  P         T   du 


this  is  called  the  differential  equation  of  the  moon's  radius 
vector. 
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Lastly, 


d9  __         f      s_ 

IT"      'dJ~efc         \dd     8d 


„    d2 
therefore    Ps-S=  H'u3 

therefore 


this  is  the  differential  equation  of  the  moon's  latitude. 

21.  If  these  three  equations  could  be  integrated  under 
these  general  forms,  the  problem  of  the  moon's  motion  would 
be  completely  solved  ;  for  as  only  four  variables  w,  0,  5,  and 
t  are  involved  (the  accelerating  forces  P,  T,  and  S  being 
functions  of  these),  the  values  of  three  of  them,  as  w,  0,  sr 
could  be  obtained  in  terms  of  the  fourth  t  ;  that  is,  the  radius 
vector,  longitude,  and  latitude  would  be  known  correspond- 
ing to  a  given  time. 

But  the  integration  has  never  yet  been  effected,  except 
for  particular  values  of  P,  T,  and  S  ;  and  the  method  which 
we  are  in  consequence  forced  to  adopt  is  that  of  successive 
approximation,  by  which  the  values  of  w,  0,  and  s  are 
obtained  in  series  in  which  the  terms  proceed  according  to 
ascending  powers  of  certain  small  fractions.  Some  one 
fraction  being  chosen  as  a  standard  with  which  all  other 
are  compared,  the  order  of  the  approximation  is  esteemed 
by  the  highest  power  of  the  small  fractions  retained. 
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It  is  usual  to  consider  -£-$  as  a  small  fraction  of  the  first  order, 
consequently  ^V  °f  ?V  =  ¥o^  ™  .......................  second  ...... 


and  so  on,  other  fractions  being  considered  as  of  the  1st, 
2nd,  &c.  orders,  according  as  they  more  nearly  coincide  with 
?V>  ^uj  &c- 

22.  It  is  necessary  therefore,  before  we  can  approximate 
at  all,  that  we  should  have  a  previous  knowledge  (a  rough 
one  is  sufficient)  of  the  values  of  some  of  the  quantities 
involved  in  our  investigations;  and  for  this  knowledge  we 
must  have  recourse  to  observation. 

We  shall  therefore  assume  as  data  the  following  results 
of  observation  : 

(1)  The  moon  moves  in  longitude  about  thirteen  times 
as  fast  as  the  sun.     The  ratio  of  the  mean  motions  in  longi- 
tude represented  by  m  is  therefore  about  y1^,  and  may  be 
considered  as  of  the  1st  order.* 

(2)  The  sun's  distance  from  the  earth  is  about  400  times 
as  great  as  the  moon's  distance. 

Hence  the  ratio  of  the  mean  distances  =  ^J^  is  of  the 
second  order.f 

(3)  The  eccentricity  e'   of  the  elliptic   orbit   which   the 
sun   approximately  describes  about  the  earth  is  about  ^ 
and  this,  approaching  nearer  in  value  to  ^  than  to   ?J-0, 
will  be  considered  as  of  the  1st  order. 

(4)  During  one  revolution,  the  moon  moves  pretty  ac- 
curately in  a  plane  inclined  to  the  plane  of  the  ecliptic  at 

*  This  approximate  value  of  m  is  easily  obtained  j  —  the  moon  is  found  to 
perform  the  tour  of  the  heavens,  returning  to  the  same  position  among  the 
fixed  stars,  in  about  27  J  days;  the  sun  takes  365£  days  to  accomplish  the 
same  journey. 

t  The  distances  of  the  luminaries  may  be  calculated  from  their  horizontal 
parallaxes,  found  by  observations  made  at  remote  geographical  stations.  (See 
the  Author's  Astronomy,  chap.  xvi.). 
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an  angle  whose  tangent  is  about  ^,  and  therefore  of  the 
1st  order.* 

(5)  Its  orbit  in  this  plane  is  very  nearly  an  ellipse  having 
the  centre  of  the  earth  in  its  focus,  and  whose  eccentricity  is 
about  equal  to  our  standard  of  small  fractions  of  the  1st  order, 
viz.  "^;  and  this  will  also  be  very  nearly  true  of  the  projec- 
tion of  the  orbit  on  the  plane  of  the  ecliptic.f 

To  calculate  the  values  of  the  forces. 

23.  We  are  now  in  possession  of  the  data  requisite  for 
beginning  our  approximations,  and  we  shall  proceed  to  the 
determination  of  the  values  of  P,  T,  and  S  in  terms  of  the 
coordinates  of  the  positions  of  the  sun  and  moon. 

Let  S,  E,  M  be  the  centres  of  the  sun,  earth,  and  moon, 


m,  E,  M  their  masses, 

E')  M',  the  projections  on  the  plane  of  the  ecliptic, 

6r  the  centre  of  gravity  of  E  and  M. 

*  That  the  moon's  orbit  during  one  revolution  is  very  nearly  a  plane  inclined 
as  we  have  stated,  will  be  found  by  noting  her  position  day  after  day  among  the 
fixed  stars;  and  the  sun's  path  having  previously  been  ascertained  in  a  similar 
way,  the  rules  of  Spherical  Trigonometry  will  easily  enable  us  to  verify  these 
statements.  (Astronomy,  p.  92). 

f  The  elliptic  nature  and  the  value  of  the  eccentricity  of  the  moon's  orbit  may 
be  found  by  daily  observation  of  her  parallax,  whence  her  distance  from  the  earth's 
centre  may  be  determined  :  corresponding  observations  of  her  place  in  the  heavens 
being  taken,  and  corrected  for  parallax  to  reduce  them  to  the  earth's  centre,  will 
determine  her  angular  motion.  Lines  proportional  to  the  distances  being  then 
drawn  from  a  point  in  the  proper  directions,  the  extremities  mark  out  the  form 
of  the  moon's  orbit. 

A  similar  method  applied  to  observations  of  the  diameter  of  the  sun  will 
determine  the  eccentricity  of  its  orbit.  (Astronomy,  p.  138). 
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The  forces  we  have  to  take  into  account  are,  the  forces 
which  act  directly  on  M,  and  forces  equal  and  opposite  to 
those  which  act  on  E' — these  last  being  applied  to  the  whole 
system  so  that  E  may  be  a  fixed  point. 


Attraction  of  E- 

attraction  of  8 
m     .     1TC 
-8M'mMS 
equivalent  to 


Forces  on  M. 
E 


mME 


Forces  on  E. 

Attraction  of  M~  -J^T. in  EM 

MJcj 


attraction  of  8  m'.EG 


equivalent  to 


Therefore,  the   whole    attraction   upon  M,   when   E  is 
brought  to  rest,  is 

E+M      ,MG-      EG\. 


and 


m'SG 


parallel  to  GS. 


These  are  the  rigorous  expressions  of  the  accelerating 
forces  on  M,  and  they  can  be  expressed  in  terms  of  the 
masses  and  coordinates  of  the  bodies. 

Approximate  values  of  P,  T,  S. 

24.   Our  investigations  will  be  carried  only  to  the  second 
order ;  it  will  be  sufficient  therefore  if,  in  the  preceding,  we 
neglect  small  quantities  of  the  fourth  and  higher  orders. 
Let  p  =  E+  M, 

MGM'  =  tan~*s  =  moon's  latitude, 

E'V  the  direction  of  the  first  point  of  Aries, 

SG  =  /  =  A ;  LVE'S=  &  =  longitude  of  sun, 

M' E'  =  r  =  i ;  LTE'M'  =6  =  longitude  of  moon, 

=  0  —  &  =  difference  of  longitude  of  sun  and  moon. 
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Now,  8M*=8M"+MM* 

=  SG*+  GMn  -  ZSG.GM'  cos  SGM'  +  MM'1 


<*l*       n@M>          an*,.       #^8\ 
=  r"  (1  -  2  —  r-  cos  S#Jf  +  —  jf-  J  ; 

therefore        -g--,  .  1,  |l  4  i_^  cos  ((9  -  0')  j  ; 

for  0-0'  or  £#'Jf'  differs   from    >S^Jf'  by   less   than  7', 


Art.  (17),  and  f — — J   is  neglected,  being  of  the  fourth  order. 

1            1     (     '    ^  f-J-  7<"  ^ 

Similarly,   ^  =  ^  j  1 ^-  cos  (6  -  6')  I ;      therefore 

the  accelerating  forces  on  the  moon  are  approximately 
u,         m' 


in  direction 


and        -(GM'+  GE')  cos((9-  6>')  .........  parallel  to  G8; 

whence  P=  (j^  4-  ^'  JfJ^  cozMGM'-  ^  M'E'  wf(0- 
u,  m'r      3m  r 


=-        M'E'  cos (0 -  ff)  sin (0 - 
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The  differential  equations  in  Art.  (20),  when  these  values 
of  the  forces  are  substituted  in  them,  would  contain  the 
variables  u  and  &  ;  but  u  is  given  in  terms  of  &  by  the 
equation  of  Art.  (18)  u  —a  \\  +e  cos(0'  —  £)},  and  we  shall 
find  means  to  establish  a  connexion  between  6  and  0',  which 
will  enable  us  to  eliminate  the  latter. 

25.  Before  proceeding  further,  it  will  be  important  to 
consider  the  order  of  the  disturbing  effect  of  the  sun's  action, 
compared  with  the  direct  action  of  the  earth.  If  we  examine 
the  values  of  P,  T,  and  $,  it  will  be  found  that  the  most 
important  of  the  terms  containing  ra',  which  are  clearly  the 
disturbing  forces  since  they  depend  upon  the  sun,  are  in- 

77?  /* 

volved  in  the  form  -73-  ,  while  those  independent  of  the  sun's 

action  enter  in  the  form  -5  . 
r 

We  must  therefore  find  the  order  of  —  ^  compared  with  -5  , 


c 

orof/° 


Now  the  orbits  being  nearly  circular,   and   m  the  ratio 
of  the  mean  motions,  Art.  (22),  we  have 

mean  motion  of  sun        periodic  time  of  moon 
~  mean  motion  of  moon  ~~   periodic  time  of  sun 


-m 


r 

therefore  -%  :  -5 : :  m2  :  1 , 


r      r 


or  the  disturbing  force  of  the   sun   is   of  the  second  order 
when  compared  with  the  direct  action  of  the  earth. 
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INTEGRATION  OF  THE  DIFFERENTIAL  EQUATIONS. 

SECTION  I. 

General  process  described. 

26.  If  in  the  differential  equations  we  substitute  the 
values  of  P,  1\  and  S  just  obtained,  we  are  led  to  the 
following  forms: 


~  +  s  = 


which  cannot  be  integrated  ;  but  it  is  found  possible  to  obtain 
a  solution  which  will  be  a  first  approximation,  if  we  retain 
only  the  important  terms  on  the  right-hand  side  of  each 
equation.  This  first  approximation  being  substituted  in  the 
equations,  the  terms  next  in  order  of  importance  may  be 
retained  and  a  new  solution  becomes  possible,  which  in  its 
turn  is  the  basis  of  a  further  approximation,  and  so  on. 

First,  neglect  the  disturbing  force  of  the  sun  which  is 
of  the  second  order,  and  also  the  moon's  latitude,  which, 
as  will  be  seen  by  referring  to  the  expressions  for  the 
forces  (Art.  24),  will  either  enter  to  the  second  power  or  else 
in  combination  with  the  disturbing  force. 
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When  this  is  done  the  equations  become  integrable,  and 
values  of  u  and  s  may  be  obtained  in  terms  of  6  correct  to 
the  first  order;  this  value  of  u  will  then  enable  us  to  get 
the  connexion  between  6  and  t  to  the  same  order. 

[Let  us,  however,  bear  in  mind  that  the  equations  thus  integrated  are  not 
the  differential  equations  of  the  moon's  motions,  but  only  approximate  forms 
of  them;  and  it  is,  therefore,  possible  that  the  results  obtained  may  be  ap- 
proximations to  the  true  solutions  only  for  a  limited  range  of  values. 

Whether  they  are  so  or  not,  will  appear  by  comparing  them  with  what  we 
already  know  of  the  motion  from  observation ;  and  this  previous  knowledge, 
in  the  event  of  their  not  being  approximations,  will  probably  suggest  such, 
modifications  of  them  as  will  render  them  so.] 

To  proceed  to  the  next  approximation: — We  must  sim- 
plify the  right-hand  members  of  our  general  equations  by 
employing  the  values  of  u  and  s  just  obtained,  retaining 
now  terms  of  the  second  order  which  we  had  previously 
neglected. 

The  integration  can  always  be  performed  when  the  right- 
hand  members  of  the  equations  are  circular  functions  (sines 
or  cosines)  of  6  or  of  its  multiples;  and  as  our  first  approxi- 
mation will  give  us  the  values  of  u  and  s  in  that  form,  the 
new  expressions  for  the  forces  will  also  be  of  the  same 
character,  and  the  equations  be  again  integrable. 

Thus  we  shall  obtain  new  values  of  u,  s,  t  correct  to  the 
second  order.  These  values,  introduced  in  the  same  manner 
in  the  second  members  and  terms  of  the  next  higher  order 
retained,  will  lead  to  a  third  approximation,  and  so  on,  to- 
any  order;  except  that  if  we  wish  to  carry  it  on  beyond 
the  third,  the  approximate  values  of  the  forces,  given  in 
Art.  (24),  would  no  longer  be  sufficiently  exact,  and  we 
must  obtain  more  correct  values  from  the  expressions  of 
Art.  (23). 

27.  We  must  here  stop  to  notice  a  peculiarity  in  these 
equations,  when  solved  by  this  process.  We  have  said  that 
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to  obtain  the  values  to  any  order,  all  terms  up  to  that  order 
must  be  retained  in  the  second  members  ;  but  it  may  happen 
that  a  term  of  an  order  beyond  that  to  which  we  are  work- 
ing would,  if  retained,  be  so  altered  by  the  integration  as 
to  come  within  the  proposed  order. 

Such  terms  must  therefore  not  be  rejected,  and  we 
shall  proceed  to  examine  by  what  means  they  may  be 
recognised. 

28.  Suppose  then  that  after  an  approximation  to  a  certain 
order,  the  substitutions  for  the  next  steps  have  brought  the 
equation  in  u  to  the  form 


where  the  coefficient  G  is  one  order  beyond  that  which  we 
intend  to  retain.  The  solution  of  this  equation  will  be  of 
the  form 


G1  being  a  constant  which  may  be  determined  by  substi- 
tuting this  value  of  u  in  the  differential  equation, 

whence  G'  =  -  -  -a  j 

from  which  we  learn  that  if  p  differs  very  little  from  1,  Gf 
will  be  at  least  one  order  lower  than  6?,  and  will  come 
within  our  proposed  approximation  ;  and  consequently  the- 
term  G  cos(p0+H)  must  be  retained  in  the  differential 
equation. 

The  equation  of  the  moon's  latitude  being  of  the  same 
form  a&  that  of  the  radius  vector,  the  same  remarks  apply 
to  it. 

29.  If  p  i»  small,  G'  =  G  very  nearly,  and  the  term? 
G  cos  (p6  +  H)  will  not  rise  in  importance  in  the  value  of  u. 
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We   should  therefore  be  led  to  reject  it  in  the    differential 
equation,  except  for  the  following  reason  : — 

In  finding  the  connexion  between  the  longitude  and  the 
time  (one  of  the  principal  objects  of  the  Theory),  we  must 
use  the  equation 


as     hu2V     "JAV 
Now,  substituting  for  u  its  value  in  terms  of  0,  we  shall 
find  that  such  a  term  as  G'  cos  (pO  -f  H)  in  u  gives  rise  to 
a  term   G"  cos(p6  +  H)  in  this  equation,  where   G  and  G' 
are  of  the  same  order,  then 

at        r 
de= ° 

G"  . 
whence  t=* —  SH 

P 

G" 

therefore,  when  p  is  a  small  quantity  of  the  first  order,  - 

will  be  one  order  lower  than  G",  and  the  term  will    have 
risen  in  importance  by  the  integration. 

T 
But  yet  further,  if  such  terms  occur  in  ^-3 ,  they  will  be 


twice  increased  in  value ;  for  they  increase  once  in  forming 

'  T 

£2—3  d0,  and  once  again,  as  above,  in  finding  t> 

30.    We  have,  therefore,  the  following  rule : — 
In  approximating  to  any  given  order  to  the  values  of  u 
and  s,  we  must,  in  their  differential  equations,  retain  periodical 
terms  ONE  ORDER  beyond  the  proposed  one,  when  the  coefficient 
of  d  in  their  argument*  is  near  unity. 

*The  angle  of  a  periodical   term  is  its  only   variable  part   and  is   called 
the  argument. 
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In  approximating  to  the  value  of  t,  we  must,  in  the  same 
equations,  also  retain  terms  ONE  ORDER  beyond  the  proposed 
one  when  the  coefficient  of  the  argument  is  near  zero  ;  and 

T 
when    such    terms  occur    in  ~-8-  they   must  be   retained  TWO 

ORDERS  beyond  the  proposed  approximation* 

31.  We  shall  here,  for  convenience  of  reference,  bring 
together  the  equations  and  the  approximate  expressions  for 
the  forces  :  — 

tfu  P  T      du  d2u  C  T 


Ps-S        T      ds  d*s       \    r  T 


p 


" 


A 


r,a-f*2)- 


-;;£  *•(•-») 


*  Instead  of  the  forces  which  really  act  on  the  moon,  we  originally  substituted 
three  equivalent  ones,   P,  T,  S ;  these  again  are,  by  the  preceding  expressions, 

p 
replaced  by  a  set  of  others.    For,  we  may  conceive  each  of  the  terms  in  j—  ,  Ac., 

to  correspond  to  a  force, — a  component  of  P,  T,  or  8;  each  force  having  the 
same  argument  as  the  term  to  which  it  corresponds,  and  therefore  going  through 
its  cycle  of  values  in  the  same  time.  Now,  by  Art.,  (29),  when  the  coefficient  of 
0  in  the  argument  is  near  unity,  the  term  becomes  important  in  the  radius  vector, 
and  when  near  zero,  in  the  longitude  :  hence,  a  force  whose  period  is  nearly  the 
same  as  that  of  the  moon,  produces  important  effects  in  the  radius  vector;  and  a 
force  whose  period  is  very  long  will  be  important  in  its  effect  on  the  longitude. 
See  Airy's  Tracts,  Planetary  Theory,  p.  78. 
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SECTION  II. 
To  solve  the  Equations  to  the  first,  order. 

32.  Neglect  all  terms  depending  on  the  disturbing  force, 
i.e.  those  which  contain  m\  such  terms  being  of  the  second 
order.  Art.  (25). 

The  latitude  s  of  the  moon  can  never  exceed  the  inclina- 
tion of  the  orbit  to  the  ecliptic  ;  but  this  inclination  is  of  the 
first  order,  Art.  (22),  therefore  s  is  at  least  of  the  first  order 
and  s2  may  be  neglected. 

„,,  P  T          Ps-S 

Whence          -,=,*;  ^  =  0;  —^-  =  0, 

and  the  differential  equations  become 


d2s 


whence      u  —  •—  (1  +  e  cos(0  -  a)}  ;  or,  writing  a  for  £  , 

w  =  a  {l  +  ecos(0-a)J  .........................  (Z7J, 

and  s  =  ksm(6  —  7)  ..................................  ($,), 

c,  a,  kj  7  being  the  four  constants  introduced  by-  integration. 

33.  These  results  are  in  perfect  agreement  with  what 
rough  observations  had  already  taught  us  concerning  the 
moon's  motion  Art.  (22);  for  u  —  a  \\  +e  cos(0-  a)}  repre- 
sents motion  in  an  ellipse  about  the  earth  as  focus. 

And  s  —  k  sin  (6  —  7)  indicates  motion  in  a  plane  inclined 
to  the  ecliptic  at  an  angle  tan"1/:. 

For,  if  TOM'  be  the  ecliptic, 
M  the  moon's  place,  MM'  an  arc 
perpendicular  to  the  ecliptic,  then 
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YJf'  =  0;   and  if  TO  be  taken  equal  to  7,  and   OM  joined 
by  an  are  of  great  circle,  we  have 

tan  MM'  =  tmMOM'  sin  OM'  ; 
or  s  =  ta,nMOM'  sin(0-  7), 

which,  compared  with  the  equation  above,  shews  that 


Therefore  the  moon  moves  in  a  plane  passing  through 
a  fixed  point  0  and  making  a  constant  angle  with  the 
ecliptic. 

The  values  of  a  and  7  introduced  in  the  above  solutions 
are  respectively  the  longtitude  of  the  apse  and  of  the  node. 

34.  What   the   equations   cannot   teach   us,  however,  is 
the  magnitude  of  the  quantities  e  and  k.     For  this  we  must 
have  recourse  to  observation,  and  by  referring  to  Art.  (22), 
we   see   that  e  is  about  -fa  and   k   about   j1^,  that   is,  both 
quantities  are  of  the  first  order.     Their  exact  values  cannot 
yet  be  obtained:    the   means   of  doing   so   from   multiplied 
observations  will  be  indicated  further  on. 

35.  Lastly,  to  find  the  connexion  between  t  and  0,  the 
equation  (C)  becomes,  T=0, 

dt_  _  J^        1  __  1  _ 
dB  ~  hit*  ~~  ha*  [l  +  e  cos(0  -  a)}'  ' 

ISow  this  is  the  very  same  equation  that  we  had  found 
connecting  t  and  6  in  the  problem  of  two  bodies.  Art.  (12),  as 
we  ought  to  expect,  since  we  have  neglected  the  sun's  action. 
Therefore,  if  p  be  the  moon's  mean  angular  velocity,  we 
should,  following  the  same  process  as  in  the  article  referred 
to,  arrive  at  the  result 


e  —  a+  ......  , 

which    is   correct    only    to   the  first   order,    since   we    have 

F 
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rejected  some  terms  of  the  second  order  by  neglecting  the 
disturbing  force. 

36.  The  arbitrary  constant  s,  introduced  in  the  process 
of  integration,  can  be  got  rid  of  by  a  proper  assumption  : 
this  assumption  is,  that  the  time  t  is  reckoned  from  the 
instant  when  the  mean  value  of  6  is  zero.* 

For,  since  the  mean  value  of  #,  found  by  rejecting  the 
periodical  terms,  is  pt  +  e  ;  if,  when  this  vanishes,  t  =  0,  we 
must  have  £  =  0  ;  therefore 


correct  to  the  first  order.f 

37.  We  have  now  obtained  three  results,  £/j,  $„  0l7 
as  solutions  to  the  first  order  of  our  differential  equations, 
and  we  must  employ  them  to  obtain  the  next  approximate 
solutions  ;  but,  before  they  can  be  so  employed,  they  will 
require  to  be  slightly  modified  —  in  such  a  manner,  however, 
as  not  to  interfere  with  their  degree  of  approximation. 

The  necessity  for  such  a  modification  will  appear  from 
the  following  considerations  : 

*  When  a  function  of  a  variable  contains  periodical  terms  which  go  through 
all  their  changes  positive  and  negative  as  the  variable  increases  continuously, 
.the  mean  value  of  the  function  is  the  part  which  is  independent  of  the  periodical 
terms. 

f  We  shall  also  employ  this  method  of  correcting  the  integral  in  our  next 
approximation  to  the  value  of  6  in  terms  of  t  ;  and  if  we  purposed  to  cany  our 
approximations  to  a  higher  order  than  the  second,  we  should  still  adopt  the 
same  value,  that  is,  zero,  for  the  arbitrary  constant  introduced  by  the  integration. 
To  shew  the  advantage  of  thus  correcting  with  respect  to  mean  values  :  suppose 
we  reckoned  the  time  from  some  definite  value  of  0,  for,  instance  when  0  —  0  ; 
then,  in  the  first  approximation, 

0  =  s  +  2e  sin  (e  —  a) 

is  the  equation  for  determining  the  constant  e,  and  in  the  second  approximation, 
«  would  be  found  from 

0  =  £+2esin  («-a)  +  $e2sin  2  («-«)+  ......  , 

giving  different  values  of  i  at  each  successive  approximation. 
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Suppose  we  proceed  with  the   values  already    obtained  ; 
we  have,  by  Art.  (24), 


2A 


=  a+ +  A  cos(0- 

and  this  being  substituted   in    the  differential  equation  (A) 
of  Art  (31;,  gives 

7*j 

—  +  u  =  a+ +  ^lcos(0-a)  + , 

the  solution  of  which  is 

w  =  a{l+ecos(0-a)}  + +  \Ad  sin(0  -a). 

Our  first  approximate  value  u  —  a  (1  +  e  cos(0  —  a)}  is  here 
corrected  by  a  term  which,  on  account  of  the  factor  0,  admits 
of  indefinite  increase,  and  thus  becomes  ultimately  more 
important  than  that  with  which  we  started.  Such  a  cor- 
rection is  inadmissible,  for  our  first  value  is  confirmed  by 
observation  (22;  to  be  very  nearly  the  true  one.  The  moon's 
distance,  as  determined  by  her  parallax,  is  never  much  less 
than  60  times  the  earth's  radius;  whereas  this  new  value 
of  w,  when  6  is  very  great,  would  make  the  distance  in- 
definitely small.  On  the  same  principle,  we  see  that  any 
solution,  which  comprises  a  term  of  the  form  Ad  sin (6  —  a), 
cannot  be  an  approximate  solution  except  for  a  small  range 
of  values  of  6. 

Such  terms,  '  if  they  really  had  an  existence  in  our  system,  must  end  in  its 
'  destruction,  or  at  least  in  the  total  subversion  of  its  present  state ;  but  when 
'they  do  occur,  they  have  their  origin,  not  in  the  nature  of  the  differential 
'  equations,  but  in  the  imperfection  of  our  analysis,  and  in  the  inadequate  repre- 
'sentation  of  the  perturbations,  and  are  to  be  got  rid  of,  or  rather  included  in 
'  more  general  expressions  of  a  periodical  nature,  by  a  more  refined  investigation 
<than  that  which  led  us  to  them.  The  nature  of  this  difficulty  will  be  easily 
'understood  from  the  following  reasoning.  Suppose  that  a  term,  such  as 
<«  sin  (AQ  +  B],  should  exist  in  the  value  of  M,  in  which  A  being  extremely 
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'  minute,  the  period  of  the  inequality  denoted  by  it  would  be  of  great  length  ; 
'  then,  whatever  might  be  the  value  of  the  coefficient  a,  the  inequality  would 
'  still  be  always  confined  within  certain  limits,  and  after  many  ages  would  return 
<  to  its  former  state. 

'Suppose  now  that  our  peculiar  mode  of   arriving  at  the  value  of  u  led  us 
'  to  this  term,  not  in  its  real  analytical  form  a  sin  (AQ  +  B),  but  by  the  way  of 
'  jts  development  in  powers  of  0,  a  +  /30  +  ^02  +  &c.  ;  and  that,  not  at  once,  but 
'  piecemeal,  as  it  were  ;   a  first   approximation   giving  us  only  the  term   a,   a 
'second  adding  the  term  /30,  and  soon.     If  we  stopped  here,  it  is  obvious  that 
'we  should  mistake  the  nature  of  this   inequality,  and  that  a  really  periodical 
'function,  from  the  effect  of  an  imperfect  approximation,  would   appear  under 
'the  form  of  one  not  periodical  .........     These  terms  in  the  value  of  u,  when  they 

'occur,  are  not  superfluous;  they  are  essential  to  its  expression,  but  they  lead 
'us  to  erroneous  conclusions  as  to  the  stability  of  our  system  and  the  general 
'  laws  of  its  perturbations,  unless  we  keep  in  view  that  they  are  only  parts  of 
'  series  ;  the  principal  parts,  it  is  true,  when  we  confine  ourselves  to  intervals 
'  of  moderate  Ipngth,  but  which  cease  to  be  so  after  the  lapse  of  very  long  times 
'the  rest  of  the  series  acquiring  ultimately  the  preponderance,  and  compensating 
'  the  want  of  periodicity  of  its  first  terms.-Sm  JOHN  HERSCHBL,  Encyclopaedia 
Metropolitans—  PHYSICAL  ASTRONOMY,  p.  679. 

38.  TO  extritcate  ourselves  from  this  difficulty,  and 
to  alter  the  solution  so  that  none  but  periodical  terms 
may  be  introduced,  let  us  again  observe  that  the  equation 

_  _  -f  w  =  —  =a,  which  gave   the    solution   £/J  and   thus   led 

to  the  difficulty,  is  only  an  approximate  form  of  the  first 
order  of  the  exact  equation.  Any  value  of  u,  therefore, 

which   satisfies  the  approximate  equation  -j^2  +  u  =  a  to  the 

first  order,  and  which  evades  the  difficulty  mentioned  above, 
may  be  taken  as  a  solution  to  the  same  order  of  the  exact 
equation. 

These  conditions  will  be  satisfied  if  we  assume 

u  =  a  {I  +  e  cos(c#  -  a)j  ; 


then  ^.,+u  =  a  +  ae(l-  c2)cos(c0-a) 

au 

=  a  to  the'first  order, 
provided  I  —  o*  be  of  the  first  order  at  least. 
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39.  The  introduction  of  the  factor  c  in  the  argument  is 
an  artifice  due  to  Clairaut ;  and  its  effect,  as  we  shall  see 
in  Art.  (66)  is  equivalent  to  supposing  that  the  apse  of  the 
moon's  orbit  is  not  fixed.  Now,  although  the  observations 
recorded  in  Art.  (22)  did  not  suggest  this,  we  must  bear  in 
mind  that  they  were  extremely  rough,  and  carried  on  only 
for  a  short  interval;  but  when  they  are  made  with  a  little 
more  accuracy,  and  extended  over  several  revolutions  of  the 
moon,  it  is  soon  found  that  both  her  apse  and  the  plane 
of  her  orbit  are  in  constant  motion. 

Clairaut  was  fully  aware  of  these  motions,  and  there  is 
no  doubt  that  he  was  led  to  the  above  form  of  the  value 
of  u  by  that  consideration,  and  by  his  acquaintance  with  the 
results  of  Newton's  ninth  section,  which,  when  translated  into 
analytical  language,  lead  at  once  to  the  same  form.* 

We  might,  therefore,  taking  for  granted  the  results  of 
observation,  have  commenced  our  approximation  at  this  step, 
and  have  at  once  written  down  u  =  a{l+e  cos(c#  —  a)} ;  but 
we  should,  in  so  doing,  have  merely  postponed  the  difficulty 
to  the  next  step,  since  there  again,  as  we  shall  find,  the 
differential  equation  is  of  the  form 

— —  -f  u  —  a  function  of  0, 

cL\j 

*  Newton  has  there  shewn,  that  if  the  angular  velocity  of  the  orbit  be  to 
that  of  the  body  as  G  -  F  to  G,  the  additional  centripetal  force  is  ^"/^  hzu*t 
the  original  force  being  juw2.  Therefore 


^1  -t  -F2*  -11 

d&      G*U~h*~  G*  h?  ~  G*  a' 


F 

where  -*  is  the  same  as  our  c. 
G 


u  =  a  j  1  -f  e  cos  f  —  0  —  a  J  j- , 
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the  correct  integral  of  which  would  be 


and  this  would  at  the  next  operation  bring  in  a  term  with 
6  for  a  coefficient,  which  we  now  know  must  not  be.  We 
shall,  therefore,  hereafter  omit  such  terms  as  A  cos(#  —  B) 
altogether,  and  merely  write 

u  =  a  [I  +e  cos(c#-  a)}  -f  ............. 

40.  The  value  of  c   will  be  found  more  and  more  cor- 
rectly at  each  successive  approximation,  by  always  assuming 
these  terms  as  the  first  two  terms  of  the  value  of  w,  then 
substituting  in  the  differential  equation  and  equating  coeffi- 
cients.    It  will  thus  be  found  that  ae(l  —  c*)  must  equal  the 
coefficient  of  cos(c#  —  a)  in  the  differential  equation  ;  and  this 
will  enable  us  to  determine  c  to  the  same  degree  of  approxi- 
mation   as    that    of    the    differential    equation    itself.      See 
Arts.  (52)  and  (103).     So  far,  all  that  we  know  about  c  is  that 
it  differs  from  unity  at  most  by  a  quantity  of  the  first  order. 

41.  In  carrying  on  the  solution  of  s,  the  same  difficulty 
arises  as  in  u^  and  it  will  be  found  necessary  to  change  it  into 

s  =  k  sin(  <7#  —  7), 

g  being  a  quantity  which  differs  from  unity  at  most  by  a 
quantity  of  the  first  order. 

The  introduction  of  g  is  connected  with  the  motion  of  the 
node  in  the  same  way  as  that  of  c  is  with  the  motion  of  the 
apse  (85)  ;  and  the  value  of  g  will  be  determined  by  a  process 
similar  to  that  explained  above  for  c.  Arts.  (51)  and  (102). 

42.  The  connexion  between  6  and  t  will  also  be  modified 
by  this  change  in  the  value  of  w, 

dt_       1  __  1 
d9~  Aail 


d^cl_  _  J. 1_ 

d»c6      hd*  {1+6  cos( 
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Here  cO  and  ct  hold  the  places  which  6  and  t  occupied 
in  (35)  ;  therefore 

c&  =  cpt  +  2e  sm(cpt  —  a), 
or  6  =  pt  +  2e  sin  (cpt  -  a), 

/> 

to  the  first  order,  since  -  =  e  to  the  first  order. 

43.  Since  the  disturbing  forces  are  to  be  taken  into 
account  in  the  next  approximation,  we  shall  have  to  use 
the  value  of  u  found  in  (18),  which  is 


but  this  introduces  6'  ;  we  must  therefore  further  modify  it 
by  substituting  for  &  its  value  in  terms  of  0,  and  it  will 
be  found  sufficient,  for  the  purpose  of  the  present  work, 
to  obtain  the  connexion  between  them  to  the  first  order, 
which  may  be  done  as  follows  : 

Let  in  be  the  ratio  of  the  mean  motions  of  the  sun  and 
moon, 

p'j     p  their  mean  angular  velocities;  .:p  =  mp, 
p't  +  fB,pt  ......  mean  longitudes   at  time  2,  /3  being  the  sun's 

longitude  when  £  =  0, 

6',      0  ......  true  longitudes  at  time  £, 

f,       a  ......  longitude  of  perigees  when  t  =  0  ; 

therefore  0'  —  f  =  sun's  true  anomaly, 

and  p'£+/3-f=  ......  mean  anomaly. 

But,  by  Art.  (13), 

true  anomaly  =  mean  anomaly  +  2e  sin  (mean  anomaly)  +  &c.  ; 
therefore  6'=  p't+  j3  +  2e  sin  (/  1  +  fi  -  f)  +  ...... 


to  the  first  order  ; 
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because      pt—Q  —  2e  sin  (c9  -  a)  to  the  first  order  by  (42). 
Whence      u  =  a  (1  +  e  cos  (m0  -f  /3  —  £)}  to  the  first  order. 

44.   The  values  of  sin  2  (0-  0')  and  cos  2  (0  -  0')  can  also 
be  readily  obtained  to  the  same  order  : 

sin  2  (6  -6') 

=  sin  {(2  -  2m)  6  -  20  -  4e'  sin  (w0  +  /3  -  f  ){ 

=  sin  {(2  -  2?n)  0  -  2/3}  -4e'sin(w0  +  £  -  f  )  cos{(2  -  2m)  0-2/3} 

=  sin  ((2  -  2m)  0  -  2/3}  -2<?'sm  {(2  -  m)  0  -  /9  -  £  } 

+  2e'sin{(2-3?7z)-3/3+£}. 
Similarly,  cos  2  (0  -  0')  =  cos  {(2  -  2m)  0  -  2/3f 


The  first  term  of  each  of  these  is  all  we  shall  require. 

SECTION  III. 

To  solve  the  equations  to  the  Second  Order. 

45*   Let  us  recapitulate  the  results  of  the  last  approxi- 
mation. 

u  =  a  {l+e  cos(c0-  a)}, 


6  -  &  -  (1  -  TO)  0  -  /3  -  2e  sin(w 

These  values  must  now  be  substituted  in  the  expressions 
for  the  forces,  retaining  terms  above  the  second  order,  when, 
according  to  the  criterion  of  Art.  (30),  they  promise  to  be- 
come of  the  second  order  after  integrating. 


SOLUTION  TO   THE   SECOND   ORDER.  41 

The  differential  equations  (A)  and  (B)  will  then  assume 
the  forms 


w 

^  +*=/(*), 

and  the  integration  of  these  will  enable  us  to  obtain  u  and  s 
to  the  second  order ;  after  which,  equation "( C)  -will  give  the 
connexion  between  0  and  t  to  the  same  order. 

/     '3 

46.   The  quantity  p-^- ,  which  we  shall  meet  with  as  a 

coefficient  of  the  terms  due  to  the  disturbing  force,  can  be 
replaced  by  mza^  m  being  the  ratio  of  the  mean  motions  of 
the  sun  and  moon. 

So  long  as  we  neglected  the  disturbing  force,  h  and  a 
had  determinate  values : — they  belonged  to  the  ellipse  which 
formed  our  first  imperfect  solution,  and  would  therefore  be 
known  from  the  circumstances  of  motion  in  that  ellipse  at 
any  instant ;  h  being  double  the  area  described  in  a  unit 
of  time,  and  a  the  reciprocal  of  the  semi-latus  rectum.  It 
would  consequently  be  impossible  to  assume  any  arbitrary 
connexion  between  them.  But,  when  we  proceed  to  a  second 
approximation  and  introduce  the  disturbing  force,  there  is 
no  longer  a  determinate  ellipse  to  which  the  h  and  a  apply : 
the  equation  /*  =  tfa  of  Art.  (32)  merely  shews  that  a  and  h 
must  refer  to  some  one  of  the  instantaneous  ellipses  which 
the  moon  could  describe  about  the  earth  if  the  disturbance 
were  to  cease,  and  we  are  at  liberty  to  select  any  one  of 
these  which  will  allow  us  to  proceed  with  our  approximation. 
The  particular  ellipse  is  determined  by  the  above  assumed 

relation    -p-^  =m*at   and   the   selection    is    suggested    and 

a 
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justified  by  the  following  reasoning: 

__p  _  average  period  of  moon  about  earth 
p  period  of  sun  about  earth 

but,  since  the  instantaneous  ellipses   are   nearly  circles,  we 
have,  as  in  Art.  (25), 

(period  of  moon  about  earth)2      ma3         , 
/ — :— i — ? —  — TTF-  =  — 3  nearly . 

(period  of  sun  about  earth)          /JLO, 

therefore  if  a  be  properly  chosen, 


fia 


AV  ' 


47.  If  we  examine  the  equations  of  Art.  (31)  it  will  be 
seen  that  we  cannot  solve  the  differential  equation  in  u  to 
the  second  order  until  we  know  the  value  of  s  accurately  to 
that  order;  for,  suppose 

s  —  k  sin(<7#  —  7)  +  F2, 

where  V2  contains  the  terms  of  the  second  order  not  yet 
found;  and  suppose  that  among  the  terms  in  F2  there  be 
one  Ls\n(l0  —  X),  then  this  value  of  s  substituted  in  the 
expression  (D)  gives 


=  a  {1  -  f*f  sin2(#0-  7)  -  3&Z  sin(#0  -  7)  sin(Z0-X).. 
and  equation  (A]  becomes,  so  far  as  this  term  is  concerned, 


Now,  if   it  should  happen    that  I  is  a  small   quantity,  then 
the  two  coefficients  g  -  1  and  g  -f  I  will  be  both  near  unity, 
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and  therefore  these  two  terms,  although  of  the  third  order 
in  the  differential  equation,  would  become  important  in  the 
value  of  u.  Again,  if  I  is  near  2,  then  I  -  g  will  be  near 
unity,  and  the  same  remark  will  apply. 

Such  terms  do  not  really  occur  in  the  value  of  5,  but 
there  is  none  the  less  a  necessity  for  verifying  the  fact. 

If  I  is  near  unity,  then  the  coefficient  g  —  I  becomes  small  ; 
and  although  the  corresponding  term  in  u  is  still  of  the  third 
order,  it  becomes  of  the  second  order,  and  therefore  impor- 
tant in  tj  Art.  (30). 

If  I  is  not  near  0,  1  or  2,  the  term  will  not  affect,  to 
that  order,  the  values  of  u  or  of  t. 

48.   Again,  it  is  necessary  to  have  the  value  of  u  correct 

T 

to  the  second  order,  before  we  can  compute  7^-3  to  the  fourth 

order,  with  a  full  assurance  that  we  have  secured  all  terms 
of   that    order    in    which    the    coefficient    of   the    argument 
is   nearly  zero*;  —  such   terms  becoming,  as  we  have  seen 
Art.  (30),  important  in  the  value  of  t. 
For,  assume  that 

u  =  a  (1  +  e  cos(c0  -  a)  +  TFJ, 


.         T  .    ,. 

tbeny;-3=-    .,  4    '  TFFTiSin    2- 

AV        2AV  {1  +  e  cos  (ct)  -  a)  +  W£ 

=  -  fro"  (1...-  4  W2)  sin  ((2  -  2m)  9  -  20}. 

Therefore  if  Wz  contain  a  term  It  cos(r0  —  p)  of  the  second 
order  where  r  is  near  2    (except  2  -  2m)  there  will  be  a 


*  This  was  pointed  out  by  Mr.  Walton  in  the  Quarterly  Journal  of  Mathematics, 
Vol.  ix.,  p.  227 ;  and  in  this  edition  the  text  had  been  modified  in  accordance 
with  his  suggestion ;  but  the  method  given  by  him  for  finding  the  value  of  u  is 
defective  on  account  of  its  not  mentioning  the  necessity  for  the  previous  deter- 
mination of  a. 
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T 

corresponding  term    in  -p— 3  of  the  fourth  order  with  a  small 

il  U 

coefficient   of  0,  and   therefore    an    important   term    in   the 
value  of  t. 

49.  For  the  sake  of  simplifying  the  expressions  we  shall 
generally  omit  cc,  /3,  and  7  in  the  arguments.  If  we  remark 
that  c  and  a  always  appear  together  in  the  form  cO  -  a,  it 
will  be  sufficient  to  write  c0,  instead  of  cO  —  a. ;  and,  for 
the  same  reason,  we  may  replace  m0  +  ft  'and  gO  —  y  by 
7W0,  and  ##,.  The  suffix  (J  will  remind  us  of  the  omission, 
and  we  may  at  any  stage  re-introduce  the  omitted  symbols. 


To  compute  s  to  the  second  order. 

50.  Since   we   must   first   determine   the  value  of  s,  it 
becomes  necessary  to  examine  the  differential  equation  (B) 
to  see  whether  any  difficulty  of  a  similar  kind  to  those  we 
have    been    here    considering    may     arise,    or   whether   the 
elements   on   which    the    computation   depends   are    already 
known  with   sufficient  accuracy.     We  fortunately  find  that 
all  the  terms  of  the  right-hand  member  of  the  equation  are 
of   the  third    order,  and   therefore  that   the  values    already 
obtained  will  ensure  accuracy  to  the  second  order,  and  we 
may  safely  begin  our  work  here. 

51.  To  the  first  order  we  have 

s  —  k  sin^, 

ds      .  Q      .  (because  a  - 1 

^=^COS^'  =  *COS^tto  first  order, 

d*s 
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In  the  expression  for 


,  which  is  of  the  third  order, 


it  will  be  sufficient  to  put  a  for  u,  and  a  for  u,  then 


j  +  f  ra*&  sin  (2  -%m—g]  01? 

the  other  term  of  the  third  order  being  rejected  because  the 
coefficient  of  6  is  not  near  unity. 

The  next  term  of  equation  (B)  is  TF^^)  an^  as  i^  *s 

T 
of   the  first  order,  it   will  be  sufficient  to  have  -7^—3  correct 

to  the  second  order, 

^  —  K™  (2  -*«)',• 
T      ds 

~tfu*'~de  =  ~  %m*k  sin  ^  ~~ 

=  -  %m*k  sin  (2  -  2m  -g]  6^ 

neglecting  the    other  term  of  the  third    order  because   the 
coefficient  of  6  is  not  near  unity. 

The   last   term   of  equation    (B)   is   2(^72+5)  /»-•<" 

.   rp 

which  =0  to  the  third  order,  since   Ij^dd   depending    on 

the  disturbing  force  is    of  the  second  order,  and  the  other 
factor  =  0  to  the  first  order. 


The  equation  in  s  becomes 


2k  sin  (2  -  2m  -  g]  6^ 


To  solve  it,  assume 


^ 
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Substitute,  and  equate  coefficients  of  like  terms, 

£  (1  -  f)  =  -  fyrfk  ;  /.  g  =  (1  +  f  m^  =  1  +  f 

Therefore,  accurately,  to  the  second  order,* 
s  =  k  sin  (gd  -  7)  4  f  mk  sin  {(2  -  2?n  -  g]  6  -  2/3  +  7}  .../%. 
The  term  of  the  second  order,  having  a  coefficient  near 
unity,  need  not  be  taken  into  account  in  finding  the  value 
of  u  (47)  ;   but  it  will  have  to  be  brought  in  subsequently 
to  determine  t. 

To  determine  u  to  the  second  order. 

52.  We  must  compute  the  right-hand  member  of  the 
equation  (A),  Art.  (31),  accurately  to  the  second  order,  and 
include  those  terms  of  the  third  order  which  have  a  coefficient 
near  unity. 

To  the  first  order  we  have 

u  —  a  (\  +  e  coscflj), 

du  (because  c  =  1 

-j-Q=-  —  aec  sin  cB.  =  —  ae  sin  c#,  •{      „          , 
d6  '   (to  first  order, 

eftt 

df+u  =  a. 

From  the  expressions  (D)  and  (£),  Art.  (31),  we  get 


f     ........ 


*  No  complementary  term  P  cos  (0  —  Q)  is  added  ;  for,  though  by  the  theory 
of  differential  equations,  this  would  form  a  necessary  part  of  the  solution,  we 
have  seen  Art.  (39)  that  it  cannot,  in  this  shape,  form  a  part  of  the  correct  value, 
but  will  be  comprised  in  the  terms  whose  argument  is  g&  —  y. 
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1  p 

ah2u? 
-  %m2  (1  +  3e'  cos(7w01  -  f )  -  3e  cosc0,}  {1+3  cos(2  -  2m)  0J, 

1   P 

-  _  =  l  -  |F  +  JA;2  cos  2#0X  -  im2  -  |m2  cos  (2  -  2m)  0, 
ci  fi  u 

+  |m2e  COSC0J  +  |m2e  cos  (2  —  2m  —  c)  0t. 
The  other  terms  of  the  third  order,  which  arise  in   the 
development    of  the  expression,  are   neglected  because   the 
coefficients  of  0  are  not  near  unity, 

T 


=  -  f w2  [1  +  3e'  co*(m0t  -  f )  -  4e  cosc^J  sin  (2  - 
=  -  |ma  sin  (2  -  2??z)  0l  +  3m*e  sin  (2  -  2m  -  c)  ^, 

fa  3  •  TZi =  |^2^e  sin  (2  —  2m)  0.  sinc^, 

7iW         6/^ 

=  f  m*ae  cos  (2  -  2m  —  c)  0, ; 
rejecting  the  other  terms  of  the  third  order  according  to  (30). 


=  |m*  cos  (2  -  2m)  0X  -  3m2e  cos  (2  -  2m  -  c)  6^ 


2 


The  equation  in  u  becomes 
1    d*u 


-  3m2  cos  (2-  2m)  ^ 

+  y  wi*«  cos  (2  -  2m  -  c)  ^. 


To  integrate  this,  assume 


1+A 

+  e  cos  c^t 


+  C  cos  (2-  2m)  0, 
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Substitute,  and  equate  the  coefficients  of  like  terms 


e  (1  -c2)  = 


C  {1-  (2-  2wi/}=  -3m2, 

D  {1  -  (2  -  2w  -  c)2}  =  y* 
Therefore 


c=(l  - 


3m2 


u  —  a 


-f  e  cos(c#  —  a) 


—  J#2cos2  (##  —  7) 
+  m2  cos  }(2  -  2m)  6  -  2/3} 
+  Vme  cos  {(2  -  2m  -  c)  0  -  2/3  +  a 
This  is  the  accurate  value  of  u  to  the  second  order. 

To  find  t  to  the  second  order. 
53.   The  equation  to  be  considered  is 


and,  making  use  of  the  values  of  s  and  u  obtained  to  the 
second  order,  we  must  retrace  our  steps  to  recompute  u  and 

T 
p-g  so  as  to  include  those  terms  of  the  third  order  in  u  in 

which  the   coefficient   of  the   argument   is   near   zero,  and 

T 

similar  terms  to  the  fourth  order  in  ,2  3  Art.  (48).     To  do 

fi  n 

this  we  need  not  repeat  the  whole  work,  but  only  so  much 
of  it  as  will  enable  us  to  pick  out  the  additional  terms  in 
question. 
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l  sin  (2-  2m-  0)  0r  f  wiV 
=  1  .  .  .-  T9Fm&2  cos  (2  —  2m  -  2#)  0t  —  f-mV  cos  (m0j  -  £), 
T          3m  Vs  (1  +  e  cosm0l  -  £)3  sin  (2  -  2m)  0t 
IV  ~  2AV  (1+6  cosc^  -  P2  cos  2^)* 

=  -  f  m2  sin  (2  -  2m)  dl  [I  +  e  cosc0l  -  ±k*  cos  2^^}"* 
=  -  f  m2  sin  (2  -  2m)  0,  {(1  +  e  cosc^)"4  +  ^  cos2#0J 
=  -fm2  sin(2  -2m)  0,  {..,10e2  cos'c^  +  A;2  cos2^0J 
=  -fm2  sin(2-2m)0t{...5e3  0082^  +  ^*  cos2#0,} 
=  _  i_5  m  V  sin  (2  -  2m  -  2c)  0t-  f  wT  sin  (2  -  2m  -  2#)  0r 
These  are  the  only  important  additional  terms. 

T 

--^  dQ=-  f-g  me3  cos  (2-  2m  -  2c)  d^-  |m^2  cos  (2-  2m  -  20)  017 

because  2-2m-2c=-2m  +  fm2,  and  2-2m-20=~2m-fm8a 
Substituting  in  the   differential   equation   for  w,  the  new 

terms  in  y^  -j-n  will  be  of  the  fourth  order   and  must  be 
h  u   au 

neglected,  but  2(jZ2+w)  Ifflf^  w^  bring  in  the   two 

terms 

—  l£me*a  cos  (2  -  2m  -  2c)  6l  and  -  f  m&2«  cos  (2  -  2m  -  20)  0t  • 

and  we  get,  so  far  as  important  additional  terms  are  concerned, 

- 

+  (-  A  +  f  )  rf  cos(2  ~  2m  -  20)  0,-  f  mV  cos  (m0r-  f), 
and  the  corresponding  terms  in  u  will  also  be  of  the  third 
order,  viz. 

-=...yW*  cos  (2-  2m  -2c)  0t  -f  -j^rf  cos  (2  -  2m-  2^)  0, 

—  |mV  cos(m01  —  f). 

H 
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The  complete  value  of  u  therefore,  so  far  as  important  terms 
are  concerned  for  determining  the  value  of  t,  is 

-  =  1  -  3£'  -  lm2  + 

-f  m*  cos  (2  -  2m)  0l  +  ^  me  cos  (2  -  2m  -  c)  0l 

+  ^m2  cos  (2  -  2m  -  2c)  Ol  +  -fymJc'  cos  (2  -  2m  -  20)  0t 


r  T 

and  the  complete  value  of  lj_^-s  dO  will  be 
r  T 

\TT-*dd  =  lm*  cos(2  -  2m)  0,  -  V™3  cos(2  -  2m  -  2c)  0t 

—  f  mF  cos  (2  —  2m  —  20)  0^ 

the  other  term  of  the  third  order  has  no  influence  on  £,  the 
coefficient  of  0  not  being  small. 

54.   We  are  now  in  a  position  to  find  the  value  of  t  to 
the  second  order: 

dd      hu*  \          J  h2u3      /         .hu2  \       J  h2u?      J 

If  -  =  1  +  e  COSC0,  +  W  where   W  contains  all  terms  of  the 
a 

second  and  third  orders  in  the  expression  above,  then 


_ 

=  ^-,{l-2(e  cosc^+TF)  +  3  (e  cosc^-f  TF)2-...} 

=  ^-.  (I  -  2e  cosc^  -  2  TF+  3e2  cos' 

ha  ^ 

=        (1  +  f  e'  -  2*  cosc^-  2  TF+  |e8 
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The  only  part  of  QeW  coscfl,  which  gives  a  term  of  the  third 
order  in  which  the  coefficient  of  the  argument  is  small,  is 

Qe  cosc^  }...+  ^me  cos(2  -  2m  —  c)  6l ...}, 
of  which  one  term  is 

cos  (2-  2m-2c)0 '. 


Therefore  T—.  is  equal  to 

tlU 

1  +  |e2  -f  f  ft  +  m*  —  2e  cosc0l  +  f  e2  cos  2^ 
-f  \ A2  cos2g0l  -  2m9  cos (2  -  2m)  0t 

cos  (2  -  2m  -  c)  0,  +  3m V  cos  (m#t  -  f) 
'  cos  (2  -  2m  -  2c)  ^—  f  m&2  cos  (2  -  2m  —  2g)  0t. 
Also  from  above 

f  T          (I  —  Jm2  cos  (2  —  2m)  0t+  y* me2  cos  (2  —  2m  —  2c)  0t 
J  h2u3         (  +  Jm&2  cos  (2  —  2m  —  2g)  0^ 

We  have  now  to  multiply  these  results  together,  and  we  see 
that  the  terms  having  for  coefficients  ^mea  and  fm&a  will 
disappear  in  the  product.*  These  were  the  terms  which 

T 

were  originally  of  the  fourth  order  in  ^,—3 . 

hu 

r  1-f  f  e3  +  f  &2  +  m*  -  2e  cosc0l  + 
—  =  j~, t  j     4  \kz  co$2g0l  —  141ma  cos (2  — 

I    —  ifme  cos  (2  —  2m  —  c)0l-\-  3mV  cos^^  —  f ). 

*  The  principle  of  the  superposition  of  small  disturbances  might  have  led  us 
a  priori  to  anticipate  the  disappearance  of  these  terms: — for,  when  approximating 
to  the  second  order  we  may,  according  to  this  principle,  calculate  separately  the 
effect  of  the  disturbing  action  of  the  sun,  which  is  of  the  second  order,  and  that 
of  any  other  independent  disturbance  of  the  same  or  of  a  higher  order,  and  then 
add  the  results.  No  term  therefore  which  bears  a  trace  of  the  action  of  both 
these  forces  should  present  itself  to  the  second  order. 

Now  since  neither  k  nor  e  enters  into  the  argument  of  any  term,  it  is  obvious 
that  a  term  whose  coefficient  contains  Jc2  or  e2  at  any  stage  will  retain  them 
through  all  subsequent  operations ;  and  if  the  terms  having  (2  —  2m  —  2g)  01 
or  (2  —  2m  —  2  c)  6,  for  argument  could  remain  in  the  result  to  the  second  order 


V  m3  sin  {(2  -  2m)  0  -  2/3} 
Ymesin{(2-2m-c)0-2y8+a} 
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Let  jJL  (1  -f  |g»  +  f  &2  +  m2)  =-  , 

therefore  ^>  =  ^a2  (1  —  f  e2  —  f  &a  —  ?w2)  to  the  third  order ; 
therefore,  multiplying  by  p  and  integrating,  we  get,  still  to 
the  second  order, 
pt=0-2e  sin  (c0  -  a)  +  f  e2  sin  2  (c0  -  a) 


,..©„ 


no  constant  is  added,  the  time  being  reckoned  from  the 
instant  when  the  mean  value  of  6  vanishes,  for  the  reasons 
explained  in  (Art.  36). 

Coordinates  expressed  in  terms  of  the  time. 

55.  The  preceding  equations  U^  S^  02  give  the  reci- 
procal of  the  radius  vector,  tbe  latitude  and  the  time  in 
terms  of  the  true  longitude  ;  but  the  principal  object  of  the 
analytical  investigations  of  the  Lunar  Tbeory  being  the 
formation  of  tables  which  give  the  coordinates  of  the  moon 
at  stated  times,  we  must  express  w,  s,  and  6  in  terms  of  t. 

To  do  this,  we  must  reverse  the  series  pt  =  6  —  &c.,  and 
then  substitute  the  value  of  6  in  the  expressions  for  u  and  s. 

Now  0  =pt  +  2e  sin  (c0  -  a)  to  the  first  order 


jbhey  would  there  have  A2  or  e2  respectively  -in  the  coefficient,  since  they  have 
them  when  they  first  appear.  Bnt  &2,  e2  and  m  being  independent,  such  terms 
must  have  arisen  from  the  combination  of  others  which  originally  (in  the 
expressions  of  the  forces)  involved  these  quantities  separately  :  —  terms  like 
k-  cos  2  (gQ  —  y)  and  e2cos2(c0  —  a)  combined  with  -terms  involving  m:  that 
is,  terms  which  are  representatives  of  disturbing  forces  of  the  second  order  (see 
note  p.  31)  combined  with  others  depending  on  the  sun's  action.  Such  terms 
must  therefore  be  of  an  order  beyond  the  second. 
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therefore  cd  —  a.  =  cpt  —  a  +  2e  sin  (cpt  —  a)  to  the  first  order, 
2e  sin  (cd  —  a)  =  2e  {sin  (cpt  —  a)  +  2e  sin  (cp£  —  a)  cos  (cpt— a)} 

to  the  second  order? 
=  2e  sm(cpt  —  0i)  +  2e*  sin2  (cpt—  a) 5 

and  as  6  and  jp£  differ  by  a  quantity  of  the  first  order,  they 
may  be  used  indiscriminately  in  terms  of  the  second  order ; 
therefere 

6=pt+2e  sin  (cpt  —  a)+  f  e*  sin  2  (cpt  —  a) 


w 


sin  {(2  -  2m)  pt  -  2/3} 


...©v 


—  3me  sin  (wp£  -f  /3  — 


56.  In  the  value  of  u  given  in  Art.  (52),  substitute  pt  for 
6  in  terms  of  the  second  order,  and  pt  +  2e  siu(cpt-  a)  in  the 
term  of  the  first  order ;  then 

.— f&2—  ^m2—  es+e  cos  (cpt— a) +e3  cos2(cpt—oi) 
—  J&2  cos  2  (gpt  —  7) 
+  ??i2  cos  {(2  -  2m) pt  —  2/3} 


or 


57.    Similarly,  the  expression  for  s  becomes 
—  7)  4  2e  sin  (cpt  -  a)} 

+  §mk  sin  {(2  —  2m  — . 
sin  (gpt  —  7) 
+  e  sin  {(g  +  c)pt  -  a  -  7} 


-2/3  +  7}; 


4-  §w  sin }(2  -2m-g')pt-2/3  +  7} 

The  expression  for  s  is  more  complex  in  this  form  than  when 
given  in  terms  of  the  true  longitude  0. 


54  LUNAR  THEORY. 

Moon's  Parallax. 

58.    If  P  be  the  moon's  mean  parallax,  and  n  the  parallax 
at  the  time  £, 

_     radius  of  earth  R 

n=  distance  of])   =  T~  ~  =  *«l  "  &  t0  the  third  °rder> 


=  Ru  {1  —  J&2  +  j&2  cos  2  (gpt  -  7)}  to  the  second  order, 

{1  -  k*  —  \m*  -  e*  -f  e  cos  (cp£  -  a)  +  e*  cos2  (cpt  -  a) 
+  w2  cos  {(2  -2^)^-2/3} 
-j-  i-gme  cos  {(2  -  2m  -  c)pt  -  2/3  -f  a)  ; 
but  P=  the  portion  which  is  independent  of  periodical  terms, 


1  +  e  cos  (cpt  —  a)  +  e2  cos  2  (cpt  —  a) 
therefore  n  =  P      +  m*  cos  J(2  -  2m)  pt  -  2/3} 

+  if  me  cos  j(2  -  2m  -c)pt-2ft  +  a} 

neglecting  terms  of  the  third  order.  We  see  that,  to  the 
second  order,  the  variable  part  of  the  parallax  is  independent 
of  the  inclination. 

59.  Here  we  terminate  our  approximations  to  the  values 
of  w,  s,  and  6.  If  we  wished  to  carry  them  to  the  third 
order,  it  would  be  necessary  to  include  some  terms  of  the 
fourth  and  fifth  orders  according  to  Art  .(30),  and  the  ap- 
proximate values  of  P,  T,  and  $,  given  in  Art.  (24),  would 
no  longer  be  sufficiently  accurate,  but  we  should  have  to 
recur  to  the  exact  values,  and  from  them  obtain  terms  of 
an  order  beyond  those  already  employed.*  The  process 
followed  in  the  preceding  pages  is  a  sufficient  clue  to  what 
would  have  to  be  done  for  a  higher  approximation. 

*  See  Parallactic  Inequality,  Art.  (105). 
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The  coordinates  u  and  6'  of  the  sun's  position  are,  by 
the  theory  of  elliptic  motion,  known  in  terms  of  the  time  £, 
and  t  is  given  in  terms  of  the  longitude  6  by  the  equation 
@2.  Hence  u  and  &  can  be  obtained  in  terms  of  0;  but 
it  will  be  necessary  to  take  into  account  the  slow  progressive 
motion  of  the  sun's  perigee,  which  we  have  hitherto  neglected. 
This  we  may  do  by  writing  c&  —  £  for  6'  —  f,  c  being  a 
quantity  which  differs  very  little  from  unity.* 

These  values  of  u,  6',  together  with  those  of  u  and  s 
given  by  U9  and  $2,  are  then  to  be  substituted  in  the 
corrected  values  of  the  forces,  and  thence  in  the  differential 
equations.  The  integrations  being  performed  as  before  will 
give  the  values  of  w,  s,  and  t  in  terms  of  9  to  the  third 
order,  and  from  these,  as  in  Arts.  (55),  (56),  and  (57),  may  be 
obtained  u,  s,  and  6  in  terms  of  t. 

60.  More  approximate  values  of  c  and  g  are  obtained  at 
the  same  time,  by  means  of  the  coefficients  of  cos  (cO  —  a) 
and  sin  (gd  —  7)  in  the  differential  equations,  (see  Appendix, 
Arts.  102  and  103). 

61.  The   values  to  the   fourth  order  are  then  obtained 
from  those  to  the  third  by  continuing  the  same  process,  and 


*  '  En  reflechissant  sur  les  termes  que  doivent  introduire  toutes  les  quantitea 
'  precedentes,  on  voit  qu'il  se  peut  glisser  des  cosinus  de  Tangle  6  dont  nous  avons 
'  vu  le  dangereux  effet  d'amener  dans  la  valeur  de  u  des  arcs  au  lieu  de  leurs 
'cosinus;  de  tels  termes  viendront,  par  exemple,  de  la  combinaison  des  cosinus 
'de  (1  —  m)  6  avec  des  cosinus  de  mQ 

' Pour  eViter  cet  inconvenient  qui  oterait  a"  la  solution  precedente 

'Tavantage  de  convenir  &  un  aussi  grand  nombre  de  revolutions  qu-'on  voudrait, 
'  et  la  priverait  de  la  simplicite  et  de  T  universalite  si  precieuses  en  mathematiques, 
<il  faut  commencer  par  en  chercher  la  cause.  Or,  on  decouvre  facilement  que 
'  ces  termes  ne  viennent  que  de  ce  qu'on  a  suppose  fixe  Tapogee  du  soleil,  ce 
'  qui  n'est  pas  permis  en  toute  rigueur,  puisque  quelque  petite  que  soit  sur  cet  astre 
'  Taction  de  la  lune,  elle  n'en  est  pas  moins  reelle  et  doit  lui  produire  un  mouvement 
'  d'apogee  quoique  tres  lent  &  la  verite'.— Clairaut,  The&i-ie  de  la  Lune,  p.  55. 
2me  Edition. 
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so  on  to  the  fifth  and  higher  orders ;  but  the  calculations  are 
so  complex  that  the  approximations  have  not  been  carried 
beyond  the  fifth  order,  and  already  the  value  of  6  in  terms 
of  t  contains  128  periodical  terms,  without  including  those 
due  to  the  disturbances  produced  by  the  planets.  The  coeffi- 
cients of  these  periodical  terms  are  functions  of  TTZ,  e,  e,  —  , 

c,  g,  k,  and  are  themselves  very  complicated  under  their 
literal  forms :  that  of  the  term  whose  argument  is  twice  the 
difference  of  the  longitude  of  the  sun  and  moon,  for  instance, 
is  itself  composed  of  46  terms,  combinations  of  the  preceding 
constants. 

See  Ponte*coulant,  Systeme  du  Monde,  torn.  IV.,  p.  572. 
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CHAPTER    V. 


NUMERICAL  VALUES  OF  THE  COEFFICIENTS. 

62.  Having  thus,  from  theory,  obtained  the  form  of  the 
developments  of  the  coordinates  of  the  moon's  position  at 
any  time,  the  next  necessary  step  is  the  determination  of  the 
numerical  values  of  the  coefficients  of  the  several  terms. 

We  here  give  three  different  methods  which  may  be 
employed  for  that  purpose,  and  these  may,  moreover,  be 
combined  according  to  circumstances. 

63.  First  Method.   The  values  of  the  constants  /?,  m,  a, 
£,  7,  f  which  enter  into  the  arguments,  and  of  the  additional 
ones  a,  &,  e,  &c.  which  enter  into  the  coefficients,   of  the 
terms  in  the  previous  developments,  may  be  obtained  with 
great  accuracy  from  observation. 

For  this  we  must  employ  observations  separated  by  very 
long  intervals,  such,  for  instance,  as  ancient  and  modern 
eclipses;  also  particular  observations  of  the  sun  and  moon 
made  when  the  bodies  occupy  certain  selected  positions; 
and  other  observations  of  a  special  character. 

When  the  values  of  the  elements  have  been  so  obtained, 
then  the  theoretical  values  of  the  coefficients  may  be  com- 
puted by  substitution  in  the  analytical  expressions. 

64.  Second  Method.   Let  the  constants  which  enter  into 
the  arguments  be  determined  as  in  the  first  method ;  and  let 
a  large  number  of  observations  be  made,  from  each  of  which 
a  value  of  the  true  longitude,  latitude,  or  parallax  is  obtained, 

I 
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together  with  the  corresponding  value  of  t  reckoned  from 
the  fixed  epoch  when  the  mean  longitude  is  zero.  Let  these 
corresponding  values  be  substituted  in  the  equations,  each 
observation  thus  giving  rise  to  a  relation  between  the  unknown 
constant  coefficients. 

A  very  great  number  of  equations  having  been  thus 
obtained,  let  them,  by  the  method  of  least  squares  or  some 
analogous  process,  be  reduced  to  as  many  as  there  are 
coefficients  to  be  determined.  The  solution  of  these  simple 
equations  will  give  the  required  values. 

This  method,  however,  would  scarcely  be  practicable  in 
a  high  order  of  approximation.  For  instance,  in  the  fifth 
order,  as  stated  in  Art.  (61),  each  observation  would  give 
rise  to  an  equation  containing  129  unknown  quantities  and 
the  immense  number  of  equations  so  obtained  would  have 
to  be  reduced  to  129  equations  of  130  terms  each. 

65.  Third  Method.  When  the  constants  which  enter  into 
the  arguments  have  been  determined  by  the  first  method,  we 
may  obtain  any  one  of  the  coefficients  independently  of  all 
the  others  by  the  following  process,  provided  the  number  of 
observations  be  very  great. 

Let  the  form  of  the  function  be 


and  let  it  be  required  to  determine  the  constants  A,  B,  (7,  &c. 

separately  ;  0,  $,  &c.,  being  known  functions  of  the  time. 
Let  the  results  of  a  great  number  of  observations  corre- 

sponding to  values  0l5  02,  03,  &c.,  <£17  <£2,  </>3,  &c.,  be'F^  7a,  F3, 

&c.  ;  so  that 

V^  =  A  +jBsin^+  Osin^-f  &c., 
F2  =  A  +  B  sin  02  +  C  sin  02  +  &c., 
F3  =  A  +  B  sin  03  +  C  sin  ft  +  &c., 

Vn  =  A  -f  B  sin  0  ,  +  C  sin  <£M  -f  &c. 
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Now,  n  being  very  great,  we  may  assume  that  the  sum  of 
the  positive  values  of  each  periodical  term  will  be  about 
counterbalanced  by  the  sum  of  its  negative  values ;  and 
therefore,  that  if  we  add  all  the  equations  together  these 
terms  will  disappear ; 

Y  4  F  +  F  +..,...+  F 

therefore  A  =  — —  — "  . 

n 

which  determines  the  non-periodic  part  of  the  function. 

To  determine  B.  Let  the  observations  be  divided  into 
two  sets  separating  the  positive  and  negative  values  of  sin0; 
then  the  other  periodical  terms,  not  having  the  same  period, 
may  be  considered  as  cancelling  themselves  in  adding  up  the 
terms  of  each  set.  Let  there  be  r  terms  in  the  first  set 

and  ,9  terms  in  the  second,  and  let   F',  F", Fr  be  the 

values  of  F  corresponding  to  positive  values  of  sin  0,  which 
values  we  may  assume  to  be  uniformly  distributed  from 

sin0  to  sin-TT,  and  therefore  to  be  sin  80,  sin  2  80, sin  r.  80, 

where  r .  80  =  TT. 

And,  again,  let  F/5  F//7  F^x, F5,  be  the  values  of  F 

corresponding  to  the  negative  values  of  sin0 ;  viz.,  sin  (-  A0), 

sin(-2A0), sin(-.s.A0),  where  s.A0  =  7r.     Then, 

V/  =A-Bsm&0     +  Osin^  4 


V  = 


Vr  = 
therefore 
F'+  7"+ 


therefore 


<7  sin 


(sin0) 


~ 


7T 


TT 


F,  =  ^1-5  sins.  A0-f 
therefore 


therefore 
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therefore* 


and  in  a  similar  manner  may   each   of  the   coefficients  be 
independently  determined.* 

66.  When  the  periods  of  two  of  the  terms  differ  hut 
slightly  —  for  instance  if  6  and  <f>  go  through  their  periodic 
variations  very  nearly  in  the  same  time,—  the  method  could 
not  then  with  safety  be  applied  ;  for,  since  the  same  values 
of  6  and  <f>  would  very  nearly  recur  together  during  a  longer 
time  than  that  through  which  the  observations  would  ex- 
tend, the  two  terms  would  be  so  blended  in  the  value  of 
F  that  they  would  enter  nearly  as  one  term  —  the  difference 
between  6  and  <f>  would  be  very  nearly  the  same  at  the 
end  as  at  the  beginning  of  the  series  of  observations. 

*  If  r  and  9  are  not  sufficiently  great  to  allow  us  to  substitute  /V71"  sin  6dQ  for 
£</""  sin  0  .  50,  we  must  proceed  as  follows  : 

V  +  V"  +  .....  H-  VT-rA  +  B  (sin  56  +  sin250  +  ......  +  sinr50) 


sin  £50 


therefore 
B  = 


/Vr  +  V"  +       +  Vr      V  +  V  + +  Vg\ 

1       1YI  (  r  ~  / 


_r/i    ^(L   IM  fv+v"+ +  F"    vt  +  vt,+ +  v.\ 

~  4  r +  24  yi + ?;;  I         r  •,       -  /  • 
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67.  Let  us  suppose  the  periods  to  be  actually  identical, 
so  that  (j>  =  6  +  a,  a  being  some  constant  angle  ;  then 


may  be  written  (B+  C  cos  a)  sin0  +  Osina  cos0, 

or         F=  A  +  (B  +C  cos  a)  sin0+  <7sinacos04-  ....... 

If  now  we  divide  the  observations,  as  before,  into  two 
sets,  corresponding  to  the  positive  and  negative  values  of 
sin#,  the  terms  involving  cos#  will  disappear  in  the  sum- 
mation of  each  set  ;  and,  following  the  process  of  the  method, 
we  shall  find  the  value  of 

B  4-  C  cos  a.  —  M  suppose. 

Dividing  again  into  two  sets  corresponding  to  the  positive 
and  negative  values  of  cos0,  the  terms  in  sin#  will  be  can- 
celled, and  the  same  process  will  give  the  value  of 

C  sin  a  =  N  suppose. 

Treating  the  observations  in  the  same  way  with  respect 
to  the  angle  $,  we  get  two  results, 


from  these  four  equations  we  easily  get 


Jf,  Nt  M'j  N'  are  connected  by  the  equation  of  condition, 

M*-M"  =  N«-N\ 

68.  When  the  periods  of  6  and  <j>  are  nearly,  but  not 
exactly,  the  same,  this  equation  of  condition  will  not  hold, 
and  the  preceding  ^values  of  B  and  C  would  not  be  exactly 
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correct  ;  but  yet  they  would  be  very  approximate,  especially 
if  the  mean  between  the  two  values  of  B  be  taken. 

Or,  we  may,  after  having  taken  one  of  these  slightly 
erroneous  values  for  J5,  make  a  further  correction  by  estab- 
lishing as  it  were  a  counterbalancing  error  in  the  value  of 
C.  Let  1?  be  the  value  so  found  for  B^  then,  from  the 
V  of  each  of  the  observations  subtract  the  value  B'  sin  #, 
the  result  27  will  be  very  nearly  equal  to  A  +  C  sin  <£  -t-  &c., 
and  from  the  n  equations 


a  value  C'  of  C  will  be  obtained,  by  the  rule,  which  will  be 
very  approximate,  and,  at  the  same  time,  agree  better  with 
B'  in  satisfying  the  equations  than  C  itself  would  do. 

69.  When  two  terms  whose  periods  are  nearly  equal  do 
occur,  it  is  plain,  by  examining  the  values  of  M  and  M  ',  that 
the  errors  which  would  be  committed  by  following  the  rule, 
without  taking  account  of  this  peculiarity,  would  be  the  taking 
B+  C  cos  a  and  (7+  B  cos  a  for  B  and  C  respectively. 


CHAPTER    VI. 


PHYSICAL  INTERPRETATION. 

70.  The  solution  of  the  problem  which  is  the  object  of 
the  Lunar  Theory  may  now  be  considered  as  effected ;  that 
is,  we  have  obtained  equations  which  enable  us  to  assign 
the  moon's  position  in  the  heavens  at  any  given  time  to 
the  second  order  of  approximation ;  we  have  explained  how 
the  numerical  values  of  the  coefficients  in  these  equations 
may  be  determined  from  observation;  and  we  have,  more- 
over, shewn  how  to  proceed  in  order  to  obtain  a  higher 
approximation.* 

It  will,  however,  be  interesting  to  discuss  the  results  we 
have  arrived  at : — to  see  whether  they  will  enable  us  to  form 
some  idea  of  the  nature  of  the  moon's  complex  motion ;  also 
how  far  they  will  explain  those  inequalities  or  departures 
from  uniform  circular  motion  which  ancient  astronomers  had 
observed ;  but  which,  until  the  time  of  Newton,  were  so 
many  unconnected  phenomena;  or,  at  least,  had  only  such 
arbitrary  connexions  as  the  astronomers  chose  to  assign,  by 
grafting  one  eccentric  or  epicycle  on  another  as  each  newly 
discovered  inequality  seemed  to  render  it  necessary. 

It  is  true  that  our  expressions,  composed  of  periodic 
terms,  are  nothing  more  than  translations  into  analytical 

*  The  means  of  taking  into  account  the  ellipsoidal  figure  of  the  earth  and  the 
disturbances  produced  by  the  planets,  are  too  complex  to  form  part  of  an  intro- 
ductory treatise.  For  information  on  these  points  reference  may  be  made  to 
Airy's  Figure  of  the  Earth.  Pontecoulant's  Systeme  du  Monde,  vol.  IV. 
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language  of  the  epicycles  of  the  ancient  ;*  but  they  are 
evolved  directly  from  the  fundamental  laws  of  force  and 
motion,  and  as  many  new  terms  as  we  please  may  be  ob- 
tained by  carrying  on  the  same  process ;  whereas  the  epicycles 
of  Hipparchus  and  his  followers  were  the  result  of  numerous 
and  laborious  observations  and  comparisons  of  observations ; 
each  epicycle  being  introduced  to  correct  its  predecessor 
when  this  one  was  found  inadequate  to  give  the  position 
of  the  body  at  all  times :  just  as  with  us,  the  terms  of  the 
second  order  correct  the  rough  results  given  by  those  of 
the  first ;  the  terms  of  the  third  order  correct  those  of  the 
second,  and  so  on.  But  it  is  impossible  to  conceive  that 
observation  alone  could  have  detected  all  those  minute 
irregularities  which  theory  makes  known  to  us  in  the  terms 
of  the  third  and  higher  orders,  even  supposing  our  instruments 
far  more  perfect  than  they  are ;  and  it  will  always  be  a 
subject  of  admiration  and  surprise,  that  Tycho,  Kepler,  and 
their  predecessors  should  have  been  able  to  feel  their  way 
so  far  among  the  Lunar  inequalities,  with  the  means  of 
observation  they  possessed. 

LONGITUDE   OF  THE   MOON. 

71.    We  shall  first  discuss  the  expression  for  the  moon's 
longitude,  as  found  Art.  (55), 

0  =pt  +  2e  sin  (cpt  -  a)  +  \e*  sin  2  (cpt  -  a) 
+  *£me  sin  {(2  —  2m  -  c)pt  -  2/3  +  a} 
4-  V  wi*  sin  {(2  -  2m)  pt  -  2/3} 
-  3  we'  sin  (mpt  +  /3  -  f) 
—  J&3  sin  2  (gpt  -  7). 

The  mean  value  of  6  is  pt;  and  in  order  to  Judge  of  the 
effect  of  any  of  the  small  terms,  we  may  consider  them 

*  See  Whewell's  History  of  the  Inductive  Sciences. 
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one  at  a  time  as  a  correction  on  this  mean  value  pt,  or 
we  may  select  a  combination  of  two  or  more  to  form  this 
correction. 

We  shall  have  instances  of  combination  in  the  elliptic 
inequality  and  the  evection.  Arts.  (73)  and  (77)  ;  but  in  the 
remaining  inequalities  each  term  of  the  expression  will  form 
a  correction  to  be  considered  by  itself. 

72.  Neglecting  all  the  periodical  terms,  we  have 

0~& 
dd 

2T* 

which  indicates  uniform  angular  velocity;  and  as,  to  the 
same  order,  the  value  of  u  is  constant,  the  two  together 
indicate  that  the  moon  moves  uniformly  in  a  circle,  the 

period   of  a   revolution  being  --  ,  which   is,  therefore,  the 

expression  for  a  mean  sidereal  month,  or  about  27  J  days.* 
The  value  of  jp  is,  according  to  Art.  (54),  given  by 


and  as  m  is  due  to  the  disturbing  action  of  the  sun,  we  see 
that  the  mean  angular  velocity  is  less,  and  therefore  the  mean 
periodic  time  greater  than  if  there  were  no  disturbance. 

Elliptic  Inequality  or  Equation  of  the  Centre. 

73.  We  shall  next  consider  the  effect  of  the  first  three 
terms  together  :  the  effect  of  the  second  alone,  as  a  correction 
of  pi,  will  be  discussed  in  the  Historical  Chapter,  Art.  (115). 


*  The  accurate  value  was  27d.  7h.  43m.  11  '26  Is.  in  the  year  1801.     See 
Art.  (108), 

K 


66  LUNAR  THEORY. 

0  =pt  +  2e  sin  (cpt  -  a)  +  \<?  sin  2  (cpt  -  a), 
which  may  be  written 


But  the  connexion  between  the  longtitude  and  the  time  in 
an  ellipse  described  about  a  centre  of  force  in  the  focus  ia, 
Art.  (13),  to  the  second  order  of  small  quantities  : 

6  s=  nt  +  2e  sin  (nt  -  a')  +  |e2  sin  2  (nt  —  a'), 

where  n  is  the  mean  motion,  e  the  eccentricity,  and  a'  the 
longitude  of  the  apse.* 

Hence,  the  terms  we  are  now  considering  indicate  motion 
in  an  ellipse  ;  the  mean  motion  being  p,  the  eccentricity  e, 
and  the  longitude  of  the  apse  a  +  (l  —  c)  pt\  that  is,  the 
apse  is  not  stationary  but  has  a  progressive  motion  in  longi- 
tude, uniform,  and  equal  to  (1  —  c)p. 

74.  The  two  terms  2e  sin  (cpt  -  a)  +  \£  sin  2  (cpt  -  a)  con- 
stitute the  elliptic  inequality,  and  their  effect  may  be  further 
illustrated  by  means  of  a  diagram. 

Let  the  full  line  AMB  re- 
present the  moon's  orbit  about 
the  earth  -E1,  when  the  time  t 
commences,  that  is,  when  the 
moon's  mean  place  is  in  the  prime 
radius  -ET,  from  which  the  lon- 
gitudes are  reckoned. 

The     angle    TEA,    the    longi- 

tude of  the  apse,  is  then  a.  At  the  time  £,  when  the  moon's 
mean  longitude  is  VEM=pt,  the  apse  line  will  have  moved 


*  The  epoch  «  which  appears  in  the  expression  of  Art.  (13)  is  here  omitted  ; 
a  proper  assumption  for  the  origin  of  *,  as  explained  in  Art.  (36),  enabling  us 
to  avoid  the  E. 
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in  the  same  direction  through  the  angle  AEA'  =  (1  —  c)  VEM, 
and  the  orbit  will  have  taken  the  position  indicated  by  the 
dotted  ellipse.  The  true  place  of  the  moon  in  this  orbit, 
so  far  as  these  two  terms  are  concerned,  will  be  wi,  where 

MEm  =  2e  sin  (cpt  —  a)  +  Je2  sin  2  (cpt  —  a) 
= 2e  sin  A' EM  +  |e2  sin  2  A 'EM  . 
=  2e  amA'EM(l  +  Je  cos  A  EM) ; 

which,  since  e  is  about  -fa,  is  positive  from  perigee  to  apogee, 
and  therefore  the  true  place  before  the  mean;  and  the 
contrary  from  apogee  to  perigee:  at  the  apses  the  places 
will  coincide. 

75.  The  angular  velocity  of  the  apse  is  (1  —  c)  p,  or,  if 
for  c  we  put  the  value  found  in  Art.  (52),  the  velocity  will 
be  f  m2p.  Hence,  while  the  moon  describes  360°,  the  apse 
should  describe  £  w*.360°  =  lf°  nearly,  m  being  about  j1^. 

But  Hipparchus  had  found,  and  all  modern  observations 
confirm  his  result,  that  the  motion  of  the  apse  is  about  3°  in 
each  revolution  of  the  moon.  See  Art.  (118). 

This  difference  arises  from  our  value  of  c  not  being  repre- 
sented with  sufficient  accuracy  by  1  —  jma. 

Newton  himself  was  aware  of  this  apparent  discrepancy 
between  his  theory  and  observation ;  and  we  are  led,  by  his 
own  expressions  (Scholium  to  Prop.  35,  lib.  ill.  in  the  first 
edition  of  the  Prmcipia)y  to  conclude  that  he  had  got  over 
the  difficulty*  This  is  rendered  highly  probable  when  we 
consider  that  he  had  solved  a  somewhat  similar  problem 
in  the  case  of  the  node ;  but  he  has  nowhere  given  a  state- 
ment of  his  method :  and  Clairaut,  to  whom  we  are  indebted 
for  the  solution,  was  on  the  point  of  publishing  a  new  hy- 
pothesis of  the  laws  of  attraction,  in  order  to  account  for 
it,  when  it  occurred  to  him  to  carry  the  approximations  to 
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the  third  order,  and  he  there  found  the  next  term  in  the 
value  of  c  nearly  as  considerable  as  the  one  already  obtained. 
See  Appendix.  The  value  is 


.-.  (1  -  c)  360°  =  (1  -f  T7^)  (value  found  previously) 
=  2f°  nearly, 

thus  reconciling  theory  and  observation,  and  removing  what 
had  proved  a  great  stumbling-block  in  the  way  of  all  as- 
tronomers.* 

When  the  value  of  c  is  carried  to  higher  orders  of  ap- 
proximation, the  most  perfect  agreement  is  obtained. 

The  motion  of  the  apse  line  is  considered  by  Newton 
in  his  Principle  lib.  I.,  Prop.  66,  Cor.  7. 

Evection. 

76.  The  next  term  +  ^mesin{(2  -  2m-  c)pt-  2/3  +  «} 
in  the  value  of  6  has  been  named  the  Evection.  We  shall 
consider  its  effect  in  two  different  ways. 

First,  by  itself,  as  forming  a  correction  on  pt, 


Let    ])  =          pt          =  moon's  mean  longitude  at  time  £, 
O  =      mpt  +  @     =sun's  ............................  .....  ? 

a!  =  (1  -  c)  pt  +  a  =  mean  longitude  of  apse    ....  ------  , 


then 


6  =pt  +  *£me  sin  [2  {pt  -  (mpt  + 13)}  -  {pt  -(l-c)pt  +  a}] 
=pt  +  lfme  sin  (2  ())  -  Q)  -  Q)  -  a')}. 

*  See  Dr.  Whewell's  Bridgewater  Treatise, 
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The  effect  of  this  term  will  therefore  be  as  follows  : 
In  syzygies 

0  =pt  -  l-£me  sin  (})  -  a')  ; 

or  the  true  place  of  the  moon  will  be  before  or  behind  the 
mean,  according  as  the  moon,  at  the  same  time,  is  between 
apogee  and  perigee  or  between  perigee  and  apogee. 
In  quadratures 

6  =pt  +  y*  me  sin  (])  —  a'), 

and  the  circumstances  will  be  exactly  reversed. 

In  both  cases,  the  correction  will  vanish  when  the  apse 
happens  to  be  in  syzygy  or  quadrature  with  regard  to  the 
sun  at  the  same  time  as  the  moon. 

In  intermediate  positions,  the  nature  of  the  correction 
is  more  complex,  but  it  will  always  vanish  when  the  sun 
is  at  the  middle  point  between  the  moon  and  the  apse, 
or  when  distant  90°  or  180°  from  that  point;  for  if 

-  r.90°,  where  r  =  0,  ±  1,  or  2, 


2 

sin  [2  (D  -  0)  -  (D  -  a')]  =  sin  ())  +  a'  -  20) 
=  sin  r.  180° 
=  0. 

77.  The  other  and  more  usual  method  of  considering 
the  effect  of  this  term  is  in  combination  with  the  two  terms 
of  the  elliptic  inequality,  as  follows  : 

To  determine  the  change  in  the  position  of  the  apse  and  in 
the  eccentricity  of  the  moon's  orbit  produced  by  the  evection. 

Taking  the  elliptic  inequality  and  the  evection  together, 
we  have 

6  =pt  +  2e  sin  (cpt  -  a)  +  f  e2  sin  2  (cpt  -  a) 

+  ^fme  sin  {(2  -  2m  -  c)pt  -  2/3  -I-  a}. 
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Let  a'  be  the  longitude  of  the  apse  at  time  t  on  supposition 
of  uniform  progression, 

©  sun ;* 

whence  a'  =  ( 1  —  c]  pt  +  a, 

0  =    mpt    -f  /3. 
And  the  above  may  be  written 

6  =*pt  +  2e  sin  (cpt  —  a)  +  {e2  sin  2  (cpt  —  a) 

+  ^me  sin  {cpt -  a  +  2  (a'  -  0)} ; 
and  the  second  and  fourth  terms  may  be  combined  into  one, 


if  we  assume    E  cos  S  =  e  +  l-gme  cos  2  (a'  —  0  ), 
and  E  sin  8  =       *^>me  sin  2  (a'  -  ©)  ; 

whence  tan  8  and  E  may  be  found  ;  and  approximately, 
E=  e  [I  +  ^m  cos2  (a  -  0)}, 


The  term  |e'  sin  2  (cpt  —  a)  may  also,  therefore,  to  the 
second  order,  be  expressed  by 


and  the  longitude  becomes] 

0  =^  +  2E  sin  (cp*  -  a  +  8)  +  %E*  sin  2  (qp*  -  a  +  S), 
or       6  =pt  +  2^  sin  (pt  -  a'  +  S)  +  f  #a  sin  2  (  /rt  -  a1  +  S)  ; 
but  the  last  two  terms  constitute  elliptic  inequality  in  an  orbit 
whose   eccentricity  is  E  and  longitude  of  the  apse  a'  -  8  ; 
therefore   the   evection,   taken    in   conjunction  with  elliptic 
inequality,  has  the  effect  of  rendering  the  eccentricity  of  the 
moon's  orbit  variable,  increasing  it  by  *£me  when  the  apse- 
line  is  in  syzygy,  and  diminishing  it  by  the  same  quantity 
when  the  apse-line  is  in  quadrature  ;  the  general  expression 
for  the  increment  being  *£me  cos  2  (a'—  ©). 
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And  another  effect  of  this  term  is,  to  diminish  the  longitude 
of  the  apse,  calculated  on  the  supposition  of  its  uniform 
progression,  by  the  quantity  8  =  ^m  sin 2  (a'-  O) ;  so  that 
the  apse  is  behind  its  mean  place  when  in  the  first  and  third 
quadrant  in  advance  of  the  sun,  and  before  its  mean  place 
in  the  second  and  fourth. 

The  cycle  of  these  changes  will  evidently  be  completed 
in  the  period  of  half  a  revolution  of  the  sun  with  respect 
to  the  apse,  or  in  about  -f^  of  a  year.* 

78.  The  period  of  the  evection  itself,  considered  inde- 
pendently of  its  effect  on  the  orbit,  is  the  time  in  which 
the  argument  (2  —  2m  —  c}pt  —  2/3  + a  will  increase  by  2?r. 

Therefore  period  of  evection 

2-7T  _  mean  sidereal  month 

(2  —  2m  —  c)  p  2  —  2m  —  c 

mean  sidereal  month 


27J  days 
=  ~Y^~  ,  nearly, 


1  -  2m  +  f  tfi2 
=  31^  days,  nearly,  f 

Newton  has  considered  the  evection,  so  far  as  it  arises 
from  the  central  disturbing  force,  in  Prop.  66,  Cor.  9,  of 
the  Principia. 

Variation. 

79.    To  explain  the  physical  meaning  of  the  term 
V m*  sin  {(2  -  2m) pt  -  2£}, 

*  The  change  of  eccentricity  and  the  variation  in  the  motion  of  the  apse 
follow  the  same  law  as  the  abscissa  and  ordinate  of  an  ellipse  referred  to 
its  centre :  for  if  E  —  e  =  x  and  d  =  y,  then 


t  The  accurate  value  is  31  •  8119  days. 
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in  the  expression  for  the  moon's  longitude, 

0=pt  +  yr/i*  sin  {(2  -  2m)pt-2/3}. 
Let  ])  represent  the  moon's  mean  longitude  at  time 
O  ..................  sun's  ...................  , 

therefore  J)  =pt, 


and  the  value  of  6  becomes 

0=pt  +  V?n8sm2(])-o), 

which  shews  that  from  syzygy  to  quadrature  the  moon's 
true  place  is  before  the  mean,  and  behind  it  from  quad- 
rature to  syzygy;  the  maximum  difference  being  y#i2  in 
the  octants. 

The  angular  velocity  of  the  moon,  so  far  as  this  term 
is  concerned,  is 

7/1 

V  (l-w)w2 


=p  {1  +  V  w8  cos  2  (})  -  o)},  nearly, 

which  exceeds  the  mean  angular  velocity  •  p  at  syzygies,  is 
equal  to  it  in  the  octants,  and  less  in  the  quadratures. 

This  inequality  has  been  called  the  Variation,  its  period 
is  the  time  in  which  the  argument  (2  —  2m)  pt  —  2/3  will 
increase  by  2w  ; 

.   n     «  27r  mean  synodical  month 

.*.  period  of  variation  =  -.  -  :—  =  -  -  - 

(2  -  2m)  p  2 

—  14|  days,  nearly.* 

80.  The  quantity  y  m2  is  only  the  first  term  of  an  endless 
series  which  constitutes  the  coefficient  of  the  variation,  the 
other  terms  being  obtained  by  carrying  the  approximation 

*  The  accurate  value  is  14-765294  days. 
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to  a  higher  order.  It  is  then  found  that  the  next  term  in 
the  coefficient  is  f  f  ra3,  which  is  about  -fj-  of  the  first  terra  ; 
and  as  there  are  several  other  important  terms,  it  is  only 
by  carrying  the  approximation  to  a  higher  order  (the  5th 
at  least)  that  the  value  of  this  coefficient  can  be  obtained 
with  sufficient  accuracy  from  theory.  In  fact,  ym2  would 
give  a  coefficient  of  26'  27"  only  ;  whereas  the  accurate  value 
is  found  to  be  39'  30". 

The  same  remark  applies  also  to  the  coefficients  of  all 
the  other  terms. 

81.  As  far  as  terms  of  the  second  order,  the  coefficient 
of  the  variation  is  independent  of  e  the  eccentricity,  and  k 
the  inclination  of  the  orbit.     It  would  therefore  be  the  same 
in  an  orbit  originally  circular,   whose  plane  coincided  with 
the  plane  of   the  ecliptic:  it  is   thus  that  Newton  has  con- 
sidered it.     Princip.  prop.  66,  cor.  3,  4,  and  5. 

Annual  Equation. 

82.  To  explain  the  physical  meaning  of  the  term 

—  Snie'  sin  (mpt  +  /3  -  f) 
in  the  expression  for  the  moon's  longitude. 
8  =  pt-  3me'  sin  (mpt  +  /3  -  £), 

=pt  —  3me'  sin  (longitude  of  sun  —  longitude  of  sun's  perigee), 

—  pt  —  ome  sin  (sun's  anomaly). 

Hence,  while  the  sun  moves  from  his  perigee  to  his 
apogee,  the  true  place  of  the  moon  will  be  behind  the  mean ; 
and  from  apogee  to  perigee,  before  it.  The  period  being  an 
anomalistic  year,  the  effect  is  called  Annual  Equation. 

Differentiating  6  we  get 

f?f) 

~T  ~P  I*  ~~  3/wV  cos  (sun's  anomaly)}. 
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Hence,  so  far  as  this  inequality  is  concerned,  the  moon's 
angular  velocity  is  least  when  the  sun  is  in  perigee,  that 
is  at  present  about  the  1st  of  January J  and  greatest  when 
the  sun  is  in  apogee,  or  about  the  1st  oijuly. 

The  annual  equation  is,  to  this  order,  independent  of  the 
eccentricity  and  inclination  of  the  moon's  orbit,  and  therefore, 
like  the  variation,  would  be  the  same  in  an  orbit  originally 
circular.  Vide  Newton,  Principia,  prop.  66,  cor.  6. 

Reduction. 

83.  Before  considering  the  effect  of  the  term 

T? 
--81112(0^-7), 

which,  as  we  shall  see  Art.  (89),  is  very  nearly  equal  to 
the  difference  between  the  longitude  in  the  orbit  and  the 
longitude  in  the  ecliptic,  it  will  be  convenient  to  examine 
the  expression  for  the  latitude  of  the  moon,  and  to  see  how 
the  motion  of  the  node  is  connected  with  the  value  of  g. 

LATITUDE  OF  THE   MOON. 

84.  The  expression  found  for  the  tangent  of  the  latitude,* 
Art.  (51),  is 

s  =  k  sin  (g6  -  7)  4  %mk  sin  {(2  -  2wi  -  #)  0  -  2/3  -f-  7}. 
If  we  reject  all  small  terms,  we  have 
s  =  0, 

or  the  orbit  of  the  moon  coinciding  with  the  ecliptic,  which 
is  a  first  rough  approximation  to  its  true  position. 


*  This  expression  for  the  tangent  of  the  latitude  is  more  convenient  than 
that  which  gives  it  in  terms  of  the  mean  longitude,  Art.  (57),  on  account  of 
the  less  number  of  terms  involved.  See  PontC-coulant,  vol.  iv.  p.  (JoO. 
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85.  Taking  the  first  term  of  the  expansion 

s  =  £sin(#<9-7), 
we  may  write  it 

s  =  k  sin  [0-{y-(ff- 
Let  TNm  be  the  ecliptic,  N  the 
moon's  node  when  her  true  longitude 
is  zero,  and  let  M  be  the  position 
of  the  moon  at  time  £,  m  her  place 
referred  to  the  ecliptic ; 

therefore  TJV=7,  Tm  =  0,   tanlfw  =  s. 

Take  NN'=(g—l)0  in  a  retrograde  direction;  and  join 
MN'  by  an  arc  of  great  circle ; 

then         tan  Mm  =  tan  MN'm  sin  N'm , 

or  a  =  tan  MN'm  sin  [6  -  {7  -  (g  -  1)  6}] ; 

which,  compared  with  the  value  of  s  given  above,  shews  that 
MN'm  =  tan'1^  is  constant,  and  therefore  the  term  k  sin  (<7#— 7) 
indicates  that  the  moon  moves  in  an  orbit  inclined  at  an 
angle  tan~V^  to  the  ecliptic,  and  whose  node  regredes  along 

the  ecliptic  with  the  velocity  (g-l}~j.,  or  with  a  mean 
velocity  (g—l)p. 

86.  Hence  the  period  of  a  revolution  of  the  nodes 

2?r       _  one  sidereal  month 

~(g-i)p~         g-i       7 

but,  from  Art.  (51),  the  value  of  g  =  1  +  f  m2 ; 

,       ,,  -in  one  sidereal  month 

therefore  period  or  revolution  of  nodes  =  - 

fm 

=  6511  days,  nearly. 

This  will,  for  the  same  reason  as  in  the  case  of  the  apse, 
Art.    (75),   be   modified   when   we  carry  the  approximation 
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to  a  higher  degree ;  this  value  of  g  is,  however,  much  more 
accurate  than  the  corresponding  value  of  c,  for  the  third 
terra  of  g  is  small ;  the  value  to  the  third  order  being  (see 
Appendix,  Art.  102) 

#=l+fm2-¥92< 

,    ,  .   ,    c         ,    .        f  ,          ,         one  sidereal  month 

and  the  period  of  revolution  of  the  nodes  = ^—* — 

im2(l-$m) 

=  6705  days,  nearly. 

This  is  not  far  from  the  accurate  value  as  given  by 
observation,  and  when  the  approximation  to  the  value  of  g 
is  carried  to  a  higher  order,  the  agreement  is  nearly  perfect. 

The  true  value  is  6793'39  days,  that  is  about  ISyrs.  7  mo. 

Evection  in  Latitude. 

87.  To  explain  the  variation  of  the  inclination  and  the 
irregularity  in  the  motion  of  the  node  expressed  by  the  term 

+  \mk  sin  {(2  -  2m  -  g)  0  -  2/3  +  7}. 

This  term,  as  a  correction  on  the  preceding,  is  analogous 
to  the  evection  as  a  correction  on  the  elliptic  inequality. 
Taking  the  two  terms  together, 

s  =  k  sin  (ad  -  7)  +  $mk  sin  {(2  -  2m  -g)  6  -  2/3  +  7}, 
Let  ])  =  longitude  of  moon  =  0, 

0  = sun    =m0  +  /3, 

Q= node   =7-(#-l)0; 

therefore    s  =  k  sin  Q)  -  &)  -f  %mk  sin  {)>-&-  2  (0  - 
Now  these  two  terms  may  be  combined  into  one, 

*  =  jrsin())-  Q  -A), 
if  /fcos  A  =  k  +  %mk  cos2  (0  -  Q), 

.AT  sin  A  =        *mk  sin2(O- 
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whence  A  and  K  may  be  found  ;  and,  approximately, 


A  =  frasin2(O  -  Q), 
but  the  equation 

s  =  JTsinQ)-  Q  -  A) 

represents  motion  in  an  orbit  inclined  at  an  angle  tan~\ZT  to 
the  ecliptic,  and  the  longitude  of  whose  node  is  Q  +  A. 

This  term  has  therefore  the  following  effects  : 

1st.  The  inclination  of  the  moon's  orbit  is  variable,  its 
tangent  increases  by  %mk  when  the  nodes  are  in  syzygies, 
relatively  to  the  sun,  and  decreases  by  the  same  quantity 
when  they  are  in  quadrature  ;  the  general  expression  for 
the  increase  being  fw&  cos2(o  —  &). 

2nd.  The  longitude  of  the  node,  calculated  on  supposition 
of  a  uniform  regression,  is  increased  by  A  =  §m  sin  2  (o  —  ££•),- 
so  that  the  node  is  before  its  mean  place  while  in  the  first 
or  third  quadrant  in  front  of  the  sun,  and  behind  it  in  the 
second  and  fourth.  Principia,  book  ill.,  props.  33  and  35. 

The  cycle  of  these  changes  will  be  completed  in  the 
period  of  half  a  revolution  of  the  sun  with  respect  to  the 
node,  that  is,  in  173'21  days,  not  quite  half-a-year. 

88*  The  tangent  of  the  latitude  has  here  been  obtained  * 
if  we  wish  to  have  the  latitude  itself  it  will  be  given  by  the 
formula 

latitude  =  s  -  J^3  +  is5  -  &c. 

which,  to  the  degree  of  approximation  adopted,  will  clearly 
be  the  same  as  s. 

Reduction. 

89.   We   may    now    consider  the  term   which    we    had 
neglected  (Art.  83)  in  the  expression  for  the  longitude,  namely, 
—  J&2  sin  2  (gpt  —  7). 
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Let   N  be   the   position   of  the  node  when  the  moon's 
longitude  is  6,  M  the  place  of   the 
moon,   m   the   place   referred   to  the 
ecliptic. 

Therefore     Tm  =  0, 


Nm  —gO  — 


The  right-angled  spherical  triangle  NMm  gives 

,T     tan  Nm 

cos  N  = 

,       „  - 

therefore  r=.     --  ^-j  — 

1  +  cos  AT      tan  NM+  ta 

2A7_  sin  (JVJf  -  JVm) 
1       ~  ' 


tan 
1  -  co$N     tan  NM  -  tan  Nm 


or,  since  both  N  and  NM—  Nm  are  small, 
tanW      JVJf  -  Nm 


therefore  -A^M—  Nm=£k*  sin  2(^^—7)  =J^2  sin  2(gpt—  7),  nearly. 

Hence  this  term,  which  is  called  the  reduction,  is  approxi- 
mately the  difference  between  the  longitude  in  the  orbit  and 
the  longitude  in  the  ecliptic. 

RADIUS  VECTOR. 

90.  To  explain  the  physical  meaning  of  the  terms  in  the 
value  of  u. 

We  shall,  for  the  explanation,  make  use  of  the  formula 
which  gives  the  value  of  u  in  terms  of  the  true  longitude, 
Art.  (52). 
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Neglecting  the  periodical  terms,   we  have  for  the  mean 
value 

w  =  a(l-p2-|w2). 

The  term  —  |mj,  which  is  a  consequence  of  the  disturbing 
effect  of  the  sun,  shews  that  the  mean  value  of  the  moon's 
radius  vector,  and  therefore  the  orbit  itself,  is  larger  than  if 
there  were  no  disturbance. 

Elliptic  Inequality. 

91.    To  explain  the  effect  of  the  term  of  the  fast  order, 
u  —  a  {1  -I-  e  cos  (cO  —  a)}, 


=  a  [1+  e  cos0- 

This  is  the  elliptic  inequality,  and  indicates  motion  in  an 
ellipse  whose  eccentricity  is  e  and  longitude  of  the  apse 
cc  +  (l—  c)0;  and  the  same  conclusion  is  drawn  with  respect 
to  the  motion  of  the  apse  as  in  Art.  (73). 

Evection. 

92.    To  explain  the  physical  meaning  of  the  term 
igmea  cos  {(2  -  2m  -  c)  6  -  2/3  +  a}. 

This,  as  in  the  case  of  the  corresponding  term  in  the  longi- 
tude, is  best  considered  in  connection  with  the  elliptic 
inequality,  and  exactly  the  same  results  will  follow. 

Thus,  calling  )),  O,  and  a'  the  true  longitudes  of  the  moon, 
sun,  and  apse  (the  last  calculated  on  supposition  of  uniform 
motion),  these  two  terms  may  be  written, 

u  =  a  [1  +  e  cos  (D  -  a')  -  ifme  cos  {D  -  a'  +  2  (a'  -  O)}] 


where  J£cos8  =  e  +  ^me  cos  2  (a'—  Q), 

EsmS=        ^fme  sin  2  (a'-  o). 
These  are  identical  with  the  equations  of  Art.  (77). 
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Variation. 

93.    To  explain  the  effect  of  the  term  m*a  cos  {(2-2w)  6-  2/3)  , 
u  =  a  [1  +  m2  cos  {(2  -  2m)  9  -  2/3}] 


As  far  as  this  term  is  concerned,  the  moon's  orbit  would 
be  an  oval  having  its  longest  diameter  in  quadratures  and 
least  in  syzygies.  Principia,  lib.  I.  prop.  66,  cor.  4. 

The  ratio  of  the  axes  of  the  oval  orbit  will  be 

—  2  =  IS  nearly,  m  being  *0748. 
See  Principia^  lib.  ill.  prop.  28. 

Reduction. 
94.   The  last  important  periodical  term  in  the  value  of  u  is 

-^  cos  2  (#0-7). 

This  term  expresses  approximately  the  variation  in  the 
difference  between  the  values  of  u  in  the  orbit  and  in  the 
ecliptic. 

For  if  wt  be  the  reciprocal  of  the  value  of  the  radius 
vector  in  the  orbit, 

u^  —  u  cos  (latitude), 


therefore          u  -  ^  =  |ws2  *=  \ak*  sin3  (gd  -  7) 

=  latf  -  latf  cos  2  (gd  -  7) 
»=  const.  —  \ak*  cos  2  (g6  —  7). 

Taking  this  result  in  connexion  with  that  of  Art.  (89) 
we  see  that  the  r  "eduction  ,  so  far  as  periodical  terms  to  the 
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second  order  are  concerned,  is  simply  a  geometrical  con- 
sequencejof  the  inclination  of  the  orbit;  and  that,  if  we 
measured  the  longitude  and  the  radius  vector  along  the 
orbit  instead  of  taking  their  projection  on  the  plane  of  the 
ecliptic,  these  periodical  terms  would  not  appear. 

95,  There  are  no  other  terms  of  the  second  order  in  the 
value  of  u.     The  annual  equation,  which,  in  the  longitude, 
is  of  the  second   order,  is  only  of  the  third   order  in   the 
radius  vector. 

Periodic  time  of  the  Moon. 

96.  We  have  seen,  Art.  (72),  that  the  periodic  time  of 
the  moon  is  greater  than  if  there  were  no  disturbing  force  ; 
but  this  refers  to  the  mean  periodic  time  established  on  an 
interval  of  a  great  number   of  years,  so   that   the  circular 
functions  in  the  expression  are  then   extremely  small  com- 
pared with  the  quantity  pt  which  has  uniformly  increased. 

When,  however,  we  consider  only  a  few  revolutions,  these 
terms  may  not  all  be  neglected.  The  elliptic  inequality  and 
the  ejection  go  through  their  values  in  about  a  month,  the 
variation  and  reduction  in  about  half-a-month ;  their  effects, 
therefore,  on  the  length  of  the  period  can  scarcely  be  con- 
sidered, as  they  will  increase  one  portion  and  then  decrease 
another  of  the  same  month. 

But  the  annual  equation  takes  one  year  to  go  through  its 
cycle,  and,  during  this  time,  the  moon  has  described  thirteen 
revolutions;  hence,  fluctuations  may,  and,  as  we  shall  now 
shew,  do  take  place  in  the  lengths  of  the  sidereal  months 
during  the  year. 

We  have,  considering  only  the  annual  equation,  Art.  (82), 

pt  =  0  +  3me  sin  (md  +  £-?). 

M 
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Let  T  be  the   length   of  the   period;   then,   when  9  is 
increased  by  2-Tr,  t  becomes  t  -f  T' 

therefore  p  (t  +  T)  =  2?r  +  Q  -f  3me'  sin  (2tmr  +  mO  +  &-  f  ), 
whence  £>  T=  2ir  +  6  W  sin  WTT  cos  (WTT  +  w#  +  ft  —  £  )  ; 


therefore  T—  mean  period  +  -  —  smmir  cos(o  - 


where  O  =m6  +  /3  +  rrnr  —  sun's  longitude  at  the  beginning 

of  the  month  -f  mir 

=  sun's  longitude  at  the  middle  of  the  month. 
Hence  Twill  be  longest  when  o  —  ?=0, 
and  shortest  when  O  -  f  =  TT  ; 

or  T  will  be  longest  when  the  sun  at  the  middle  of  the  month 
is  in  perigee,  and  shortest  when  in  apogee  ;  but,  at  present, 
the  sun  is  in  perigee  about  the  1st  of  January,  and  in  apogee 
about  the  1st  of  July;  therefore,  owing  to  annual  equation, 
the  winter  months  will  be  longer  than  the  summer  months, 
the  difference  between  a  sidereal  month  in  January  and  July, 
from  this  cause,  being  about  20  minutes. 

97.  All  the  inequalities  or  equations,  which  our  expres- 
sions contain,  have  thus  received  a  physical  interpretation. 
They  were  the  only  ones  known  before  Newton  had  estab- 
lished his  theory,  but  the  necessity  for  such  corrections  was 
fully  recognized,  and  the  values  of  the  coefficients  had  already 
been  pretty  accurately  determined  ;  still,  with  the  exception 
of  the  reduction,  which  is  geometrically  necessary,  they  were 
corrections  empirically  made,  and  it  was  scarcely  to  be  ex- 
pected that  any  but  the  larger  inequalities,  viz.  those  of  the 
first  and  second  orders  which  we  have  here  discussed,  could 
-be  detected  by  observation:  we  find,  -however,  that  three 
others,  have,  since  Newton's  time,  been  indicated  by  obser- 
vation before  theory  had  explained  their  cause.  These  are  — 
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the  secular  acceleration^  discovered  by  Halley ;  an  inequality, 
found  by  Mayer,  in  the  longitude  of  the  moon,  and  of  which 
the  longitude  of  the  ascending  node  is  the  argument;  and 
finally  an  inequality  discovered  by  Burg,  which  has  only  of 
late  years  obtained  a  solution.  For  a  further  account  of 
these,  as  also  of  some  other  inequalities  which  theory  has 
made  known,  see  Appendix,  Arts.  (108)  to  (112). 


CHAPTER    VII. 

APPENDIX. 

In  this  chapter  will  be  found  collected  a  few  propositions 
intimately  connected  with  the  results  or  the  processes  of  the 
Lunar  Theory  as  explained  in  the  previous  pages.  Refer- 
ence has  been  made  to  some  of  them  in  the  course  of  the 
work,  and  the  interest  and  importance  of  the  others  are 
sufficient  to  justify  their  introduction  here. 

98.  The  moon  is  retained  in  her  orbit  by  the  force  of 
gravity,  that  is±  by  the  same  force  which  acts  on  bodies  at 
the  surface  of  the  earth. 

The  proof  of  this  is  merely  a  numerical  verification  ;  the 
data  required  from  observation  are, 

the  space  fallen  through  from  rest  in  1"  by  bodies  at   the 

earth's  surface =  161  feet, 

the  radius  of  the  earth =  4000  miles, 

the  periodic  time  of  the  moon  =  27j  days, 

the  distance  of  the  moon  from  the  earth's  centre=60x4000  miles. 

The  force  of  the  earth's  attraction  oc  —. — ^ .     Therefore,  the 

(dist.)2 

space  fallen  through  in  1"  at  distance  of  moon  by  a  body 
moving  from  rest  under  the  earth's  action  =  -^p  feet 

=  -00447  feet. 
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But  the  moon  in  one  second  describes  an  angle  — -^ ^=0), 

during  which  the  approach  to  the  earth 

=  60  x  4000  x  5280  (vers.  CD)  feet 
60  x  4000  x  5280  x 


(27J/.(24)2.(60)4 
=  -00448  feet. 


feet 


Therefore,  the  space  through  which  the  moon  is  deflected 
in  one  second  from  her  straight  path,  is  just  the  quantity 
through  which  she  would  fall  towards  the  earth,  supposing 
her  to  be  subject  to  the  earth's  attraction;  and  we  may, 
therefore,  conclude  that  she  is  retained  in  her  orbit  by  the 
force  of  gravity. 

When  first  Newton,  in  1666,  attempted  to  verify  this 
result,  he  found  a  difference  between  the  two  values  equal 
to  one-sixth  of  the  less:  the  reason  of  his  failure  was  the 
incorrect  measures  of  the  earth,  which  he  made  use  of  in 
his  computation;  and  it  was  not  till  about  16  years  later 
that  he  was  led  to  the  true  result,  by  using  the  more 
correct  value  of  the  earth's  radius  obtained  by  Picart. 

Prmcipia,  lib.  in.,  prop.  4. 

99.    The  moon's  orbit  is  everywhere  concave  to  the  sun. 

Let  S9  Ej  and  M  be  the  centres  of  the  sun,  earth,  and 
moon.  We  ought  first  to  apply 
to  each  body  forces  equal  and 
opposite  to  those  which  act  on 
the  sun  in  order  to  bring  him  to  rest.  These  forces  are, 
however,  so  small  that  we  may  neglect  them,  and  we  shall 
consider  the  moon  to  be  moving  about  the  sun  fixed,  and  to 
be  disturbed  by  the  earth  alone. 
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The  forces  on  M  are,  therefore,  -am  ........  >  .....  -in 


E 
and     -.,.  ..............  in  ME. 


This  last  must  be  resolved  into  two,  one  in  MS,  the  other 
perpendicular  to  it. 

Therefore,  the  whole  central  force  on  the  moon  in  MS 

m'  E 


and  the  proposition  will  be  proved  if  we  shew  that  this  force 
is  always  positive.     Now,  Art.  (25), 


.  27T.SM* 

periodic  time  of  sun     =  -  —  ^  —  =  --  ^  —  nearly, 

*  * 


and  .....................  m00n 


m  E+M 

therefore  -         =  Tfo  ~          nearly  ; 


therefore 


SM3         9  EM* ' 
m'  SM      E  E 


SM*         9  EM'  EM*         *  EM 

4.1.      f  m'          E 

therefore  -         -  is  positive  : 


but  this  is  the  value  of  the  central  force  corresponding  to 
eos^f=—l,  and  is  therefore  its  least  value.  Hence  the 
central  force  always  tends  to  the  sun,  and  the  path  is  always 
concave. 

At  new  moon  the  force  with  which  the  moon  tends  to  the 
sun  is,  therefore,  greater  than  that  with  which  she  tends  to 
the  earth:  the  earth  being  itself  in  motion  in  the  same 
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direction,  and,  at  that  instant,  with  greater  velocity,  will 
easily  explain  how,  notwithstanding  this,  the  moon  still 
revolves  about  it.* 

Central  and  Tangential  Disturbing  Forces. 

100.  We  have  hitherto  considered  the  effects  of  the 
central  and  tangential  disturbing  forces  in  combination ;  but 
it  will  be  interesting  to  determine  to  which  of  them  the 
several  inequalities  principally  owe  their  existence. 

(1)    To  determine  the  effect  of  the  central  disturbing  force. 
Make  T=0; 

72  T) 

therefore  ^2  +u-  jj-a  =  0, 

du  II  u 

p 

or  substituting  for  j^-r,  from  Art.  (52), 


,f 


—  J/t2  —  Jm2  +  |  m*e  cos  c0l  —  f  m2  cos  (2  —  2m) 
+  |m2e  cos  (2  -  2m  -  c)  0 


u  =  a  i+ 196  we  cos  K2  -  2m  ~  °)  6 

-1£2  cos 2(00-7). 

If  we  compare  this  with  the  value  of  w  found  Art.  (52),  we 
see  that  the  elliptic  inequality  and  the  reduction  are  due 
to  the  central  or  radial  force,  as  also  one-half  of  the  varia- 
tion and  three-tenths  of  the  evection. 

It  would  perhaps  be  proper  to  separate  the  absolute  cen- 
tral force  from  the  central  disturbing  force;  the  terms  due 
to  the  latter  are  those  which  contain  m ;  therefore,  the  elliptic 

*  See  the  Author's  Astronomy,  p.  233. 
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inequality  and  the  reduction  are  the  effects  of  the  former, 
except  that  in  the  elliptic  inequality  the  introduction  of  c,  or 
the  motion  of  the  apse,  is  due  to  the  disturbing  force. 

(2)    To  determine  the  effect  of  the  tangential  disturbing  force. 

Let  the  central  disturbing  force  be  zero  ; 

then   Y2—2  =  y-2  (1  —  f  s'2)  =  a,  neglecting  the  inclination,  and, 

ft   it  it 

substituting  for  7^-3  •  JT>  and  for  2  (  -^  +  u  }  I^r3dd  their 
n  u      au  \au         J  j  h,  u 

values  from  Art.  (52),  the  equation  in  u  becomes 


f  1  -f  e  cos  (c0  -  a)  -1-  %m2  cos  {(2  -  2w)  0  -  2/3} 
wiience  w — a  \  i/^  -^  n     <-»^       i 


We  have  here  the  remaining  half  of  the  variation  and  the 
remaining  seven-tenths  of  the  evection  as  the  effects  of  the 
tangential  disturbance.  Also  c  =  l,  or,  to  the  second  order, 
the  tangential  force  has  no  effect  on  the  motion  of  the  apse. 

101.  The  separate  effects  of  the  central  and  tangential 
forces  in  producing  the  inequalities  in  the  longitude  may  be 
traced   in   a   similar   manner.      We  shall  leave  this  as  an 
exercise  for  the  student.     The  annual  equation  will  be  found 
to  be  due  to  the  central  force.     The  other  inequalities  will 
be  divided  as  in  the  radius  vector;   except  as  regards  the 
variation,   of  which   four-elevenths   are   due  to  the  central 
and  the  remaining  seven-elevenths  to  the  tangential  force. 

To  calculate  the  value  of  g  to  the  third  order. 

102.  We  must  here  make  use  of  the  results  which  the 
approximations  to  the  second  order  have  furnished  ;  but  as 
the  value  of  g  is  determined  by  that  term  of  the  differential 
equation  whose  argument  is  gd  —  7,  we  need  only  consider 
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those  terms  which  by  their  combinations  will  lead  to  it  with- 
out rising  to  a  higher  order  than  the  fourth. 
From  Arts.  (51),  (31),  we  obtain 

s  =  k  smg0l  +  %mk  sin  (2  -  Zm-g]  0^ 

Ps—  S 

--  =  -  f  nfs  (1  +  cos  (2  -  2m)  0J, 


From  these  we  obtain 

T)   _    or 


T    ds 


< 

-7-  +  s  =  0,  to  the  second  order  : 
ctu 

therefore  2  (-=£  +  s^  £m  dO  =  0,  to  the  fourth  order. 

Therefore  the  equation  in  s  becomes,  so  far  as  these  terms 
are  concerned, 


assume  s 

therefore  *(!-/)  =  -  |  m*k  4- 


To  find  the  value  of  c  to  the  third  order. 

103.   Proceeding  in  a  similar  manner  from  the  equations, 
Arts.  (52),  (31),  we  get 


N 
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=  a  ~  i«*      {1  +  3  cos  (2  -  3m)  0J, 


we  obtain 

p 

™-5  =  o  -  Jm*a  {1  -  3e  coac0l...—*gme  cos  (2  -  2m  -  c)  0J 

x  {1  +  3  cos  (2  -201)0,} 
=  a  +  %m*a  {3e  +  ^we}  cosc^, 

T 

_s  =-f  w2  {l-4e  cosc^-4w2  cos(2-2m)  6r*fme  cos(2-2m-c)^| 

x  sin  (2  -  2m).  ^ 
=  -f  m2  {sin  (2  -  2m)  6l  -  2e  sin  (2  -  2m  -  c)  Ol  -  lfme  sin  c0j, 

^  =-a{e  sinc01  +  2m2  sin(2  -2m)  0X+  ^5m6  sin  (2-2m-c)01)  ; 
therefore  y^-g  .  -^  =.  .  .-ff  m*ae  cos  C017 

the  other  term  is  of  the  fifth  order. 
T 


f  T 

Jyg-s  J0=+fm2 


therefore   2  f  ^  +  wj  I«ni  dQ  =  —  *£m3ae  cosc01} 
and  the  equation  in  w  becomes 
^+  M  =  a  {1  +  (f  m2 

=  a  {1  +  (f 
Assume 

therefore  ae  (1  —  c2)  =  (f  m2e  +  *££ 

whence  c  =  1  -  jm2  -  *£f  m\ 
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104.   Hence,  to  the  third  order  of  approximation, 
mean  motion  of  apse  _  1  —  c  _  f  m2  +  Wm8  __  8  + 


mean  motion  of  node      a  —  1       f  w*  —  Aw3  ""  8  —  3m  y 

J  4  o  & 

and  since  m  —  ^  nearly,  we  see  that  the  moon's  apse 
progredes  nearly  twice  as  fast  as  the  node  regredes. 

In  the  case  of  one  of  Jupiter's  satellites,  the  periodic 
time  round  Jupiter  is  only  a  few  of  our  days,  and  the 
periodic  time  of  Jupiter  round  the  sun  is  12  of  our  years, 
therefore  w,  the  ratio  of  these  periods,  is  very  small. 

Hence,  the  apse  of  one  of  Jupiter's  satellites  progredes 
along  Jupiter's  ecliptic,  with  pretty  nearly  the  same  velocity 
as  the  node  regredes,  assuming  these  motions  to  be  due  to 
the  sun's  disturbing  force;  they  are,  however,  principally 
due  to  the  oblateness  of  the  planet. 

Parallactic  Inequality. 

105.  In  carrying  on  the  approximations  to  a  higher  order, 
it  is  found  that  the  expressions  for  the  forces  contain  terms 
whose  argument  is  the  moon's  elongation,  i.e.  the  difference 
of  longitude  of  the  sun  and  moon.  These  terms  appear 
in  P  and  in  T  of  the  fourth  order,  and  since  1  —  m,  the 
coefficient  of  9,  is  near  unity,  the  terms  will  become  of  the 
third  order,  and  therefore  of  considerable  importance  in  the 
values  of  u  and  of  L  We  shall  work  out  these  expressions 
from  the  earliest  steps  ;  and  we  select  these  terms  on  account 
of  the  peculiar  use  which  has  been  made  of  them  to  determine 
the  sun's  parallax,  whence  they  have  received  the  name  of 
parallactic  inequality. 

Let  us  go  back  to  the  expressions  for  the  forces 
Arts.  (23),  (24).  In  the  first  place  we  may  remark  that 
the  terms  we  are  about  to  compute  must  be  independent 
of  s  ;  because  any  term  into  the  composition  of  which  s  enters 
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will  necessarily  have  g  in  the  coefficient  of  the  argument. 
We  may  therefore  at  the  outset  suppose  the  orbit  to  coincide 
with  the  ecliptic.  We  have  then 

8M 2  =  r'2  -  Vr'MG  cos  (0  -  0')  +  M G\ 
1          1 


I  O.      I  VJL~<  V^f  .    ft  fLf.  J-J  ^-^ 

80  --c^  =  ir3|l--7-cos(0-0)4— ;2- 


We  must  now  substitute  these  values  in 
E+M       ,(MG       EG\          ,/ 


ME* 


and  pick  out  the  terms  which  have  for  argument  (0  —  0'). 

,  (MG       EG\ 
In  m  l-^3  +  ggjj 


,  'r2  E-M 

3m1-   --  -cos  0-0'  =  ——    cos  0- 


In    wV-g  -  cos  (0-0')   we   shall  have  the  terms 


^j  +  i_s  COs(0-0')  cos2(0-0')}, 

which  produce 

,  {If  £2  -  ^6^2}  ,  . 

s  0-0  = 


and  finally  m'r  f-^m  -  -s*sj  sin(0-  0')  will  produce 
\ftoz        o^/  / 

,  MGP-EG*  n  o  wiV  -E--W  . 

^    -     ~7r-    -   (I-V)  810(0-0')=!^ 
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Therefore  the  new  terms  will  be 


Since  —  (=  -5  Jo  nearly)  is  of  the  second  order,  these  terms 
are  of  the  fourth  order. 

106.  On  the  principle  of  the  superposition  of  small 
disturbances,  we  may  compute  the  effects  of  these  terms  by 
themselves  as  disturbers  of  the  elliptic  path. 


T 


u  —  a  {1  +e  cos(c#  —  a) 

du  •    f  a       \ 

—  =  —  ae  sm(c0—  a) 


therefore 

YT-  „  —  -.  =  0.  to  the  fourth  order, 
A  IT  dtf 

~jm  +  M  =  a,  to  the  second  order, 


Substituting  in  the  differential  equation  for  w,  we  get 


Assume  u  =  a\  1+  .........  +  -^>    M~  cos{(1~7?2)  ^""^1  +••••    5 
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therefore  A  =     ~  fm—p  =  - 

i  —    i  —  7/t 

therefore 


w=a    1+  ......  ~if^  W~M~  cos{(l-wi)  9—  /3}+...    . 


107.   The  corresponding  term  in  the  value  of  0  will  also 
be  of  the  third  order, 


r 

/ 
j 


T 

T2~3  ^  is  of  the  fourth  order  and  will  not  rise  in  t  : 
lfw 


therefore 


Mayer  was  the  first  who  applied  this  term  to  the  deter- 
mination of  the  sun's  parallax,  by  comparing  the  analytical 
expression  of  the  coefficient  with  its  value  as  deduced  from 

W 
observation.     The  values  of  m  and  of  -p.,  and  therefore  of 

^  —  ^  being  pretty  accurately  known,  —  will  be  determined, 

that  is,  the  ratio  of  the  sun's  parallax  to  that  of  the  moon  ; 
but  the  moon's  parallax  is  well  known  ;  'therefore,  also,  that 
of  the  sun  can  be  calculated.  The  value  so  obtained  by 
Laplace  was  8'6";  but  Hansen,  by  the  discussion  of  a  great 
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number  of  modern  observations,  has  found  8'92",  a  value 
which  agrees  very  closely  with  recent  determinations  by  other 
methods. 

Secular  Acceleration. 

108.  Halley,  about  1693,  found,  by  the  comparison  of 
ancient  and  modern  eclipses,  that  the  moon's  mean  revolution 
is  now  performed  in  a  shorter  time  than  at  the  epoch  of  the 
recorded  Chaldean  and  Babylonian  eclipses.  This  pheno- 
menon, called  the  secular  acceleration  of  the  moon's  mean 
motion,  has  not  yet  been  fully  accounted  for.  For  a  long 
time  it  was  altogether  unexplained,  but  in  1787  Laplace 
gave  what  seemed  to  be  a  satisfactory  and  complete  account 
of  it. 

The  value  of  p,  Art.  (54),  on  which  the  length  of  the 
mean  period  depends,  is  found,  when  the  approximation  is 
carried  to  a  higher  order,  to  contain  the  quantity  e  the 
excentricity  of  the  earth's  orbit.  Now,  this  excentricity  is 
undergoing  a  slow  but  continual  change  from  the  action  of 
the  planets,  and  therefore  p,  as  deduced  from  observations 
made  in  different  centuries,  will  have  different  values. 

The  value  of  p  is  at  present  increasing,  or  the  mean 
motion  is  being  accelerated,  and  it  will  continue  thus  to 
increase  for  a  period  of  immense,  but  not  infinite  duration ; 
for,  as  shewn  by  Lagrange,  the  actions  of  the  planets  on 
the  excentricity  of  the  earth's  orbit  will  be  ultimately 
reversed,  e  will  cease  to  diminish  and  begin  to  increase,  and 
consequently  p  will  begin  to  decrease,  and  the  secular 
acceleration  will,  so  far  as  this  cause  is  concerned,  become  a 
secular  retardation. 

It  is  worthy  of  remark  that  the  action  of  the  planets  on 
the  moon,  thus  transmitted  through  the  earth's  orbit,  is  more 
considerable  than  their  direct  action. 
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109.  In  his  investigation,  Laplace  treated  the  excentricity 
e'  of  the  earth's  orbit  as  a  constant  in  the  differential  equa- 
tions, and  only  considered  its  variation  in  the  final  result. 
This  of  course  greatly  simplified  the  work,  but  ought  to  have 
been  looked  upon  only  as  a  first  approximation ;  yet,  strange 
to  say,  the  result  so  obtained  agreed  almost  perfectly  with 
that  which  had  been  deduced  from  a  comparison  of  ancient 
and   modern   observations.      "  Malheureusement   cette   mer- 
veilleuse  concordance  a  laquelle  Laplace  a  du  peut-6tre  une 
partie  de  1'eclat  de  sa  plus  belle  decouverte,  est  sensiblement 
alte're'e  par  les  approximations  suivantes."* 

Some  years  ago  Professor  Adams  took  up  the  question 
ab  initio  without  any  limitation,  and  introduced  the  variability 
of  e  into  the  differential  equations  themselves.  This  strictly 
correct  process  gave  a  value  of  the  secular  variation  only 
about  half  of  Laplace's  which  had  agreed  so  well  with 
observation. 

The  variation  of  excentricity  does  not  therefore  account  for 
the  whole  observed  change,  and  the  mean  motion  of  the  moon 
is  affected  by  some  other  cause  or  causes  at  present  unknown. 
Perhaps  the  remaining  acceleration  is  only  apparent,  and 
arises  from  a  gradual  retardation  of  the  earth's  diurnal  motion. 
The  effect  would  be  precisely  the  same,  and  such  a  retar- 
dation may  be  due  to  the  action  of  the  tides.  Whatever  be 
the  cause,  the  explanation  of  this  unexplained  phenomenon  is 
a  question  worthy  of  the  mathematician's  most  serious  efforts. 

Inequalities  depending  on  the  Figure  of  the  Earth. 

110.  The  earth,  not  being  a  perfect  sphere,  will  not  at- 
tract as  if  the  whole  of  its  mass  were  collected  at  its  centre : 
hence,  some   correction   must    be   introduced    to   take   into 
account  this  want  of  sphericity,  and  some  relation  must  exist 

*  Pontecoulant,  Sytteme  du  Monde.    Supplement,  1860. 
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between  the  oblateness  and  the  disturbance  it  produces. 
Laplace  in  examining  its  effect  found  that  it  satisfactorily 
explained  the  introduction  of  a  term  in  the  longitude  of  the 
moon,  which  Mayer  had  discovered  by  observation,  and  the 
argument  of  which  is  the  true  longitude  of  the  moon's 
ascending  node. 

By  a  comparison  of  the  observed  and  theoretical  values 
of  the  coefficient  of  this  term,  we  may  determine  the 
oblateness  of  the  earth  with  as  great  accuracy  as  by  actual 
measures  on  the  surface. 

111.  By  pursuing  his  investigations,   with  reference   to 
the  oblateness,  in  the   expression   for  the   moon's  latitude, 
Laplace  found  that  it  would   there   give  rise  to   a  term  in 
which  the  argument  was  the  true  longitude  of  the  moon. 

This  term,  which  was  unsuspected  before,  will  also  serve 
to  determine  the  earth's  oblateness,  and  the  agreement  with 
the  result  of  the  preceding  is  almost  perfect,  giving  the 
compression  -gfa*  which  is  about  a  mean  between  the 
different  values  obtained  by  other  methods. 

Perturbations  due  to  Venus. 

112.  After  the  expression  for  the  moon's  longitude  had 
been  obtained  by  theory,  it  was  found  that  there  was  still 
a  slight  deviation  between  her  calculated  and  observed  places, 
and  Burg,  who  discovered  it  by  a  discussion  of  the  observa- 
tions of  Lahire,  Flamsteed,  Bradley,  and  Maskelyne,  thought 
it  could  be  represented  by  an  inequality  whose  period  would 
be   184  years  and  coefficient   15".     This  was  entirely  con- 
jectural, and  though  several  attempts  were  made,  it  was  not 
accounted  for  by  theory. 

*  Pontecoulant,  Systems  du  Monde,  vol.  iv. 
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About  1848,  Professor  Hansen,  of  Seeberg,  in  Gotha, 
having  commenced  a  revision  of  the  Lunar  Theory,  found 
two  terms,  which  had  hitherto  been  neglected,  due  to  the 
action  of  Venus.  One  of  them  is  direct  and  arises  from 
a  '  remarkable  numerical  relation  between  the  anomalistic 
1  motions  of  the  moon  and  the  sidereal  motions  of  Venus 
'  and  the  earth ;  the  other  is  an  indirect  effect  of  an  inequality 
'of  long  period  in  the  motions  of  Venus  and  the  earth, 
'which  was  discovered  some  years  ago  by  the  Astronomer 
1  Eoyal.'* 

The  periods  of  these  two  inequalities  are  extremely  long, 
one  being  273  and  the  other  239  years,  and  their  coefficients 
are  respectively  27*4"  and  23'2".  '  These  are  considerable 
1  quantities  in  comparison  with  some  of  the  inequalities  already 
'recognised  in  the  moon's  motion,  and,  when  applied,  they 
'  are  found  to  account  for  the  chief,  indeed  the  only 
'remaining,  empirical  portion  of  the  moon's  motion  in 
1  longitude  of  any  consequence ;  so  that  their  discovery  may 
'  be  considered  as  a  practical  completion  of  the  Lunar  Theory, 
'  at  least  for  the  present  astronomical  age,  and  as  establishing 
'the  entire  dominion  of  the  Newtonian  Theory  and  its 
'  analytical  application  over  that  refractory  satellite.'t 

Motion  of  the  Ecliptic. 

113.  We  have  seen,  Art.  (14),  that  our  plane  of  reference 
is  not  a  fixed  plane,  but  its  change  of  position  is  so  slow 
that  we  have  been  able  to  neglect  it,  and  it  is  only  when 
the  approximation  is  carried  to  a  higher  order  that  the 
necessity  arises  for  taking  account  of  its  motion. 

*  Report  to  the  Annual  General  Meeting  of  the  Royal  Astronomical  Society, 
Feb.  11,  1848. 

t  Address  of  Sir  John  Herschel  to  the  Meeting  of  the  Royal  Astronomical 
Society. 
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It  has  been  found  to  have  an  angular  velocity,  about  an 
axis  in  its  own  plane,  of  48"  in  a  century,  and  the  correction 
thus  introduced  produces  in  the  latitude  of  the  moon  a  term 

-  CCO  COS  (6  —  <£), 

where  co  is  the  angular  velocity  of  the  ecliptic,  -  the  angular 

c 

velocity  with  which  the  ascending  node  of  the  moon's  orbit 
recedes  from  the  instantaneous  axis  about  which  the  ecliptic 
rotates,  <£  the  longitude  of  this  axis  at  time  £,  and  6  the 
longitude  of  the  moon  at  the  same  instant. 


Let  YAm  be  the  position  of  the  ecliptic  at  time  £, 

A  the  point  about  which  it  is  turning,  YA  =  <£, 
MN  the  moon's  orbit,  M  the  moon,  and  Mm  a  perpen- 
dicular to  the  ecliptic ;  Tm  =  6 ;  Mm  =  lat.  =  /3. 
i  the  inclination   of    the  orbit,   and  N  the  longi- 
tude of  the  node. 

Let  APN'm  be  the  ecliptic  after  a  time  St. 
Any  point  whose  longitude  is  L  may  be  considered  as 
moving  perpendicularly  to  the  ecliptic  with  a  velocity 
co  sin  (£  —  <£). 

Hence,  the  point  N  will  move  in  the  direction  PN  with 
a  velocity  co  sin  (N—  $).  And  N'  will  move  along  PN'  with 
a  velocity  co  sin  (N—  <f>)  cott ; 

dN 
therefore  -^-  =  co  sin  (N—  <£)  coU". 

Again,  the  point   of  the  ecliptic  90°  in   advance  of  N 
will  move  towards  the  moon's  orbit  with  a  velocity 
co  sin  (90°  +  N-  </>) ; 
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therefore  —  =  —  o>  cos  (N—  </>). 

Now,  cot  ij  G>,  and  —  —  -,  —  -  '  =  —  may  be  considered 


stant  in  integrating  ; 

therefore  BN=  co>  cos  (N—  </>)  cot  t, 

fo'  =  c«  sin  (N-  <£), 
and  if  NM=  i/r,  we  have 

"  ceo  x 

>==  _.  __  —  .  cos  r\r_  6). 


cos  z         sin 

Now,  sin  ft  =  sin  i.  sin  i/r  ; 
cos  /3  .  S^?  =  cos  t  .  sin  ty  .  8^'  +  sin  *'.  cos  ^  . 


but  cosz.sini|r=cos/3sin(^—  JV^)  and  cosi|r  =  cos/3  cos(^—  N)  ; 
therefore  S/3  =  —  co>  cos  (0  —  0). 

The  discovery  of  this  term  is  due  to  Professor  Hansen; 
its  coefficient  is  extremely  small,  about  1/5"  ;  but,  being  of 
a  totally  different  nature  from  those  due  to  successive 
approximations,  it  was  thought  desirable  to  examine  it,  and 
the  above  investigation,  which  was  communicated  to  me  by 
J.  C.  Adams,  Esq.,*  will  be  read  with  interest  on  account 
of  its  elegance. 

With  respect  to  the    foregoing  investigation,  perhaps  the  following  remarks 

dN 
will  not  be   superfluous:  —  The  value  of       -  is  not  the  actual  velocity  of  N, 

but  its  velocity  relatively  to  the  position  of  the  node  as  determined  when  the 
motion  of  the  ecliptic  is  neglected  ;  its  integral  is  therefore  SN  the  change  of 
longitude  due  to  this  motion,  and  in  this  integration  no  constant  is  added,  zero 

being  taken  for  the  mean  value.  The  periodic  forms  of  both  ~^-  and  SN  shew 
that  they  oscillate  about  mean  values,  the  time  of  a  complete  oscillation  being 

*  Now  Lowndean  Professor  of  Astronomy  in  the  University  of  Cambridge. 
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that  required  by  sin  (N-(f>)  and  cos  (N-(j>)  to  go  through  their  cycle.  This 
relative  motion  of  the  node  is  analogous  to  that  of  a  particle  moving  in  a  straight 
line  under  the  action  of  a  force  varying  directly  as  the  distance. 

Similar  remarks  apply  to  the  inclination. 

It  must  also  be  borne  in  mind  that  the  result  obtained  8 (3  is  not  the  difference 
between  the  latitude  referred  to  the  actual  ecliptic  and  that  referred  to  a  fixed 
plane,  but  is  the  difference  between  the  calculated  latitude  referred  to  the  actual 
ecliptic  on  supposition  of  its  being  fixed,  and  the  correct  latitude  referred  to  the 
same  actual  ecliptic  when  its  motion  is  taken  into  account. 

We  may  obtain  an  approximate  value  of  the  coefficient 
ceo  by  substituting  for  it  —  ,  where  a  is  the  number  of  seconds 

through  which  the  ecliptic  is  deflected  in  one  year  =  0'48"> 
and  n  is  the  number  of  years  in  which  the  node  of  the  moon's 

orbit  makes  a  complete  revolution  =  18*6;    for  then,  -  -  is 

the  angle  described  by  the   node   in  one  year;     therefore 

na  .  .. .    -         d  (N-  0) 

—  is  the  ratio  of  &> :  -   —^ — — ' ,  supposing  $  to  remain  con- 

A7T  (it 

stant,  which  is  nearly  the  case  ;  therefore 
na       9'3  x  0-48" 


*  This   affords  the  solution    of   a  problem  proposed  in   the  Senate-House 
in  the  January  Examination  of  1852.    Question  21,  Jan.  22. 
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CHAPTER    VIII. 

HISTORY  OF  THE  LUNAR  PROBLEM  BEFORE  NEWTON. 

114.  The  idea  which  most  probably  suggested  itself  to 
the  minds  of  those  men  who  first  considered  the  motion  of 
the  moon  among  the  stars,  was  that  it  described  a  circle  with 
uniform  velocity  about  the  earth  as  a  centre ;  and  this  first 
rough  result  is  represented  in  our  equations  by  neglecting 
all  small  terms  and  writing  6  =pt^  u  =  a. 

It  must,  however,  have  been  very  soon  perceived  that 
the  actual  motion  is  far  from  being  so  simple,  and  that 
the  moon  moves  with  very  different  velocities  at  different 
times. 

115.  The  earliest  recorded  attempts  to  take  into  account 
the    irregularites    of   the   moon's    motion    were    made    by 
Hipparchus  (140  B.C.).     He  imagined  the  moon  to  move  with 
uniform  velocity  in  a  circle,  of  which  the  earth  occupied, 
not  the  centre,  but  a  point  nearer  to  one  side.     By  a  similar 
hypothesis  he   had  accounted  for  the  irregularities   in  the 
sun's  motion,  and  his  success  in  this  led  him  to  apply  it  also 
to  the  moon. 

It  is  clear  that,  on  this  supposition,  the  moon  would  seem 
to  move  faster  when  nearest  the  earth  or  in  perigee,  and 
slower  when  in  apogee,  than  at  any  other  points  of  her 
orbit,  and  thus  an  apparent  unequal  motion  would  be 
produced. 
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Let  BAM  be  a  circle,    CA  a  radius,  E  a  point  in  A  G 
near    (7;     CB,  ED  two    parallel 
lines   making    an    angle    a    with 
CA. 

Suppose  a  body  M  to  describe 
this  circle  uniformly  with  an 
angular  velocity^,  the  time  being 
reckoned  from  the  instant  when 
the  body  was  at  B,  and  the  longi- 
tude as  seen  from  E  being  reckoned 
from  the  line  ED  ; 
therefore  DEM  =  0,  BCM  =pt, 

AEM  =  6-a.    ACM=pt-a. 

EC 
Now  -:  is  a  small  fraction,  and  if  we  represent  it  by 


EC 
sin  M=   m  sin  AEM 


we  shall  have 


=  e  sin  (6  —  a.) 

—  e  sin(^  +  M-  a), 

ir        e  sin  (pt-  a) 
tanJf=—         v\    M  ]   .  : 
1  —  ecos  (pt  —  a) 

this  would  give  If,  and  then  0  by  the  formula  0=pt  +  M. 

This  was  called  an  excentric,  and  the  value  of  e  was  called 
the  excentricity,  which,  for  the  moon,  Hipparchus  fixed  at 
sin5°l'. 

116.  Another  method  of  considering  the  motion  was  by 
means  of  an  epicycle,  which  led  to  the  same  result. 

A  small  circle  PM,  with  a  radius  equal  to  EC  of  previous 
figure,  has  its  centre  in  the  circumference  of  the  circle  EPD 
(which  has  the  same  radius  as  that  of  the  excentric),  and 
moves  round  E  with  the  uniform  angular  velocity  p,  the 
body  M  being  carried  in  the  circumference  of  the  smaller 
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circle,  the  radius  PM  remaining  parallel  to  itself,  or,  which 

is  the   same    thing,    revolving 

from  the  radius  PE  with   the 

same   angular    velocity   p,    so 

that    the    angle   EPM   equals 

PEA. 

Now,  when  the  angle  AEP 
equals  the  angle  ACM  of  the 
former  figure,  it  is  easily  seen 
that  the  two  triangles  EPM, 
EGM  are  equal,  and  there- 
fore the  distance  EM  and  the  angle  AEM  will  be  the  same 
in  both,  and  the  two  motions  will  be  identical. 

117.   The  value  of  e  being  small,  we  find,  rejecting  ea,  &c  , 

M=es\n  (pt  —  a), 
therefore  6  =pt  +  e  sin  (pt  -  a). 

If  we  reject  terms  of  the  second  order  in  our  expression 
for  the  longitude,  and  make  c  =  1,  we  get,  Art.  (55), 

6=pt  +  Ze  sin  (pt-a), 

which  will  be  identical  with  the  above  if  we  suppose  the 
excentricity  of  the  excentric  to  be  double  that  of  the  elliptic 
orbit. 

Ptolemy  (A.D.  140)  calculated  the  excentricity  of  the 
moon's  orbit,  and  found  for  it  the  same  value  as  Hipparchus, 

viz. 

sin  5°  1'  =  -j^,  nearly. 

The  excentricity  in  the  elliptic  orbit  is,  we  know,  about  ^0, 
and  y1^  and  -^Q  will  pretty  nearly  reconcile  the  two  values  of  6 
given  above.  This  shews  us,  that  for  a  few  revolutions  the 
moon  may  be  considered  as  moving  in  an  excentric,  and  her 
positions  in  longitude  calculated  on  this  supposition  will  be 
correct  to  the  first  order. 
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Her  distances  from  the  earth  will  not  however  agree ;  for 
the  ratio  of  the  calculated  greatest  and  least  distances  would 

be  - — ^p  or  £f ,  while  that  of  the  true  ones  would  be  — 
1  —  ra  *~jk 

orfi. 

It  would,  therefore,  have  required  two  different  excentrics 
to  account  for  the  changes  in  the  moon's  longitude  and  in 
her  radius  vector.  Changes  in  the  latter  could  not,  however, 
be  easily  observed  with  the  rude  instruments  the  ancients 
possessed,  and  it  was  very  long  before  this  inconsistency  was 
detected. 

118.  We  have  said  that  the  moon's  longitude,  calculated 
on  the  hypothesis  of  an  excentric,  will  be  pretty  accurate  for 
^  few  revolutions. 

The  data  requisite  for  this  calculation  are,  the  mean 
angular  motion  of  the  moon,  the  position  of  the  apogee,  and 
the  magnitude  of  the  excentricity. 

But  it  was  known  to  Hipparchus  and  to  the  astronomers  of 
his  time,  that  the  point  of  the  moon's  orbit  where  she  seems 
to  move  slowest,  is  constantly  changing  its  position  among 
the  stars.  Now  this  point  is  the  apogee  of  Hipparchus's  ex- 
centric,  and  he  found  that  he  could  very  conveniently  take 
account  of  this  further  change  by  supposing  the  excentric 
itself  to  have  an  angular  motion  about  the  earth  in  the  same 
direction  as  the  moon,  so  as  to  make  a  complete  revolution 
in  about  nine  years,  or  about  3°  in  each  revolution.* 

This  motion  of  the  apsidal  line  follows  also  from  our 
expression  for  the  longitude,  as  shewn  in  Art.  (75).  It  is 
there,  however,  connected  with  an  ellipse  instead  of  an 

*  On  the  supposition  of  an  epicycle,  this  motion  of  the  apse  could  as  easily 
be  represented  by  supposing  the  radius  which  connects  the  moon  with  the  centre 
of  the  epicycle  to  have  this  uniform  angular  velocity  of  about  3°  in  each  revolu- 
tion, and  also  in  the  same  direction. 

P 
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excentric ;  and  though  the  discovery  that  the  elliptic  is  the 
true  form  of  the  fundamental  orbit  was  not  the  next  in  the 
order  of  time  after  those  of  Hipparchus,  yet,  as  all  the  irre- 
gularities which  were  discovered  in  the  intervening  seventeen 
centuries  are  common  both  to  Hipparchus's  excentric  and  to 
Kepler's  ellipse,  it  will  be  as  well  for  us  to  consider  at  once 
this  new  form  of  the  orbit. 


Elliptic  Form  of  the  Orbit. 

119.  We  need  not  dwell  on  the  steps  which  led  to  this 
great  and  important  discovery.  Kepler,  finding  that  the 
predicted  places  of  the  planet  Mars',  as  given  by  the  circular 
theories  then  in  use,  did  not  always  agree  with  the  computed 
ones,  sought  to  reconcile  these  variances  by  other  combina- 
tions of  circular  orbits,  and  after  a  great  number  of  attempts 
and  failures,  and  eight  years  of  patient  investigation,  he 
found  it  necessary  to  discard  the  excentrics  and  epicycles 
altogether,  and  to  adopt  some  new  supposition.  An  ellipse 
with  the  sun  in  the  focus  was  at  last  his  fortunate  hypothesis? 
which  was  found  to  give  results  in  accordance  with  obser- 
vation; and  this  form  of  the  orbit  was,  with  equal  success, 
afterwards  extended  to  the  moon :  but  the  departures  from 
elliptic  motion,  due  to  the  disturbing  force  of  the  sun,  are, 
in  the  case  of  the  moon,  much  greater  than  the  disturbances 
of  the  planet  Mars  by  the  other  planets. 

In  Keplers's  hypothesis,  then,  the  earth  is  to  be  considered 
as  occupying  the  focus  of  an  ellipse,  in  the  perimeter  of  which 
the  moon  is  moving,  no  longer  with  either  uniform  linear  or 
angular  velocity,  but  in  such  a  manner  that  the  radius  vector 
sweeps  over  equal  areas  in  equal  times.- 

This  agrees  with  our  investigation  of  the  motion  of  two 
bodies,  Art.  (10). 
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Evection. 

120.  The  hypothesis  of  an  excentric,  whose  apse  line  has 
a  progressive  motion,  as  conceived  by  Hipparchus,  served  to 
calculate  with  considerable  accuracy  the  circumstances  or 
eclipses  5  and  observations  of  eclipses,  requiring  no  instru- 
ments, were  then  the  only  ones  which  could  be  made  with 
sufficient  exactness  to  test  the  truth  or  fallacy  of  the  sup- 
position. 

Ptolemy  (A.D.  140)  having  constructed  an  instrument,  by 
means  of  which  the  positions  of  the  moon  could  be  observed 
in  other  parts  of  her  orbit,  found  that  they  sometimes  agreed, 
but  were  more  frequently  at  variance  with  the  calculated 
places ;  the  greatest  amount  of  error  always  taking  place  at 
quadrature  and  vanishing  altogether  at  syzygy. 

What  must,  however,  have  been  a  source  of  great  per- 
plexity to  Ptolemy,  when  he  attempted  to  investigate  the 
law  of  this  new  irregularity,  was  to  find  that  it  did  not 
return  in  every  quadrature, — in  some  quadratures  it  totally 
disappeared,  and  in  others  amounted  to  2°  39',  which  was 
its  maximum  value. 

By  dint  of  careful  comparison  of  observations,  he  found 
that  the  value  of  this  second  inequality  in  quadrature  was 
always  proportional  to  that  of  the  first  in  the  same  place, 
and  was  additive  or  subtractive  according  as  the  first  was 
so ;  and  thus,  when  the  first  inequality  in  quadrature  was  at 
its  maximum  or  5°  1',  the  second  increased  it  to  7°  40',  which 
was  the  case  when  the  apse  line  happened  to  be  in  syzygy 
at  the  same  time.* 


*  It  would  seem  as  if  Hipparchus  had  felt  the  necessity  for  some  further 
modification  of  his  first  hypothesis,  though  he  was  unable  to  determine  it;  for 
there  is  an  observation  made  by  him  on  the  moon  in  the  position  here  specified 
when  the  error  of  his  tables  would  be  greatest ;  and  at  a  time  also  when  she 
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But  if  the  apse  line  was  in  quadrature  at  the  same  time  as 
the  moon,  the  second  inequality  vanished  as  well  as  the  first. 

The  mean  value  of  the  two  inequalities  combined  waa 
therefore  fixed  at  6°  20J'. 

121.  To  represent  this  new  inequality,  which  was  sub- 
sequently called  the  Evection,  Ptolemy  imagined  an  excentric 
in  the  circumference  of  which  the  centre  of  an  epicycle  moved 
while  the  moon  moved  in  the  circumference  of  the  epicycle. 

The  centre  of  the  excentric  and  of  the  epicycle  he  supposed 
in  syzygy  at  the  same  time,  and  both  on  the  same  side  of 
the  earth. 


Thus,  if  E  represent  the  earth, 

S sun, 

M moon, 

c  the  centre  of  the  excentric  EKT  in  syzygy, 
,5,  the  centre  of  the  epicycle,  would  also  be  in 

syzygy. 

Now  conceive  c,  the  centre  of  the  excentric,  to  describe  a 
small  circle  about  E  in  a  retrograde  direction  cc',  while  R^ 
the  centre  of  the  epicycle,  moves  in  the  opposite  direction, 


was  in  the  nonagesimal,  so  that  any  error  of  longitude,  arising  from  her  yet 
uncertain  parallax,  would  be  avoided.  Ptolemy,  who  records  the  observation, 
employs  it  to  calculate  the  evection,  and  obtains  a  result  agreeing  with  that 
of  his  own  observations.  (See  Delambre,  Ast,  Ancienne.) 
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in  such  a  manner  that  each  of  the  angles  S'Ec,  S'ER'  may 
be  equal  to  the  synodical  motion  of  the  moon,  that  is,  her 
mean  angular  motion  from  the  sun ;  SES'  being  the  motion 
of  the  sun  in  the  same  time. 

Now  we  have  seen,  Art.  (116),  that  the  first  inequality 
was  accounted  for  by  supposing  the  epicycle  EM  to  move 
into  the  position  rm,  r  and  R  being  at  the  same  distance 
from  Uj  and  rm  parallel  to  BM*  the  first  inequality  being 
the  angle  rEm.  But  when  the  centre  of  the  epicycle  is  at 
E'j  and  R'M'  is  parallel  to  rm^  the  inequality  becomes  R'EM', 
and  we  have  a  second  correction  or  inequality  mEM'. 

122.  That  this  hypothesis  will  account  for  the  phenomena 
observed  by  Ptolemy,  will  be  readily  understood. 

At  syzygies,  whether  conjunction  or  opposition,  the  centres 
of  the  excentric  and  epicycle  are  in  one  line  with  the  earth 
and  on  the  same  side  of  it ;  the  points  r  and  R'  coincide,  as 
also  m  and  M '.  Hence  mEM'  =  0. 

At  quadratures  (figs.  1  and  2)  c  and  R  are  in  a  straight 
line  on  opposite  sides  of  the  earth,  and  therefore  R  and  r  at 

Fig.  1.  Fig.  2. 


their  furthest    distance.     If,  however,  M1  and  m  be  at  the 
same  time  in  this  line,  or,  in  other  words,  if  the  apse  line 


*  For  simplicity  we  leave  out  of  consideration  the  motion  of  the  apse. 
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be  in  quadratures  (fig.  1),  the  angle  mEM'  will  still  be  zero, 
or  there  will  be  no  error  in  the  longitude.  But,  if  the  apse 
line  is  in  syzygy  (fig.  2),  the  angle  mEM'  attains  its  greatest 
value.* 

Ptolemy,  as  we  have  said,  found  this  greatest  value  to  be 
2°  39',  the  angle  mEr  being  then  5°  1'. 

123.  Copernicus  (A.D.  1543),  having  seen  that  Ptolemy's 
hypothesis  gave  distances  totally  at  variance  with  the  obser- 
vations on  the  changes  of  apparent  diameter,f  made  another 
and  a  simpler  one  which  accounted  equally  well  for  the  in- 
equality in  longtitude,  and  was  at  the  same  time  more  correct 
in  its  representation  of  the  distances. 

Let  E  be  the  earth,  OD  an  epicycle  whose  centre  C  de- 
scribes the  circle  C'  GO"  about  E  with  the  moon's  mean 
angular  velocity. 


Let  CO,  a  radius  of  this  epicycle,  be  parallel  to  the  apse 

*  If  Ptolemy  had  used  the  hypothesis  of  an  excentric  instead  of  an  epicycle 
for  the  first  inequality  of  the  moon,  an  epicycle  would  have  represented  the 
second  inequality  more  simply  than  his  method  did.  Dr.  Whewell's  History  of  the 
Inductive  Sciences,  vol.  I.,  p.  230. 

f  See  Delambre,  Ast.  Moderne,  vol.  I.,  p.  116.  Whewell's  History  of  Inductive 
Sciences,  vol.  I.,  p.  395. 


COPERNICUS'S   HYPOTHESIS.  Ill 

line  EA)  and  about  0  as  centre  let  a  second  small  epicycle 
be  described,  the  radii  GO  and  OM  be  so  taken  that 

CO-  OM  ,  CO+OM      .     .  4 

sins   l,and—  =sm7   40'. 


- 

The  radius  OM  must  now  be  made  to  revolve  from  the 
radius  OC  twice  as  rapidly  as  EG  moves  from  ES,  so  that 
the  angle  COM  may  be  always  double  of  the  angle  CES. 

From  this  construction,  it  follows  that  in  syzygies  the 
angle  CES  being  0°  or  180°,  the  angle  COM  is  0°  or  360°; 
and  therefore  C  and  M  are  at  their  nearest  distances,  as  in 
the  positions  C1  and  C'"  in  the  figure.  Then  CM=  CO  -  OM, 
and  the  angle  GEM  will  range  between  0°  and  5°  1',  the 
greatest  value  being  attained  when  the  apse  line  is  in  quad- 
rature. 

When  the  moon  is  in  quadrature  CES=$Q°  or  270°, 
and,  therefore,  CGM=1SO°  or  540°  and^  C  and  M  are  at 
their  greatest  distance  apart,  as  in  the  position  C"  ';  then, 
CM  =  CO  +  OM,  and  the  angle  GEM  will  range  between 
0°  and  7°  40',  the  former  value  when  the  apse  line  is  itself 
in  quadrature,  and  the  latter  when  it  is  in  syzygy. 

124.  Thus  the  results  attained  by  Ptolemy's  construction 
are,  as  far  as  the  longitudes  at  syzygies  and  quadratures  are 
concerned,  as  well  represented  by  that  of  Copernicus;  and 
the  variations  in  the  distances  of  the  moon  will  be  far  more 
exact,  the  least  apparent  diameter  being  28'  45"  and  the 
greatest  37'  33"  ;  whereas,  Ptolemy's  would  make  the  greatest 
diameter  1°.* 

The  values  which  modern  observations  give  vary  between 
28'  48"  and  33'  32". 


Delambre,  A st,  Moderne. 
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125.  It  will  not  now  be  difficult  to  shew  that  the  intro- 
duction of  this  small  epicycle  corresponds  with  that  of  the 
term  l-fme  sin  {(2  —  2m  -  c)  pt  —  2/8  +  a}  in  our  value  of  6. 

For,  referring  to  the  preceding  figure,  we  have 

OEM=  sin  OEM  =  ~  sin  OME 

UJ.ii 

sin  (  COM-  AEM) 


OM  .    .    . 
==TTPsm  '    (moon  s  wean  long.  -  sun  s  long.) 

—  (moon's  true  long.  -  long,  of  apse)}? 

and  OEM  being  a  small  angle  whose  maximum  is  1°  19J', 
we  may  write  moon's  mean  longitude  instead  of  the  true 
in  the  argument,  and  also  EC  for  OE  ;  therefore, 

OEM=  -=^  sin  {2  (moon's  mean  longitude  -  sun's  longitude) 

—  (moon's  mean  longitude  —  longitude  of  apse)} 
=  79  y  sin  [2  {pt-  (mpt  +  j3)}  -  {pt-(l-c)pt  +  a}'] 
=  4770"  sin  j(2  -  2m  -  c}pt  -  2/3  +  a}. 

The  value  of  the  coefficient  is  from  modern  observations 
found  to  be  4589'61". 

126.  In  Art.  (77),  we  have  considered  the  effect  of  this 
second  inequality  in  another  light,  not  simply  as  a  small 
quantity  additional  to  the  first  or  elliptic  inequality,  but  as 
forming  a  part  of  this  first;  and  therefore,  modifying  and 
constantly  altering  the  excentricity  and  the  uniform  pro- 
gression of  the  apse  line. 

Boulliaud  (A.  D.  1645),  by  whom  the  term  Evection  was 
first  applied  to  the  second  inequality,  seems  to  hint  at  some- 
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thing  of  this  kind  in  the  rather  obscure  explanations  of  Ms 
lunar  hypothesis,  which,  never  having  been  accepted,  it  would 
be  useless  to  give  an  account  of  it.* 

In  Ptolemy's  theory,  Art.  (121),  the  evection  was  the 
result  of  an  apparent  increase  of  the  first  lunar  epicycle 
caused  by  its  approaching  the  earth  at  quadratures ;  but,  in 
this  second  method,  it  is  the  result  of  an  actual  change  in 
the  elements  of  the  elliptic  orbit. 

D'Arzachel,  an  Arabian  astronomer,  who  observed  in 
Spain  about  the  year  1080,  seems  to  have  discovered  the 
unequal  motion  of  the  apsides,  but  his  discovery  must  have 
been  lost  sight  of,  for  Horrocks,  about  1640,  re-discovered  it 
4  in  consequence  of  his  attentive  observations  of  the  lunar 
1  diameter :  he  found  that  when  the  distance  of  the  sun  from 
4  the  moon's  apogee  was  about  45°  or  225°,  the  apogee  was 
{  more  advanced  by  25°  than  when  that  distance  was  about 
1 135°  or  315°.  The  apsides,  therefore,  of  the  moon's  orbit 
1  were  sometimes  progressive  and  sometimes  regressive,  and 
'required  an  equation  of  12°  30',  sometimes  additive  to  their 
1  mean  place  and  sometimes  subtractive  from  it.'f 

Horrocks  also  made  the  excentricity  variable  between  the 
limits  -06686  and  -04362, 

The  combination  of  these  two  suppositions  was  a  means 
of  avoiding  the  introduction  of  Ptolemy's  excentric  or  the 
second  epicycle  of  Copernicus :  their  joint  effect  constitutes 
the  evection. 


*  Apres  avoir  etabli  les  mouvements  et  les  epoques  de  la  lune,  Boulliaud  revient 
a  1'explication  de  1'evection  ou  de  la  secoude  inegalite.  Si  sa  theorie  n'a  pas  fait 
fortune,  le  nom  du  moins  est  reste.  '  En  meme  temps  que  la  lune  avance  sur  son 
'  cone  autour  de  la  terre,  tout  le  systeme  de  la  lune  est  deplace ;  la  terre  emportant  la 
'lune,  rejette  loin  d'elle  1'apogee,  et  rapproche  d'autant  le  perigee  j  mais  cette 
'  evection  a  des  bornes  fixees.'  Delambre,  Hist,  de  I'  Ast.  Mod.,  torn.  II.  p.  157. 

f  Small's  Astronomical  Discoveries  of  Kepler,  p.  307. 

Q 
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Variation. 

]  27.  After  the  discovery  of  the  evection  by  Ptolemy,  a 
period  of  fourteen  centuries  elapsed  before  any  further  ad- 
dition was  made  to  our  knowledge  of  the  moon's  motions. 
Hipparchus's  hypothesis  was  found  sufficient  for  eclipses,  and 
when  corrected  by  Ptolemy's  discovery,  the  agreement  be- 
tween the  calculated  and  observed  places  was  found  to 
extend  also  to  quadratures;  any  slight  discrepancy  being 
attributed  to  errors  of  observation  or  to  the  imperfection 
of  instruments. 

But  when  Tycho  Brahe  (A.  D.  1580)  with  superior  instru- 
ments extended  the  range  of  his  observations  to  all  inter- 
mediate points,  he  found  that  another  inequality  manifested 
itself.  Having  computed  the  places  of  the  moon  for  differ- 
ent parts  of  her  orbit  and  compared  them  with  observation, 
he  perceived  that  she  was  always  in  advance  of  her  com- 
puted place  from  syzygy  to  quadrature,  and  behind  it  from 
quadrature  to  syzygy ;  the  maximum  of  this  variation  taking 
place  in  the  octants,  that  is,  in  the  points  equally  distant 
from  syzygy  and  quadrature.  The  moon's  velocity  therefore, 
so  far  as  this  inequality  was  concerned,  was  greatest  at  new 
and  full  moon,  and  least  at  the  first  and  third  quarters.* 


*  'It  appears  that  Mohammed- Aboul-Wefa-al-Bouzdjani,  an  Arabian  astro- 
nomer of  the  tenth  century,  who  resided  at  Cairo,  and  observed  at  Bagdad  in  975, 
'  discovered  a  third  inequality  of  the  moon,  in  addition  to  the  two  expounded  by 
'  Ptolemy,  the  equation  of  the  centre  and  the  evection.  This  third  inequality,  the 
'variation,  is  usually  supposed  to  have  been  discovered  by  Tycho  Brahe,  six 

'centures  later In  an  almagest  of  Aboul-Wefa,   a  part  of  which   exists  in 

1  the  Royal  Library  at  Paris,  after  describing  the  two  inequalities  of  the  moon, 
'he  has  a  Section  IX.,  "Of  the  third  anomaly  of  the  moon  called  Muhazal  or 

'  Prosneusis " But   this   discovery  of  Aboul-Wefa  appears   to  have   excited 

'no  notice  among  his  contemporaries  and  followers;  at  least  it  had  been  long 
*  quite  forgotten,  when  Tycho  Brahe  re-discovered  the  same  lunar  inequality.' 
Whewell's  Hist,  of  Inductive  Sciences,  vol.  I.  p.  243. 
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Tycho  fixed  the  maximum  of  this  inequality  at  40'  30". 
The  value  which  results  from  modern  observations  is 
39'  30". 

128.  We  have  already  two  epicycles,  or  one  epicycle  and 
an  excentric,  to  explain  the  first  two  inequalities:    by  the 
introduction  of  another  epicycle  or  excentric,  the  variation 
also  might  have  been  brought  into  the  system ;  but  Tycho 
adopted   a  different  method:*    like  Ptolemy,  he  employed 
an  excentric  for  the  evection,  but  for  the  first  or  elliptic 
inequality    be    employed    a    couple    of   epicycles,   and   this 
complicated    combination,   which   it   is   needless    further   to 
describe,   represented   the   change   of    distance   better  than 
Ptolemy's. 

To  introduce  the  variation,  he  imagined  the  centre  of  the 
larger  epicycle  to  librate  backwards  and  forwards  on  the  ex- 
centric,  to  an  extent  of  40^'  on  each  side  of  its  mean  position ; 
this  mean  place  itself  advancing  uniformly  along  the  excentric 
with  the  moon's  mean  motion  in  anomaly ;  and  the  libration 
was  so  adjusted,  that  the  moon  was  in  her  mean  place  at 
syzygy  and  quadrature,  and  at  her  furthest  distance  from  it 
in  the  octants,  the  period  of  a  complete  libration  being  half 
a  synodical  revolution. 

Annual  Equation. 

129.  Tycho  Brahe*  was  also  the  discoverer  of  the  fourth 
inequality,  called  the  annual  equation.     This  was  connected 
with  the  anomalistic  motion  of  the  sun,  and  did  not,  like  the 
previous  inequalities,  depend  on  the  position  of  the  moon  in 
her  orbit. 


*  For  a  full  description  of  Tycho's  hypothesis,  see  Delambre,  Hist,  de  H  Ast. 
Mod.  torn.  I.  p.  162,  and  An  Account  of  the  Astronomical  Discoveries  of  Kepler 
by  Robert  Small,  p.  135. 
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Having  calculated  the  position  of  the  moon  corresponding 
to  any  given  time,  he  found  that  the  observed  place  was 
behind  her  computed  one  while  the  sun  moved  from  perigee 
to  apogee,  and  before  it  in  the  other  half  year. 

Tycho  did  not  state  this  distinctly,  but  he  made  a  cor- 
rection which,  though  wrong  in  quantity  and  applied  in  an 
indirect  manner,  shewed  that  he  had  seen  the  necessity  and 
understood  the  law  of  this  inequality. 

He  did  not  try  to  represent  it  by  any  new  excentric  or 
epicycle,  but  he  increased  by  (8m.  13s.)  sin(sim's  anomaly] 
the  time  which  had  served  to  calculate  the  moon's  place;* 
thus  assuming  that  the  true  place,  after  that  interval,  would 
agree  with  the  calculated  one.  Now,  as  the  moon  moves 
through  4'  30"  in  8m.  13s.,  it  is  clear  that  adding  (8m.  13s.) 
sin  (sun's  anomaly]  to  the  time  is  the  same  thing  as  sub- 
tracting (4'  30")  sin  (suns  anomaly]  from  the  calculated  lon- 
gitude, which  was  therefore  the  correction  virtually  intro- 
duced by  Tycho.f  Modern  observations  shew  the  coefficient 
to  be  11'  9". 

We  have  seen,  Art.  (82),  how  this  inequality  may  be 
inferred  from  our  equations. 

Reduction. 

130.  The  next  inequality  in  longitude  which  we  have  to 
consider  is  not  an  inequality  in  the  same  sense  as  the  fore- 
going ;  that  is,  it  does  not  arise  from  any  irregularity  in  the 
motion  of  the  moon  herself  in  her  orbit,  but  simply  because 
that  orbit  is  not  in  the  same  plane  as  that  in  which  the 


*  That  is,  the  equation  of  time  which  he  used  for  the  moon  differed  by  that 
quantity  from  that  used  for  the  sun. 

t  Horrocks  (1639)  made  the  correction  in  the  same  manner  as  Tycho,  but  so 
increased  it  that  the  corresponding  coefficient  was  11'  51"  instead  of  4'  30". 
Flamsteed  was  the  first  to  apply  the  correction  to  the  longtitude  instead  of  the  time. 
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longitudes  are  reckoned,  so  that  even  a  regular  motion  in 
the  one  would  be  necessarily  irregular  when  referred  to  the 
other.  Thus  if  NMn  be  the  moon's  orbit  and  TJVm  the 
ecliptic,  and  if  M  the  moon  be  re- 
ferred to  the  ecliptic  by  the  great 
circle  Mm  perpendicular  to  it,  then 
MN  and  mN  are  0°,  90°,  180°,  270°, 
and  360°  simultaneously,  but  they 
differ  for  all  intermediate  values:  the  difference  between 
them  is  called  the  reduction. 

The  difference  between  the  longitude  of  the  node  and  that 
of  the  moon  in  her  orbit  being  known,  that  is  the  side  NM  of 
the  right-angled  spherical  triangle  NMm,  and  also  the  angle 
N  the  inclination  of  the  two  orbits,  the  side  Nm  may  be 
calculated  by  the  rules  of  spherical  trigonometry,  and  the 
difference  between  it  and  NM^  applied  with  a  proper  sign 
to  the  longitude  in  the  orbit,  gives  the  longitude  in  the 
ecliptic. 

Tycho  was  the  first  to  make  a  table  of  the  reduction 
instead  of  calculating  the  spherical  triangle.  His  formula 
was 

reduction  =  tan2^/  sin2L  —  \  tan4^/  sin4L, 

where  /  is  the  inclination  of  the  orbit  and  L  the  longitude  of 
the  moon  diminished  by  that  of  the  node. 

The  first  term  corresponds  with  the  term  —  %k2  sin  2  (gpt  —  7) 
of  the  expression  for  6. 

Latitude  of  the  Moon. 

131.  That  the  moon's  orbit  is  inclined  to  the  ecliptic  wa& 
known  to  the  earliest  astronomers,  from  the  non-recurrence 
of  eclipses  at  every  new  and  full  moon;  and  it  was  also 
known,  sinces  the  eclipses  did  not  always  take  place  in  the 
same  parts  of  the  heavens,  that  the  line  of  nodes  represented 
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by  Nn,  in  the  preceding  figure,  has  a  retrograde  motion  on 
the  ecliptic,  N  moving  towards  T. 

Hipparchus  fixed  the  inclination  of  the  moon's  orbit  to 
the  ecliptic  at  5°,  which  value  he  obtained  by  observing  the 
greatest  distance  at  which  she  passes  to  the  north  or  south 
of  some  star  known  to  be  in  or  very  near  the  ecliptic,  as  for 
instance  the  bright  star  Regulus ;  and  by  comparing  the  re- 
corded eclipses  from  the  times  of  the  Chaldean  astronomers 
down  to  his  own,  he  found  that  the  line  of  nodes  goes  round 
the  ecliptic  in  a  retrograde  direction  in  about  18J  years. 

This  result  is  indicated  in  our  expression  for  the  value 
of  the  latitude  by  the  term  k  sin  (g6  —  7),  as  we  have  shewn 
Art.  (85). 

132.  Tycho  Brane*  further  discovered  that  the  inclination 
of  the  lunar  orbit  to  the  ecliptic  was  not  a  constant  quantity 
of  5°  as  Hipparchus  had  supposed,  but  that  it  had  a  mean 
value  of  5°  8',  and  ranged  through  9'  30"  on  each  side  of 
this,  the  least  inclination  4°  .58J'  occurring  when  the  node 
was  in  quadrature,  and  the  greatest  5°  17J'  being  attained 
when  the  node  was  in  syzygy.* 

He  also  found  that  the  retrograde  motion  of  the  node  was 
not  uniform:  the  mean  and  the  true  positions  agreed  very 

*  Ebn  Jounis,  an  Arabian  astronomer  (died  A.D.  1008),  whose  works  were 
translated  about  50  years  since  by  Mons.  Sedillot,  states  that  the  inclination  of  >the 
moon's  orbit  had  been  often  observed  by  Aboul-Hassan-Aly-ben-Amajour  about 
the  year  918,  and  that  the  results  he  had  obtained  were  generally  greater  than 
the  5°  of  Hipparchus,  but  that  they  varied  considerably. 

Ebn  Jounis  adds,  however,  that  he  himself  had  observed  the  inclination 
several  times  and  found  it  5°  3',  which  leads  us  to  infer  that  he  always  observed 
in  similar  circumstances,  for  otherwise  a  variation  of  nearly  23'  could  scarcely 
have  escaped  him.  See  Delambre,  Hist,  de  V  Ast.  du  Moyen  Age,  p.  139. 

The  mean  value  of  the  inclination  is  5°  8'  55:46",— the  extreme  values  are 
4°  57'  22"  and  5°  20'  6".  4 

The  mean  daily  motion  of  the  line  of  nodes  is  3'  10  '64",  or  one  revolution  in 
6793-29  days,  or  18  y.  218  d.  21  h.  22m.  46  s. 
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well  when  they  were  in  syzygy  or  quadrature,  but  they 
were  1°  46'  apart  in  the  octants. 

By  referring  to  Art.  (87)  we  shall  see  that  these  correc- 
tions, introduced  by  Tycho  Brahe,  correspond  to  the  second 
term  of  our  expression  for  s. 

Since  Hipparchus  could  observe  the  moon  with  accuracy 
only  in  the  eclipses,  at  which  time  the  node  is  in  or  near 
syzygy,  we  see  why  he  was  unable  to  detect  the  want  of 
uniformity  in  the  motion  of  the  node. 

133.  To  represent  these  changes  in  the  position  of  the 
moon's  orbit,  Tycho  made  the  following  hypothesis. 

Let  ENF  be  the  ecliptic,  K  its  pole,  BA  C  a  small  circle, 


having  also  K  for  pole  and  at  a  distance  from  it  equal 
to  5°  8'.  Then,  if  we  suppose  A  the  pole  of  the  moon's 
orbit  to  move  uniformly  in  the  small  circle  and  in  the 
direction  BAG,  the  node  N,  which  is  at  90°  from  both  A 
and  K,  will  retrograde  uniformly  on  the  ecliptic,  and  the 
inclination  of  the  two  orbits  will  be  constant  and  equal  to  AK. 
But  instead  of  supposing  the  pole  of  the  moon's  orbit  to 
be  at  Ay  let  a  small  circle  abed  be  described  with  A  as  pole 
and  a  radius  of  9'  30" ;  and  suppose  the  pole  of  the  moon's 
orbit  to  describe  this  small  circle  with  double  the  velocity 
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of  the  node  in  ,  its  synodical  revolution  which  is  accomplished 
in  about  346  days,  in  such  a  manner  that  when  the  node  is 
in  quadrature  the  pole  may  be  at  a,  the  nearest  point  to  7T, 
and  at  c  the  most  distant  point  when  the  node  comes  to 
syzygy,  at  d  in  the  first  and  third  octants,  and  at  b  in  the 
second  and  fourth,  so  as  to  describe  the  small  circle  in  about 
173  days,  the  centre  A  of  the  small  circle  retrograding  mean- 
while with  its  uniform  motion. 

By  this  method  of  representing  the  motion,  we  see  that 

when  node  is  in  quadrature) 


while  at  the  octants  it  has  its  mean  value 


Again,  with  respect  to  the  motion  of  the  node,  since  N  is 
the  pole  of  KaAc^  it  follows  that  when  in  syzygy  and  quad- 
rature, the  node  occupies  its  mean  place;  in  the  first  and 
third  octants,  the  pole  being  at  d,  the  node  is  before  its 
mean  place  by  the  angle  dKA  =  (9'  30")  cosec5°  8'  =  1°  46', 
nearly,  and  it  is  as  much  behind  in  the  second  and  fourth 
octants. 

So  that  the  whole  motion  of  the  node,  and  the  correction 
which  Tycho  had  discovered,  were  properly  represented  by  this 
hypothesis,  which  is  exactly  similar  to  that  which  Copernicus 
had  imagined  to  explain  the  precession  of  the  Equinox. 
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SELECTION  OF  EXAMINATION  QUESTIONS 

FROM    COLLEGE  AND  SENATE-HOUSE   PAPERS   AND  FROM   THE 
MODERATORSHIP  AND  FELLOWSHIP  EXAMINATIONS  AT 
TRINITY  COLLEGE,  DUBLIN. 

1.  Define  the  plane  of  the  ecliptic  and  prove  that,  as 
seen  from  the  earth,  sin  ^  =  -g-^J  o"o"  sinV  nearly  5  A,  being 
the  latitude  of  the  sun  and  V  that  of  the  moon. 

2.  Obtain  the  differential  equation  of  the  moon's  radius 
vector. 

3.  What  is  the  principle  on  which  the  successive  approxi- 
mations to  the  moon's  path  are  obtained  ? 

Compare  the  coefficients  of  the  principal  term  in  the  cen- 
tral disturbing  force  with  the  principal  term  in  the  central  force. 

4.  If  in  solving  the  equation 

fu  P         T  du         tfu 


we  tiake  as  our  first  approximation  u  =  a  {1+e  cos(0-a)}, 
we  obtain  at  our  higher  approximations  terms  containing  6  in 
their  coefficient.  How  is  this  defect  avoided  by  taking 
u  =  a  }1  +  e  cos  (cO  —  a)}  as  our  first  approximation  ? 

Upon  what  principle  do  we  approximate  to  the  value  of  c  ? 

T 
5.   Shew  that  in  the  expansion  of  73—3  there  will  be  a 

/&.% 

term  -  ^mV  sin  {(2  -  2m  -  2o)  6  -  20  +  2oc}  which  would 
rise  to  the  second  order  in  the  longitude. 

Shew  further  that  the  term  with  the  same  argument  does 
not  appear  actually  in  the  longitude  —  assuming 

u=a  [l  +  ecos(c0-  a)+...+  */mecos{(2-  2m-  c)  0-2/3  +  a] 
+...+  yme1  cos  {(2  -  2m  -  2c)  d  -  2/3  +  2a}+...]. 

R 
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6.  The  longitude  of  the  moon  contains  a  term  of  the  form 

A  sin  {(2  -  2m  -  c)  pt-  2/3  -f  a] 

where   A   is   of  the   second  order;  find  how  it  must  have 
arisen  and  determine  A. 

7.  Investigate   the  following  expression  for  the  moon's 
longitude  as  far  as  the  second  order,  supposing  the  orbits 
of  the  moon  relative  to  the  earth,  and  of  the  earth  relative 
to  the  sun,  to  be  originally  circles  in  the  plane  of  the  ecliptic, 


8.  Suppose  that  whatever  constants  are  involved  in  the 
arguments  0,  $,  >/r,  &c.  in  the  equation 

L  =  A  +  B  sin  6  +  C  sin  $  +  D  sin  ty  +  &c. 
are  perfectly  known,  and  that  observation  gives  the  numerical 
values  of  a  considerable  number  of  longitudes  L^  L^  L^  &c., 
corresponding  to  the  known  angles  fl,,^,-^,,  &c.,  02,02,^o,  &c., 
03,  <£s,  o/rg,  &c.,  &c.,  shew  how  the  constants  A,  B,  u,  &c. 
may  be  numerically  determined  in  the  two  distinct  cases 
where  the  observations  at  our  disposal  are  and  are  not 
unlimited  in  number. 

9.  Explain  how  this  process  fails  and  may  be  modified 
when  two  of  the  terms  in  question  B  sin  0,  C  sin  <£  are  nearly 
synchronous  in  their  periods. 

10.  Shew  that  the  evection  in  longitude,  viz. 

i-fme  sin  {(2  -  2m  -  c)  pt  -  20  +  a), 

may  be  represented  as  the  joint  effect  of  certain  periodic 
changes  in  the  excentricity  of  the  lunar  orbit  and  in  the 
mean  longitude  of  its  apse. 

11.  Assuming  the  usual  notation  of  the  Lunar  Theory, 
explain  the  physical  meaning  of  the  following  equations  : 

s  =  k  sin(0  —  7), 


s  =  k  sin  (#0  -  7)  +  f  w&  sin  {(2  -  2m  -g)  0-2/3  +  7}. 
1  2.   Explain  the  effect  of  the  term  mza  cos  {(2-2/n)  0-2/9)} 
in  the  moon's  radius  vector;   find  the  number  of  days  in 
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the  period  of  the  resulting  inequality  and  the  ratio  of  the 
axes  of  the  oval  orbit. 

13.  Investigate  the  effect  of  the  annual  equation  on  the 
length  of  the  lunar  month. 

Assuming  w='075,  e='017,  sin  13°30'='23345,  and  the  mean 
period  of  a  sidereal  revolution  of  the  moon  =  27d.  7h.  43m.  12s., 
find  the  difference  between  a  winter  and  a  summer  month. 

14.  Shew  that  the  moon's  orbit  in  space  is  always  concave 
towards  the  sun. 

Assuming  that  the  earth  describes  a  circle  about  the  sun, 
and  the  moon  a  circle  about  the  earth  and  in  the  same  plane 
with  the  earth's  orbit,  compare  the  curvatures  of  the  moon's 
orbit  in  space  at  her  perihelion  and  aphelion. 

15.  Prove  that  if  we  go  to  the  third  order  of  approxi- 
mation, the  motion   of  the  moon's  apse  in  one  revolution 
of  the  moon  equals  fra2  (1  -  7fm)  360.° 

P  T 

16.  The  expressions  for  y^  and  ^3-5  contain  terms  of 

AV  h*u3 

the  form  ZA  cos  (9  —  6'}  and  A  sin  (6  —  6')  respectively ;  com- 
pute the  terms  resulting  from  them  in  the  values  of  u  and  6. 

17.  Copernicus  represented  the  evection  by  an  epicycle 
superposed   on  the  epicycle   which   represented   the   elliptic 
inequality.      From    the   construction    given    by   him    shew 
geometrically  that  the  evection  vanishes  when  the  moon's 
mean  elongation  from   the  sun  is   half  her  true   anomaly, 
and,   in   general,  varies  nearly    as  the   sine   of   twice   the 
difference  of  these  angles. 

18.  The  secular  change  of  inclination  of  the  actual  ecliptic 
to  a  fixed  ecliptic  being  a"  annually,  shew  that  this  will 
give  rise  to  an  inequality  in  the  moon's  latitude  whose  type 

71 

is  —  —  a"cos<£,  nearly;    n  being  the  number  of   years  in 

which  the  moon's  nodes  make  a  complete  revolution,  and 
<j>  the  difference  of  longitude  of  the  moon  and  of  the  ascending 
node  of  the  actual  ecliptic. 
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